Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,158,370 G→A V241I (GTT→ATT)  yqhG → DUF3828 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,158,3700GA97.8% 141.6 / ‑4.0 46V241I (GTT→ATT) yqhGDUF3828 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (25/20);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

TTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAA  >  NC_000913/3158248‑3158496
                                                                                                                          |                                                                                                                              
ttAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTc                                                                                                                    <  1:900277/135‑1 (MQ=255)
ttAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTc                                                                                                                    >  2:900277/1‑135 (MQ=255)
          cTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTc                                                                                                                    >  1:1108571/1‑125 (MQ=255)
          cTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGaaa                                                                                                      >  1:321966/1‑139 (MQ=255)
          cTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGAAGAAATTc                                                                                                                    <  2:1108571/125‑1 (MQ=255)
                   caccTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGAt                                                                                             <  2:321966/139‑1 (MQ=255)
                     ccTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTc                                                                                                               <  2:234207/119‑1 (MQ=255)
                     ccTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTc                                                                                                               >  1:234207/1‑119 (MQ=255)
                     ccTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTc                                                                                           >  2:294358/1‑139 (MQ=255)
                     ccTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTc                                                                                           >  2:217265/1‑139 (MQ=255)
                       ttCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTccg                                                                                         >  2:676536/1‑139 (MQ=255)
                             aCGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCa                                                                                   <  1:217265/139‑1 (MQ=255)
                                gCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGa                                                                                                                  <  2:1669862/105‑1 (MQ=255)
                                gCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGa                                                                                                                  >  1:1669862/1‑105 (MQ=255)
                                          ccGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTccgc                                                                                        >  2:1381774/1‑121 (MQ=255)
                                          ccGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTccgc                                                                                        <  1:1381774/121‑1 (MQ=255)
                                                  gATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGa                                                              >  1:844044/1‑139 (MQ=255)
                                                       gCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAAt                                                         <  1:294358/139‑1 (MQ=255)
                                                             aGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGcaac                                                   >  2:857301/1‑139 (MQ=255)
                                                               aaCGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGaa                                                                                                       >  2:1421087/1‑85 (MQ=255)
                                                               aaCGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGaa                                                                                                       <  1:1421087/85‑1 (MQ=255)
                                                                  ggCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGa                                              <  1:857301/139‑1 (MQ=255)
                                                                       tttatttaCGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCt                                                                >  1:1445153/1‑116 (MQ=255)
                                                                       tttatttaCGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCt                                                                <  2:1445153/116‑1 (MQ=255)
                                                                                  cAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGttt                              >  1:32038/1‑139 (MQ=255)
                                                                                        gATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAg                        <  1:676536/139‑1 (MQ=255)
                                                                                         aTCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAg                                                                                                                 >  1:186222/1‑49 (MQ=255)
                                                                                         aTCTGTTTTGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAg                                                                                                                 <  2:186222/49‑1 (MQ=255)
                                                                                                ttGGCGAAGAGTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGa                >  2:270722/1‑139 (MQ=255)
                                                                                                        gagTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAGCAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGt        >  1:107752/1‑139 (MQ=255)
                                                                                                        gagTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGAt                                                                                             >  2:968768/1‑54 (MQ=255)
                                                                                                        gagTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGAt                                                                                             <  1:968768/54‑1 (MQ=255)
                                                                                                        gagTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTccgcc                                                                                       >  2:245215/1‑60 (MQ=255)
                                                                                                        gagTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTccgcc                                                                                       <  1:245215/60‑1 (MQ=255)
                                                                                                        gagTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCAt                                                                                  >  2:200657/1‑65 (MQ=255)
                                                                                                        gagTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCAt                                                                                  <  1:200657/65‑1 (MQ=255)
                                                                                                        gagTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGca                                                     <  1:936157/94‑1 (MQ=255)
                                                                                                        gagTCGCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGca                                                     >  2:936157/1‑94 (MQ=255)
                                                                                                             gCTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGaa   >  1:1399781/1‑139 (MQ=255)
                                                                                                              cTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAgg             <  2:141785/128‑1 (MQ=255)
                                                                                                              cTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAgg             >  1:141785/1‑128 (MQ=255)
                                                                                                              cTATTTTCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGaaa  >  2:1435156/1‑139 (MQ=255)
                                                                                                                   ttCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGa                                              >  1:370093/1‑90 (MQ=255)
                                                                                                                   ttCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGa                                              <  2:370093/90‑1 (MQ=255)
                                                                                                                    tCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGa                                     <  1:704522/98‑1 (MQ=255)
                                                                                                                    tCACGCATTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGa                                     >  2:704522/1‑98 (MQ=255)
                                                                                                                          |                                                                                                                              
TTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAA  >  NC_000913/3158248‑3158496

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: