Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,574,961 T→C 58.4% R212R (CGT→CGC glyS → glycine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,574,9610TC58.4% 24.5 / 45.4 68R212R (CGT→CGCglySglycine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (14/14);  new base C (9/30);  total (24/44)
Fisher's exact test for biased strand distribution p-value = 3.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGCGAGCCGGAATTCACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCGATGCCGAAGAAGC  >  minE/2574900‑2575023
                                                             |                                                              
tGGGCGAGCCGGAATTCACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGaaa                                                          <  1:1767300/68‑1 (MQ=255)
tGGGCGAGCCGGAATTCACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGaaa                                                          <  1:2138781/68‑1 (MQ=255)
tGGGCGAGCCGGAATTCACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGaaa                                                          <  1:2818792/68‑1 (MQ=255)
                cACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGATAGTCATCGCCGATTa                                            >  1:2276609/1‑66 (MQ=255)
                cACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTa                                            >  1:1996328/1‑66 (MQ=255)
                cACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTa                                            >  1:84266/1‑66 (MQ=255)
                cACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTa                                            >  1:2058439/1‑66 (MQ=255)
                cACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTa                                            >  1:2079316/1‑66 (MQ=255)
                cACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTa                                            >  1:323498/1‑66 (MQ=255)
                cACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTa                                            >  1:2122155/1‑66 (MQ=255)
                cACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTa                                            >  1:3085984/1‑66 (MQ=255)
                cACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTa                                            >  1:1478665/1‑66 (MQ=255)
                       gATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAAc                                    >  1:1012317/1‑67 (MQ=255)
                       gATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAAc                                    >  1:543760/1‑67 (MQ=255)
                       gATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAAc                                    >  1:1153500/1‑67 (MQ=255)
                       gATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGt                                  <  1:139361/69‑1 (MQ=255)
                       gATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGt                                  <  1:913388/69‑1 (MQ=255)
                          aatgCCGATCAGTATCCGGAAATTCTGCGTGAGCGGGGGAAAGTCATCGCCGATTACGAAGAACGTAAg                               >  1:1193824/4‑69 (MQ=255)
                          aatgCCGATCAGTATCCGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTaa                                >  1:3070867/4‑68 (MQ=255)
                          aatgCCGATCAGTATCCGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTaa                                >  1:1038953/4‑68 (MQ=255)
                          aatgCCGATCAGTATCCGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAg                               >  1:1570856/4‑69 (MQ=255)
                          aatgCCGATCAGTATCCGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAg                               >  1:2400634/4‑69 (MQ=255)
                          aatgCCGATCAGTATCCGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAg                               >  1:2466092/4‑69 (MQ=255)
                          aatgCCGATCAGTATCCGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAg                               >  1:2593085/4‑69 (MQ=255)
                          aatgCCGATCAGTATCCGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAg                               >  1:2982305/4‑69 (MQ=255)
                          aatgCCGATCAGTATCCGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAg                               >  1:1780939/4‑69 (MQ=255)
                          aatgCCGATCAGTATCCGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAg                               >  1:3393487/4‑69 (MQ=255)
                               cGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCTTCGCCGATTACGAAGAACGTAAGGCGaa                          <  1:397070/69‑1 (MQ=255)
                               cGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGaa                          <  1:3211272/69‑1 (MQ=255)
                               cGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGaa                          <  1:2981195/69‑1 (MQ=255)
                               cGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGaa                          <  1:1537603/69‑1 (MQ=255)
                               cGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGaa                          <  1:2133105/69‑1 (MQ=255)
                               cGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGaa                          <  1:1186369/69‑1 (MQ=255)
                               cGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGaa                          <  1:2522712/69‑1 (MQ=255)
                               cGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGaa                          <  1:2438825/69‑1 (MQ=255)
                                          cGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGcc                <  1:1293934/68‑1 (MQ=255)
                                          cGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGcc                <  1:3058523/68‑1 (MQ=255)
                                          cGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGcc                <  1:3227001/68‑1 (MQ=255)
                                          cGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGcc                <  1:1460484/68‑1 (MQ=255)
                                          cGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGcc                <  1:361088/68‑1 (MQ=255)
                                          cGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGcc                <  1:2156694/68‑1 (MQ=255)
                                          cGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGcc                <  1:11599/68‑1 (MQ=255)
                                          cGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGcc                <  1:1558110/68‑1 (MQ=255)
                                          cGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGcc                <  1:689850/68‑1 (MQ=255)
                                          cGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGcc                <  1:1981128/68‑1 (MQ=255)
                                          cGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGcc                <  1:1866031/68‑1 (MQ=255)
                                          cGGAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGcc                <  1:1825869/68‑1 (MQ=255)
                                           ggAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGcc                <  1:1102920/67‑1 (MQ=255)
                                           ggAAATTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:1992296/68‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:404378/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:1640581/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:447178/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:344989/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:741977/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:1142664/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:3323889/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:326199/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:1367964/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:1402870/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:2383071/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:1464231/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:1600650/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:1937108/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:187158/66‑1 (MQ=255)
                                             aaaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:1733154/66‑1 (MQ=255)
                                              aaTTCTGCGTGAGCGCGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCg               <  1:2181883/65‑1 (MQ=255)
                                                       tGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCGATGCCGaagaag   >  1:2436848/1‑68 (MQ=255)
                                                       tGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCGATGCCGAAGAagc  >  1:9733/1‑69 (MQ=255)
                                                             |                                                              
TGGGCGAGCCGGAATTCACCATCGATAACGCCGATCAGTATCCGGAAATTCTGCGTGAGCGTGGGAAAGTCATCGCCGATTACGAAGAACGTAAGGCGAAGATTAAAGCCGATGCCGAAGAAGC  >  minE/2574900‑2575023

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: