Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 626,252 T→C 34.6% S55S (AGT→AGC ycaQ → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE626,2520TC34.6% 50.8 / 42.4 75S55S (AGT→AGCycaQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (41/8);  new base C (21/5);  total (62/13)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGATCTCCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCT  >  minE/626195‑626314
                                                         |                                                              
aaCGATCTCCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCata                                                           >  1:1669969/1‑63 (MQ=255)
aaCGATCTCCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCata                                                           >  1:3321397/1‑63 (MQ=255)
aaCGATCTCCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCata                                                           >  1:1110899/1‑63 (MQ=255)
aaCGATCTCCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCata                                                           >  1:1136334/1‑63 (MQ=255)
aaCGATCTCCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCata                                                           >  1:1580546/1‑63 (MQ=255)
aaCGATCTCCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCata                                                           >  1:814954/1‑63 (MQ=255)
aaCGATCTCCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCata                                                           >  1:2045079/1‑63 (MQ=255)
aaCGATCTCCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCata                                                           >  1:2369220/1‑63 (MQ=255)
aaCGATCTCCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCata                                                           >  1:3061850/1‑63 (MQ=255)
      ctcCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCtt                                                 >  1:298074/1‑67 (MQ=255)
      ctcCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCtt                                                 >  1:3319018/1‑67 (MQ=255)
      ctcCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCtt                                                 >  1:3182070/1‑67 (MQ=255)
      ctcCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCtt                                                 >  1:984983/1‑67 (MQ=255)
      ctcCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCtt                                                 >  1:2471243/1‑67 (MQ=255)
      ctcCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCtt                                                 >  1:2915016/1‑67 (MQ=255)
      ctcCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCtt                                                 >  1:2808785/1‑67 (MQ=255)
      ctcCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCtt                                                 >  1:1004560/1‑67 (MQ=255)
      ctcCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCtt                                                 >  1:1624778/1‑67 (MQ=255)
      ctcCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCtt                                                 >  1:436497/1‑67 (MQ=255)
      ctcCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCtt                                                 >  1:1452750/1‑67 (MQ=255)
      ctcCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCtt                                                 >  1:889318/1‑67 (MQ=255)
              tGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAgtcgtc                                       >  1:3257880/1‑69 (MQ=255)
              tGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAgtcgtc                                       >  1:133103/1‑69 (MQ=255)
              tGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAgtcgt                                        >  1:564396/1‑68 (MQ=255)
                tCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCt                                      >  1:3436353/1‑68 (MQ=255)
                tCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCt                                      >  1:966354/1‑68 (MQ=255)
                tCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCt                                      >  1:2862215/1‑68 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAgccg                                         >  1:1685817/1‑55 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTAt                            >  1:1746876/1‑68 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:2738872/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:594451/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:2997019/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:3017598/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:516147/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:484176/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:3362860/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:3420740/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:2237045/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:131070/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:1819781/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:1523200/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:1731285/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:1393987/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:2234190/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:2307125/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:1093865/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           >  1:2577171/1‑69 (MQ=255)
                          aaaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTCTCAGCCGTCTGGGAGATTATc                           >  1:2102096/1‑69 (MQ=255)
                           aaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGTGAGATTATc                           <  1:3004942/68‑1 (MQ=255)
                           aaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           <  1:885545/68‑1 (MQ=255)
                           aaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           <  1:533575/68‑1 (MQ=255)
                           aaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           <  1:2498651/68‑1 (MQ=255)
                           aaTCGATACCATCAATATTGTTGTCCGTAGCCCATATCTGGTGCTTTTCAGCCGTCTGGGAGATTATc                           <  1:5181/68‑1 (MQ=255)
                             tCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTg                        >  1:3360344/1‑69 (MQ=255)
                             tCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTg                        >  1:867053/1‑69 (MQ=255)
                             tCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTg                        >  1:1427952/1‑69 (MQ=255)
                                        aaTATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTg             >  1:349397/1‑69 (MQ=255)
                                        aaTATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTg             >  1:757102/1‑69 (MQ=255)
                                        aaTATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTg             >  1:3272872/1‑69 (MQ=255)
                                        aaTATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTg             >  1:2216587/1‑69 (MQ=255)
                                            ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg         >  1:1306717/1‑69 (MQ=255)
                                            ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg         >  1:1066277/1‑69 (MQ=255)
                                            ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg         >  1:2730976/1‑69 (MQ=255)
                                            ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg         >  1:1653515/1‑69 (MQ=255)
                                            ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg         >  1:3324505/1‑69 (MQ=255)
                                            ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg         >  1:3293259/1‑69 (MQ=255)
                                            ttgttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATg         >  1:2808472/1‑69 (MQ=255)
                                              gttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGa        <  1:1999733/68‑1 (MQ=255)
                                              gttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGa        <  1:2149632/68‑1 (MQ=255)
                                              gttgCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGa        <  1:2574915/68‑1 (MQ=255)
                                                 gCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGtc     <  1:3352669/68‑1 (MQ=255)
                                                 gCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGtc     <  1:3325009/68‑1 (MQ=255)
                                                 gCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGtc     <  1:1357729/68‑1 (MQ=255)
                                                 gCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGtc     <  1:2382793/68‑1 (MQ=255)
                                                    cGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTctct  <  1:1269484/68‑1 (MQ=255)
                                                         |                                                              
AACGATCTCCCGCATGTCCTTGCTGCAAATCGATACCATCAATATTGTTGCCCGTAGTCCATATCTGGTGCTTTTCAGTCGTCTGGGAAATTATCCTGCCCAGTGGCTGGATGAGTCTCT  >  minE/626195‑626314

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: