Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,400,733 A→G I70I (ATT→ATC mreB ← cell wall structural complex MreBCD, actin‑like component MreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,400,7330AG92.7% 128.9 / 0.8 41I70I (ATT→ATCmreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base G (22/16);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 9.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCG  >  W3110S.gb/3400665‑3400797
                                                                    |                                                                
tGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGat                                                                <  1:1365950/71‑1 (MQ=255)
tGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGat                                                                <  1:2838337/71‑1 (MQ=255)
tGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGat                                                                <  1:4065547/71‑1 (MQ=255)
tGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGat                                                                <  1:401277/71‑1 (MQ=255)
     gAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTg                                                            <  1:270402/70‑1 (MQ=255)
     gAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTg                                                            <  1:2267851/70‑1 (MQ=255)
     gAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTg                                                            <  1:337770/70‑1 (MQ=255)
     gAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTg                                                            <  1:1683523/70‑1 (MQ=255)
     gAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTg                                                            <  1:1367116/70‑1 (MQ=255)
       gCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTGccc                                                         >  1:3958884/1‑71 (MQ=255)
            ttttCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCg                                                     >  1:1877269/1‑70 (MQ=255)
                     aCGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCa                                            >  1:2193382/1‑70 (MQ=255)
                     aCGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAg                                           >  1:3981581/1‑71 (MQ=255)
                     aCGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAg                                           >  1:250085/1‑71 (MQ=255)
                                   gATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGt                             <  1:3370461/71‑1 (MQ=255)
                                   gATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGt                             <  1:2914335/71‑1 (MQ=255)
                                   gATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGt                             <  1:2130496/71‑1 (MQ=255)
                                   gATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGt                             <  1:3669317/71‑1 (MQ=255)
                                   gATAACGCCGTCTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGt                             <  1:3998081/71‑1 (MQ=255)
                                       aCGCCGTCTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATg                         <  1:2897348/71‑1 (MQ=255)
                                       aCGCCGTCTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATg                         <  1:1858897/71‑1 (MQ=255)
                                          ccGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGAcc                      <  1:96933/71‑1 (MQ=255)
                                          ccGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGAcc                      <  1:333172/71‑1 (MQ=255)
                                          ccGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGAcc                      <  1:1754267/71‑1 (MQ=255)
                                              cTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCtt                                  >  1:3553527/1‑55 (MQ=255)
                                              cTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTAc                   >  1:1891993/1‑70 (MQ=255)
                                              cTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                  >  1:3305929/1‑71 (MQ=255)
                                              cTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                  >  1:83111/1‑71 (MQ=255)
                                              cTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                  >  1:1449405/1‑71 (MQ=255)
                                              cTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                  >  1:1792794/1‑71 (MQ=255)
                                              cTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                  >  1:1793658/1‑71 (MQ=255)
                                              cTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                  >  1:3407205/1‑71 (MQ=255)
                                              cTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                  >  1:3152914/1‑71 (MQ=255)
                                              cTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                  >  1:3374507/1‑71 (MQ=255)
                                              cTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                  >  1:3214096/1‑71 (MQ=255)
                                              cTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                  >  1:3257258/1‑71 (MQ=255)
                                              cTTTCATTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                  >  1:3276787/1‑71 (MQ=255)
                                              cTTTCAGTGGGCGAATGGCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACt                  >  1:526286/1‑71 (MQ=255)
                                                              ggCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTcgcg                              >  1:1057167/1‑43 (MQ=255)
                                                              ggCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCg  >  1:2390394/1‑71 (MQ=255)
                                                              ggCAGCGATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCg  >  1:817598/1‑71 (MQ=255)
                                                                    |                                                                
TGCTGGAGCATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCG  >  W3110S.gb/3400665‑3400797

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