Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,624,357 (T)8→7 coding (141/393 nt) ydeI ← hydrogen peroxide resistance OB fold protein, putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,624,3500T.97.7% 172.6 / ‑2.9 44coding (148/393 nt)ydeIhydrogen peroxide resistance OB fold protein, putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (22/21);  minor base C (0/1);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTAATACA  >  NC_000913/1624239‑1624480
                                                                                                               |                                                                                                                                  
cACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAg                                                                                                        <  1:1442729/139‑1 (MQ=255)
cACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAg                                                                                                        <  1:1494755/139‑1 (MQ=255)
    gCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCAt                                                                                                    <  2:227246/139‑1 (MQ=255)
     ctctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                                   >  1:1221509/1‑139 (MQ=255)
          aTCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTa                                                                                              >  2:787622/1‑139 (MQ=255)
                aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGcatcat                                                                                                                      >  2:563199/1‑109 (MQ=255)
                aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGcatcat                                                                                                                      <  1:563199/109‑1 (MQ=255)
                aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCcgcg                                                                                                            <  2:983368/119‑1 (MQ=255)
                aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCcgcg                                                                                                            >  1:983368/1‑119 (MQ=255)
                aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCAtct                                                                                                  >  2:830428/1‑129 (MQ=255)
                aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCAtct                                                                                                  <  1:830428/129‑1 (MQ=255)
                             tCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                      >  1:1059929/1‑112 (MQ=255)
                             tCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                      <  2:1059929/112‑1 (MQ=255)
                                      ttGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATc                                                                  <  1:787622/139‑1 (MQ=255)
                                      ttGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATc                                                                  <  1:1513993/139‑1 (MQ=255)
                                       tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATcc                                                                 >  2:777064/1‑139 (MQ=255)
                                       tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATcc                                                                 >  1:774185/1‑139 (MQ=255)
                                       tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATcc                                                                 >  2:1313122/1‑139 (MQ=255)
                                           agaAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGtt                                                             >  2:381449/1‑139 (MQ=255)
                                            gaAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCgg                                                                                                                   <  1:844264/84‑1 (MQ=255)
                                            gaAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCATGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCgg                                                                                                                   >  2:844264/1‑84 (MQ=255)
                                              aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCAt                                                                   <  1:457611/130‑1 (MQ=255)
                                              aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCAt                                                                   >  2:457611/1‑130 (MQ=255)
                                                tCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTaaa                                                        <  2:774185/139‑1 (MQ=255)
                                                tCAAATTACCGCGTAATGAAACCGTAACGCCTTCGGGCATGGTTTTTGCGAAATCAGCGGGCCTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTaaa                                                        <  1:620800/139‑1 (MQ=255)
                                                             tAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCAt                                           <  1:381449/139‑1 (MQ=255)
                                                                        aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGcatcat                                                                                                                      >  1:945982/1‑53 (MQ=255)
                                                                        aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGcatcat                                                                                                                      <  2:945982/53‑1 (MQ=255)
                                                                        aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTg                                                                            <  2:1340275/95‑1 (MQ=255)
                                                                        aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTg                                                                            >  1:1340275/1‑95 (MQ=255)
                                                                        aaGCGCCATCGTGCATGGATTTTGCGAAATCAACGGTCATTTTTTTGCAACATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa                                >  1:1450431/1‑139 (MQ=255)
                                                                          gcgcCATCGTGCATGGTTTTTGCG‑‑ATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAg                            >  2:1441909/1‑139 (MQ=255)
                                                                                tCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAt                        >  1:941999/1‑139 (MQ=255)
                                                                                            tttGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGa                                                                                                           <  2:476695/44‑1 (MQ=38)
                                                                                            tttGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGa                                                                                                           >  1:476695/1‑44 (MQ=38)
                                                                                            tttGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaaa              <  1:137317/137‑1 (MQ=255)
                                                                                            tttGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaaa              >  2:137317/1‑137 (MQ=255)
                                                                                            tttGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACt            <  1:696195/139‑1 (MQ=255)
                                                                                             ttGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAggcgg                                                                                 <  2:524811/69‑1 (MQ=255)
                                                                                             ttGCGAAATCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAggcgg                                                                                 >  1:524811/1‑69 (MQ=255)
                                                                                                  aaaTCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAggc                                                                                   >  1:196476/1‑62 (MQ=255)
                                                                                                  aaaTCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAggc                                                                                   <  2:196476/62‑1 (MQ=255)
                                                                                                  aaaTCAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTaa      >  2:874974/1‑139 (MQ=255)
                                                                                                      cAACGGTCATTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTAATACa  <  1:874974/139‑1 (MQ=255)
                                                                                                               |                                                                                                                                  
CACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCATTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTAATACA  >  NC_000913/1624239‑1624480

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