Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 100,925 C→G 55.8% A305A (GCG→GCC hofB ← conserved hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE100,9250CG55.8% 11.8 / 55.7 52A305A (GCG→GCChofBconserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base C (7/16);  new base G (14/15);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGATCTCTGCTGTTTCGCCATCGCGGATCTCTCCGATCATGATGACGTCAGGAT  >  minE/100858‑100989
                                                                   |                                                                
aaCGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGcgccg                                                                 >  1:3198966/1‑69 (MQ=255)
aaCGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGcgccg                                                                 >  1:2568892/1‑69 (MQ=255)
aaCGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGcgccg                                                                 >  1:331217/1‑69 (MQ=255)
aaCGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGcgccg                                                                 >  1:1433344/1‑69 (MQ=255)
aaCGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGcgccg                                                                 >  1:345757/1‑69 (MQ=255)
aaCGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGcgccg                                                                 >  1:506944/1‑69 (MQ=255)
aaCGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGcgcc                                                                  >  1:35778/1‑68 (MQ=255)
    tACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCttt                                                             <  1:3397730/69‑1 (MQ=255)
    tACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCttt                                                             <  1:2454764/69‑1 (MQ=255)
    tACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCttt                                                             <  1:517216/69‑1 (MQ=255)
    tACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCttt                                                             <  1:2282538/69‑1 (MQ=255)
    tACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCttt                                                             <  1:2767818/69‑1 (MQ=255)
    tACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCttt                                                             <  1:3050163/69‑1 (MQ=255)
    tACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCttt                                                             <  1:2926524/69‑1 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                                             >  1:1710850/1‑63 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                                             >  1:963472/1‑63 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                                             >  1:2249621/1‑63 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                                             >  1:2160890/1‑63 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                                             >  1:286603/1‑63 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                                             >  1:2943886/1‑63 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                                             >  1:3032615/1‑63 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                                             >  1:1188463/1‑63 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                                             >  1:1165873/1‑63 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                                             >  1:3363881/1‑63 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                                             >  1:14690/1‑63 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                                             >  1:3372810/1‑63 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                                             >  1:486503/1‑63 (MQ=255)
          cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCtt                                                              >  1:2429996/1‑62 (MQ=255)
            tttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACGCCAGGTGACCAGTTTGCGCCGCTTTAATAGCGa                                                     >  1:2747287/1‑69 (MQ=255)
            tttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACTAGTTTGCGCCGCTTTAATAGCGa                                                     >  1:2560765/1‑69 (MQ=255)
            tttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGa                                                     >  1:1487417/1‑69 (MQ=255)
            tttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGa                                                     >  1:1317629/1‑69 (MQ=255)
            tttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGa                                                     >  1:2724292/1‑69 (MQ=255)
            tttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGa                                                     >  1:1177518/1‑69 (MQ=255)
            tttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGa                                                     >  1:1179684/1‑69 (MQ=255)
             ttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGAt                                                    <  1:2943869/69‑1 (MQ=255)
             ttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGAt                                                    <  1:1744569/69‑1 (MQ=255)
             ttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGAt                                                    <  1:301933/69‑1 (MQ=255)
             ttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGAt                                                    <  1:397864/69‑1 (MQ=255)
             ttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGAt                                                    <  1:1878140/69‑1 (MQ=255)
             ttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGAt                                                    <  1:261515/69‑1 (MQ=255)
                gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                                  <  1:1641945/68‑3 (MQ=255)
                gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                                  <  1:1194674/68‑3 (MQ=255)
                gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                                  <  1:599672/68‑3 (MQ=255)
                gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                                  <  1:1262631/68‑3 (MQ=255)
                gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                                  <  1:1306135/68‑3 (MQ=255)
                gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                                  <  1:883034/68‑3 (MQ=255)
                gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                                  <  1:1158697/68‑3 (MQ=255)
                gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                                  <  1:1334957/68‑3 (MQ=255)
                gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                                  <  1:15095/68‑3 (MQ=255)
                gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                                  <  1:164918/68‑3 (MQ=255)
                gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                                  <  1:1805812/68‑3 (MQ=255)
                gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                                  <  1:2190730/68‑3 (MQ=255)
                 cAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCTCGGCTTTAATGGCAAttt                                                  <  1:626443/67‑3 (MQ=255)
                 cAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                                  <  1:2987855/67‑3 (MQ=255)
                                   aGGGTAGACAACACCAGGTGACCGGTTTGCGCGTCTTTAATGGCAAttt                                                  <  1:75487/49‑3 (MQ=21)
                                             acacCAGGTGACCAGTTTGCGCCGCTTTAATAGCGATCTCTGCTGTTTCGCCATCGCGGATCTCTCCg                     <  1:2491074/68‑1 (MQ=255)
                                                 cAGGTGACCAGTTTGCGCCGCTTTAATAGCGATCTCTGCTGTTTCGCCATCGCGGATCTCTCCGATCat                <  1:2011048/69‑1 (MQ=255)
                                                 cAGGTGACCAGTTTGCGCCGCTTTAATAGCGATCTCTGCTGTTTCGCCATCGCGGATCTCTCCGATCat                <  1:2240691/69‑1 (MQ=255)
                                                                cgccgcTTTAATAGCGATCTCTGCTGTTTCGCCATCGCGGATCTCTCCGATCATGATGACGTCAGGAt  <  1:1123555/68‑1 (MQ=255)
                                                                cgccgcTTTAATAGCGATCTCTGCTGTTTCGCCATCGCGGATCTCTCCGATCATGATGACGTCAGGAt  <  1:2575510/68‑1 (MQ=255)
                                                                   |                                                                
AACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGATCTCTGCTGTTTCGCCATCGCGGATCTCTCCGATCATGATGACGTCAGGAT  >  minE/100858‑100989

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: