Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,195,442 T→C 40.0% L570L (CTA→CTG yoaA ← conserved hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,195,4420TC40.0% 27.8 / 47.4 70L570L (CTA→CTGyoaAconserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base T (23/19);  new base C (13/15);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

ACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAG  >  minE/1195379‑1195507
                                                               |                                                                 
aCGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTg                                                              <  1:1872893/69‑1 (MQ=255)
     cgtcgGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACtt                                                         <  1:1028066/69‑1 (MQ=255)
     cgtcgGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACtt                                                         <  1:2419128/69‑1 (MQ=255)
     cgtcgGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACtt                                                         <  1:1079845/69‑1 (MQ=255)
               gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCa                                                       >  1:1156197/1‑61 (MQ=255)
               gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCCTCGCACgg                                               >  1:2063533/1‑69 (MQ=255)
               gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACgg                                               >  1:971574/1‑69 (MQ=255)
               gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACgg                                               >  1:167126/1‑69 (MQ=255)
               gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACgg                                               >  1:2330523/1‑69 (MQ=255)
               gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACgg                                               >  1:2945335/1‑69 (MQ=255)
               gcgAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACgg                                               >  1:3149768/1‑69 (MQ=255)
                 gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:2927747/68‑2 (MQ=255)
                 gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:2795261/68‑2 (MQ=255)
                 gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:273599/68‑2 (MQ=255)
                 gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:2719335/68‑2 (MQ=255)
                 gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:2511617/68‑2 (MQ=255)
                 gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:1081096/68‑2 (MQ=255)
                 gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:2063211/68‑2 (MQ=255)
                 gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:1954520/68‑2 (MQ=255)
                 gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:3177260/68‑2 (MQ=255)
                 gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:2442992/68‑2 (MQ=255)
                 gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:1279117/68‑2 (MQ=255)
                 gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:1521142/68‑2 (MQ=255)
                 gAATCAATCGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:1559398/68‑2 (MQ=255)
                  aatcaatcGCCCTACCCCCTGCTTGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:1550896/67‑2 (MQ=255)
                                       cttGAGAGTAATGACCGCATCCGGCAGTTGTACTTCATCGAACgga                                              <  1:2088320/45‑2 (MQ=25)
                                         tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                      >  1:3281235/1‑68 (MQ=255)
                                         tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                      >  1:2753097/1‑68 (MQ=255)
                                         tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                      >  1:1388854/1‑68 (MQ=255)
                                         tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                      >  1:1309215/1‑68 (MQ=255)
                                         tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                      >  1:2606907/1‑68 (MQ=255)
                                         tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                      >  1:2279381/1‑68 (MQ=255)
                                         tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                      >  1:2046670/1‑68 (MQ=255)
                                         tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                      >  1:1572431/1‑68 (MQ=255)
                                         tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                      >  1:3343862/1‑68 (MQ=255)
                                         tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                      >  1:369364/1‑68 (MQ=255)
                                         tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                      >  1:792803/1‑68 (MQ=255)
                                         tGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCt                      >  1:1109759/1‑68 (MQ=255)
                                           agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGTTCGCCACCGCGCAAACGACAATCTTc                    <  1:1027920/68‑1 (MQ=255)
                                           agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                    <  1:1737149/68‑1 (MQ=255)
                                           agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                    <  1:3299977/68‑1 (MQ=255)
                                           agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                    <  1:1205203/68‑1 (MQ=255)
                                           agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                    <  1:3168337/68‑1 (MQ=255)
                                           agagTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                    <  1:1497078/68‑1 (MQ=255)
                                                   gACGGCATCTGGTAGTTGCACTTCATCTAACGGGTCGCCACCGCGCAAACGACAATCTTCCATgcgcg            <  1:1240699/68‑1 (MQ=255)
                                                   gACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATgcgcg            <  1:592706/68‑1 (MQ=255)
                                                   gACGGCATCTGGTAGTTGCACTTCATCGAACGGGGCGCCACCGCGCAAACGACAATCTTc                    <  1:1649760/60‑1 (MQ=255)
                                                       gCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCtttt       <  1:2651409/69‑1 (MQ=255)
                                                       gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTcc                   >  1:2750770/1‑57 (MQ=255)
                                                       gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt       >  1:1385147/1‑69 (MQ=255)
                                                       gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt       >  1:1766677/1‑69 (MQ=255)
                                                       gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt       >  1:765249/1‑69 (MQ=255)
                                                       gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt       >  1:646772/1‑69 (MQ=255)
                                                       gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt       >  1:420987/1‑69 (MQ=255)
                                                       gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt       >  1:371337/1‑69 (MQ=255)
                                                       gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt       >  1:369199/1‑69 (MQ=255)
                                                       gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt       >  1:3382154/1‑69 (MQ=255)
                                                       gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt       >  1:3059343/1‑69 (MQ=255)
                                                       gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt       >  1:2619957/1‑69 (MQ=255)
                                                       gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt       >  1:2940174/1‑69 (MQ=255)
                                                       gCATCCGGCAGTTGTACTTCATCGAACGGATCACCACCGCGCAAACGGCAATCTTCCATGCGTGCtttt       >  1:1372538/1‑69 (MQ=255)
                                                          tCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAAc    <  1:3023428/69‑1 (MQ=255)
                                                            tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACaa                         >  1:297784/1‑46 (MQ=255)
                                                            tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                    <  1:221187/51‑1 (MQ=255)
                                                            tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                    <  1:2266175/51‑1 (MQ=255)
                                                            tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                    <  1:3388037/51‑1 (MQ=255)
                                                            tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTc                    <  1:2777270/51‑1 (MQ=255)
                                                            tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg  >  1:1580254/1‑69 (MQ=255)
                                                            tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg  >  1:150283/1‑69 (MQ=255)
                                                            tGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAg  >  1:1848657/1‑69 (MQ=255)
                                                               |                                                                 
ACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAG  >  minE/1195379‑1195507

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: