Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 100,916 A→G 44.6% T308T (ACT→ACC 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,9160AG44.6% 9.6 / 51.5 65T308T (ACT→ACChofBconserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base A (20/16);  new base G (14/15);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGATCTCTGCTGTTTCGCCATCGCGGATCTCTCCGATCAT  >  minE/100851‑100975
                                                                 |                                                           
tGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGtt                                                          >  1:2421003/1‑69 (MQ=255)
tGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGtt                                                          >  1:3165026/1‑69 (MQ=255)
tGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGtt                                                          >  1:3313313/1‑69 (MQ=255)
tGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGt                                                           >  1:2274809/1‑68 (MQ=255)
tGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTATAGGGTAGACAACACCAGGTGACCAGtt                                                          >  1:1233744/1‑69 (MQ=255)
tGCTGTAAACATACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGtt                                                          >  1:1519222/1‑69 (MQ=255)
       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:2282538/69‑1 (MQ=255)
           tACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCttt                                               <  1:517216/69‑1 (MQ=255)
           tACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCttt                                               <  1:3050163/69‑1 (MQ=255)
           tACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCttt                                               <  1:2454764/69‑1 (MQ=255)
           tACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCttt                                               <  1:2926524/69‑1 (MQ=255)
           tACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCttt                                               <  1:2767818/69‑1 (MQ=255)
                 cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                               >  1:2943886/1‑63 (MQ=255)
                 cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                               >  1:963472/1‑63 (MQ=255)
                 cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                               >  1:286603/1‑63 (MQ=255)
                 cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                               >  1:3363881/1‑63 (MQ=255)
                 cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                               >  1:3032615/1‑63 (MQ=255)
                 cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                               >  1:486503/1‑63 (MQ=255)
                 cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                               >  1:1165873/1‑63 (MQ=255)
                 cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                               >  1:3372810/1‑63 (MQ=255)
                 cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                               >  1:14690/1‑63 (MQ=255)
                 cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                               >  1:2160890/1‑63 (MQ=255)
                 cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                               >  1:2249621/1‑63 (MQ=255)
                 cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                               >  1:1710850/1‑63 (MQ=255)
                 cGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCttt                                               >  1:1188463/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:1179684/1‑69 (MQ=255)
                   tttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGa                                       >  1:1317629/1‑69 (MQ=255)
                   tttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGa                                       >  1:1487417/1‑69 (MQ=255)
                   tttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGa                                       >  1:1177518/1‑69 (MQ=255)
                   tttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGa                                       >  1:2724292/1‑69 (MQ=255)
                    ttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGAt                                      <  1:1878140/69‑1 (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:261515/69‑1 (MQ=255)
                    ttCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGAt                                      <  1:397864/69‑1 (MQ=255)
                       gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                    <  1:164918/68‑3 (MQ=255)
                       gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                    <  1:1194674/68‑3 (MQ=255)
                       gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                    <  1:599672/68‑3 (MQ=255)
                       gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                    <  1:883034/68‑3 (MQ=255)
                       gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                    <  1:1262631/68‑3 (MQ=255)
                       gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                    <  1:1158697/68‑3 (MQ=255)
                       gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                    <  1:1306135/68‑3 (MQ=255)
                       gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                    <  1:1334957/68‑3 (MQ=255)
                       gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                    <  1:15095/68‑3 (MQ=255)
                       gCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCGGTTTGCGCGGCTTTAATGGCAAttt                                    <  1:1641945/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)
                                                                 |                                                           
TGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTGCGCCGCTTTAATAGCGATCTCTGCTGTTTCGCCATCGCGGATCTCTCCGATCAT  >  minE/100851‑100975

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: