Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,218,928 T→C 29.6% I243I (ATT→ATC mlrA → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,218,9280TC29.6% 33.3 / 18.9 27I243I (ATT→ATCmlrADNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (6/13);  new base C (2/6);  total (8/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

AACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCAT  >  W3110S.gb/2218859‑2218989
                                                                     |                                                             
aaCCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCAttt                                                              <  1:1489339/71‑1 (MQ=255)
           aaCAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGt                                                    >  1:1006297/1‑70 (MQ=255)
           aaCAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGt                                                    >  1:467979/1‑70 (MQ=255)
             cagcagCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATCTAATGATTCGTTaa                                                 >  1:1095658/1‑71 (MQ=255)
                    aaCTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTa                                                             <  1:280472/52‑1 (MQ=255)
                      cTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATCTAATGATTCGTTAACAAATgcgc                                        >  1:138629/1‑71 (MQ=255)
                         aCCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATCTAATGATTCGTTAACAAATGCGCtt                                      <  1:1677388/70‑1 (MQ=255)
                          ccAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACa                                               <  1:256522/60‑1 (MQ=255)
                             gTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTg                                 <  1:1060116/71‑1 (MQ=255)
                              ttGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAAt                                                           <  1:317472/44‑1 (MQ=255)
                              ttGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAAt                                                           <  1:608630/44‑1 (MQ=255)
                                ggCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTa                               <  1:1754794/70‑1 (MQ=255)
                                ggCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTa                               >  1:32675/1‑70 (MQ=255)
                                  ctaGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGAtt                                                       <  1:905312/42‑1 (MQ=255)
                                  cAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACaa                            >  1:1520997/1‑71 (MQ=255)
                                   aaGAACAAGGCCATGATATTTTCCCACTCGGCATCTAATGATTCGTTa                                                  <  1:274700/48‑1 (MQ=255)
                                   aaGAACAAGGCCATGATATTTTCCCACTCGGCATCTAATGATTCGTTa                                                  <  1:704313/48‑1 (MQ=255)
                                         aaGGCCATGATATTTTCCCACTCGGCATCTAATGATTCGTTAACaa                                              <  1:1466404/46‑1 (MQ=255)
                                         aaGGCCATGATATTTTCCCACTCGGCATCTAATGATTCGTTAACAAATGCGCTTTACTGTAc                              <  1:1732811/62‑1 (MQ=255)
                                            gCCATGATATTTTCCCACTCGGCATCTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGtt                  <  1:647228/71‑1 (MQ=255)
                                                tGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCtt                                      <  1:1771332/47‑1 (MQ=255)
                                                  atatTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACAt             <  1:35399/70‑1 (MQ=255)
                                                  atatTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATa            <  1:948171/71‑1 (MQ=255)
                                                    atTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACAt             >  1:129296/1‑68 (MQ=255)
                                                     ttttCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCtt                                      <  1:634838/42‑1 (MQ=255)
                                                          ccACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTg     <  1:586485/70‑1 (MQ=255)
                                                            aCTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCAt  <  1:1230152/71‑1 (MQ=255)
                                                             cTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCAt  >  1:1488040/1‑70 (MQ=255)
                                                                     |                                                             
AACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCAT  >  W3110S.gb/2218859‑2218989

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