Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 71,175 A→G 34.8% R263R (CGA→CGG araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP00904871,1750AG34.8% 71.6 / 67.2 69R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (24/21);  new base G (11/13);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATG  >  AP009048/71106‑71245
                                                                     |                                                                      
cACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                       >  1:1853943/1‑71 (MQ=255)
 aCTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGt                                                                      <  1:1168361/71‑1 (MQ=255)
    aCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                   >  1:1901644/1‑71 (MQ=255)
     cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                       >  1:2291890/1‑66 (MQ=255)
     cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                   >  1:3768573/1‑70 (MQ=255)
     cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAt                                                                    <  1:4159497/69‑1 (MQ=255)
      ccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                   >  1:2682779/1‑69 (MQ=255)
      ccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTa                                                                 <  1:3288642/71‑1 (MQ=255)
       cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTa                                                                 >  1:3476986/1‑70 (MQ=255)
         gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:3428484/71‑1 (MQ=255)
         gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:4242040/71‑1 (MQ=255)
         gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaa                                                               <  1:259290/70‑1 (MQ=255)
          aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              >  1:1787021/1‑70 (MQ=255)
           tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            >  1:3394762/1‑71 (MQ=255)
           tGCCTATCGCCACCGTCGGTCGAAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:3063069/71‑1 (MQ=255)
             ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGa                                                                        <  1:2600136/57‑1 (MQ=255)
             ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAAt                                                           <  1:2027107/70‑1 (MQ=255)
                aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            >  1:491247/1‑66 (MQ=255)
                aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                       >  1:3755881/1‑71 (MQ=255)
                aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAc                                                       >  1:4533673/1‑71 (MQ=255)
                 tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                      >  1:4469742/1‑71 (MQ=255)
                 tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                      >  1:2913670/1‑71 (MQ=255)
                   gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                     <  1:1542018/70‑1 (MQ=255)
                   gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCgg                                                    <  1:2139828/71‑1 (MQ=255)
                   gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCg                                                     <  1:3684374/70‑1 (MQ=255)
                      aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaa                                                                >  1:3966632/1‑56 (MQ=255)
                       ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                            <  1:95880/59‑1 (MQ=255)
                        cGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCa                                               <  1:4768938/71‑1 (MQ=255)
                         gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAg                                              >  1:3764765/1‑71 (MQ=255)
                          tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                       >  1:628707/1‑45 (MQ=255)
                          tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaa                                                             >  1:2891430/1‑55 (MQ=255)
                          tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAAt                                                           <  1:2407380/57‑1 (MQ=255)
                            ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGccc                                           <  1:3202595/71‑1 (MQ=255)
                            ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                            >  1:3368076/1‑54 (MQ=255)
                              tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGccc                                           >  1:3738755/1‑69 (MQ=255)
                                gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                       <  1:3707175/71‑1 (MQ=255)
                                gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                       <  1:906083/71‑1 (MQ=255)
                                 cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                      <  1:3300268/71‑1 (MQ=255)
                                  aaTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                       >  1:2993005/1‑69 (MQ=255)
                                   aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                    >  1:1492410/1‑71 (MQ=255)
                                    tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGt                                   >  1:1991085/1‑71 (MQ=255)
                                      ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                                  <  1:4991623/38‑1 (MQ=255)
                                      ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                                  <  1:2900364/70‑1 (MQ=255)
                                         gTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCg                               <  1:4748361/70‑1 (MQ=255)
                                          ttttGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                             >  1:4696798/1‑39 (MQ=255)
                                           tttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                         <  1:614555/42‑1 (MQ=255)
                                            ttGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:1806847/38‑1 (MQ=255)
                                            ttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCggg                                                   >  1:4452126/1‑47 (MQ=255)
                                            ttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGc                            <  1:5128710/70‑1 (MQ=255)
                                            ttGAAGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGc                            <  1:2964427/70‑1 (MQ=255)
                                              gACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                          >  1:4327307/1‑38 (MQ=255)
                                              gACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                          >  1:5054246/1‑38 (MQ=255)
                                              gACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGc                                                 >  1:4840205/1‑47 (MQ=255)
                                                  aTCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGcc                                                <  1:1609334/44‑1 (MQ=255)
                                                   tCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                     <  1:4354683/38‑1 (MQ=255)
                                                   tCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgt                     >  1:317514/1‑70 (MQ=255)
                                                   tCAACTCTATTTCTCGCGAGTATTTAAAAAATGCa                                                        >  1:837961/1‑35 (MQ=255)
                                                   tCAACTCTATTTCTCGCGAGTATTTAAAAAATGCa                                                        >  1:3915687/1‑35 (MQ=255)
                                                   tCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGc                                             >  1:989224/1‑46 (MQ=255)
                                                   tCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGc                                             >  1:2713731/1‑46 (MQ=255)
                                                    cAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                    <  1:4105111/54‑1 (MQ=255)
                                                       ctctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcga                                         <  1:450258/46‑1 (MQ=255)
                                                            tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa              <  1:4615251/68‑1 (MQ=255)
                                                            tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa              <  1:1400938/68‑1 (MQ=255)
                                                              tctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGt          >  1:549115/1‑70 (MQ=255)
                                                              tctcGCGAGTATTTAAAAAATGCACAGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGt          <  1:2729000/70‑1 (MQ=255)
                                                               ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg         <  1:583383/70‑1 (MQ=255)
                                                               ctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTg         <  1:2639173/70‑1 (MQ=255)
                                                                   cgAGAATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgtg                    >  1:5367163/1‑55 (MQ=255)
                                                                    gAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgat   >  1:1397413/1‑71 (MQ=255)
                                                                     aGTATTTAAAAAATGCACCGGGGCCAGCACGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgatg  <  1:1097584/71‑1 (MQ=255)
                                                                     |                                                                      
CACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATG  >  AP009048/71106‑71245

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: