Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 42.9% 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
*W3110S.gb71,1750AG42.9% 19.3 / 62.4 49R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (12/16);  new base G (10/11);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 7.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATG  >  W3110S.gb/71110‑71242
                                                                 |                                                                   
aCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                >  1:1750663/1‑71 (MQ=255)
aCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGAGAGTAtt                                                                >  1:518419/1‑71 (MQ=255)
     gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                           <  1:1787339/71‑1 (MQ=255)
       tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         >  1:1083469/1‑71 (MQ=255)
       tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                         <  1:1537762/71‑1 (MQ=255)
         ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                       >  1:526864/1‑71 (MQ=255)
         ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                       >  1:150117/1‑71 (MQ=255)
         ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                       >  1:1490740/1‑71 (MQ=255)
           tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                       <  1:551496/69‑1 (MQ=255)
           tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                       <  1:1839505/69‑1 (MQ=255)
           tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                       <  1:918553/69‑1 (MQ=255)
             tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAcc                                                   >  1:1612389/1‑71 (MQ=255)
             tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAcc                                                   >  1:510005/1‑71 (MQ=255)
               gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCgg                                                 <  1:1110010/71‑1 (MQ=255)
                   ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                  <  1:143453/50‑1 (MQ=255)
                   ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                  <  1:551546/50‑1 (MQ=255)
                   ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                  <  1:165070/50‑1 (MQ=255)
                   ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                  <  1:543484/50‑1 (MQ=255)
                               aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                 >  1:1542371/1‑71 (MQ=255)
                               aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                 >  1:267294/1‑71 (MQ=255)
                                tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            <  1:947599/43‑1 (MQ=255)
                                tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            <  1:1774848/43‑1 (MQ=255)
                                  ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                              <  1:735374/71‑2 (MQ=255)
                                  ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                              >  1:971182/1‑70 (MQ=255)
                                  ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                         >  1:1222112/1‑44 (MQ=255)
                                  ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                         >  1:1307617/1‑44 (MQ=255)
                                  ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                           <  1:1299293/42‑1 (MQ=255)
                                  ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                           <  1:521714/42‑1 (MQ=255)
                                  ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                           <  1:764148/42‑1 (MQ=255)
                                  ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                           <  1:1319137/42‑1 (MQ=255)
                                  ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                           <  1:926595/42‑1 (MQ=255)
                                    ggTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCg                            >  1:1108031/1‑71 (MQ=255)
                                     gTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCg                            <  1:600954/70‑1 (MQ=255)
                                               tCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGc                                          >  1:1307406/1‑46 (MQ=255)
                                               tCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                 <  1:805096/71‑1 (MQ=255)
                                                cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                 >  1:312466/1‑70 (MQ=255)
                                                cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                 >  1:1279002/1‑70 (MQ=255)
                                                cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgt                  >  1:1092223/1‑69 (MQ=255)
                                                    tctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcga                                  <  1:760220/49‑1 (MQ=255)
                                                    tctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcga                                  <  1:154075/49‑1 (MQ=255)
                                                    tctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgtg                 <  1:1332187/66‑1 (MQ=255)
                                                      tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa           >  1:1710151/1‑70 (MQ=255)
                                                      tATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtg                   >  1:355384/1‑62 (MQ=255)
                                                       aTTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa           >  1:1552660/1‑69 (MQ=255)
                                                        tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa           <  1:1376961/68‑1 (MQ=255)
                                                           ctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTg      >  1:1253016/1‑70 (MQ=255)
                                                              gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg  <  1:953538/71‑1 (MQ=255)
                                                              gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAat   <  1:1489509/70‑1 (MQ=255)
                                                              gcgAGTATTTAAAAAATGCACCGGGGCCAGCCAGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatg  <  1:1096750/71‑1 (MQ=255)
                                                                 |                                                                   
ACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATG  >  W3110S.gb/71110‑71242

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