Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,332,182 A→G 44.4% A315A (GCT→GCC cpsB ← mannose‑1‑phosphate guanyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,332,1820AG44.4% 11.7 / 57.6 72A315A (GCT→GCCcpsBmannose‑1‑phosphate guanyltransferase
Reads supporting (aligned to +/- strand):  ref base A (26/14);  new base G (12/20);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 3.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGG  >  minE/1332117‑1332246
                                                                 |                                                                
aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCa                                                               <  1:242577/69‑1 (MQ=255)
aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCa                                                               <  1:246823/69‑1 (MQ=255)
aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCa                                                               <  1:3174710/69‑1 (MQ=255)
aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCa                                                               <  1:3238005/69‑1 (MQ=255)
aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCa                                                               <  1:1309103/69‑1 (MQ=255)
aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCa                                                               <  1:1310401/69‑1 (MQ=255)
aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCa                                                               <  1:3295030/69‑1 (MQ=255)
aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCa                                                               <  1:1332024/69‑1 (MQ=255)
aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCa                                                               <  1:3363261/69‑1 (MQ=255)
aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCa                                                               <  1:451879/69‑1 (MQ=255)
aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCa                                                               <  1:1447924/69‑1 (MQ=255)
      tCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGCCGGTGGTGACCAGGCCAGATTCAGCATaca                                                           >  1:2173841/1‑67 (MQ=255)
      tCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATaca                                                           >  1:570908/1‑67 (MQ=255)
      tCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATaca                                                           >  1:1497621/1‑67 (MQ=255)
      tCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATaca                                                           >  1:1396865/1‑67 (MQ=255)
      tCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATaca                                                           >  1:1730724/1‑67 (MQ=255)
      tCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATaca                                                           >  1:677271/1‑67 (MQ=255)
      tCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATaca                                                           >  1:72466/1‑67 (MQ=255)
      tCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATaca                                                           >  1:2812765/1‑67 (MQ=255)
      tCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATaca                                                           >  1:979084/1‑67 (MQ=255)
      tCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATaca                                                           >  1:1056690/1‑67 (MQ=255)
      tCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATa                                                             >  1:3033337/1‑65 (MQ=255)
      tCTTTGGTCTGCACCACTACCAGATCTCTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATaca                                                           >  1:688102/1‑67 (MQ=255)
                              tCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTa                                 <  1:3310594/69‑1 (MQ=255)
                               cTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGAtt                                  >  1:1070326/1‑67 (MQ=255)
                               cTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTa                                 >  1:1153550/1‑68 (MQ=255)
                               cTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTa                                 >  1:2860756/1‑68 (MQ=255)
                               cTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTa                                 >  1:1817981/1‑68 (MQ=255)
                               cTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTa                                 >  1:465881/1‑68 (MQ=255)
                               cTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTa                                 >  1:2501133/1‑68 (MQ=255)
                                                     cAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCa                <  1:2707866/63‑1 (MQ=255)
                                                     cAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCa                <  1:2731328/63‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:3349373/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:3339547/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:1676947/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:1879754/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:1586155/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:1416056/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:3250367/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:1302209/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:1857478/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:887591/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:3216439/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:1874733/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:3164921/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:2994381/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:2354326/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:2029975/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:2556618/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:2568547/46‑1 (MQ=255)
                                                        gCCAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:1895147/46‑1 (MQ=255)
                                                         ccAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATc                              <  1:755074/45‑1 (MQ=255)
                                                         ccAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGccc      >  1:2880747/1‑69 (MQ=255)
                                                         ccAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGccc      >  1:417395/1‑69 (MQ=255)
                                                         ccAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGccc      >  1:145579/1‑69 (MQ=255)
                                                         ccAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGccc      >  1:2017526/1‑69 (MQ=255)
                                                         ccAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGccc      >  1:605054/1‑69 (MQ=255)
                                                          cAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCATGGCAAACGTTGccgt     >  1:2884101/1‑67 (MQ=255)
                                                          cAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCATGGCAAACGTTGccc      >  1:282730/1‑68 (MQ=255)
                                                          cAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCATGGCAAACGTTGccc      >  1:2888048/1‑68 (MQ=255)
                                                          cAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCATGGCAAACGTTGCCCt     >  1:891927/1‑69 (MQ=255)
                                                          cAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCATGGCAAACGTTGCCCt     >  1:894413/1‑69 (MQ=255)
                                                          cAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCATGGCAAACGTTGCCCt     >  1:1245/1‑69 (MQ=255)
                                                          cAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCATGGCAAACGTTGCCCt     >  1:2256002/1‑69 (MQ=255)
                                                          cAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCATGGCAAACGTTGCCCt     >  1:1317233/1‑69 (MQ=255)
                                                          cAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCATGGCAAACGTTGCCCt     >  1:590602/1‑69 (MQ=255)
                                                          cAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCATGGCAAACGTTGCCCt     >  1:372471/1‑69 (MQ=255)
                                                          cAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCATGGCAAACGTTGCCCt     >  1:3360490/1‑69 (MQ=255)
                                                          cAGATTCGGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCATGGCAAACGTTGCCCt     >  1:2127714/1‑69 (MQ=255)
                                                              ttCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGtt          >  1:2638752/1‑60 (MQ=255)
                                                              ttCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTcgg  >  1:322672/1‑68 (MQ=255)
                                                              ttCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGAGGCAAACGTTGCCCTcgg  >  1:2584829/1‑68 (MQ=255)
                                                                 |                                                                
ACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGG  >  minE/1332117‑1332246

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