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
                                                                                                               |                                                                                                                                  
cACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                                        <  1:1442729/139‑1 (MQ=255)
cACCGCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAg                                                                                                        <  1:1494755/139‑1 (MQ=255)
    gCTCTTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAt                                                                                                    <  2:227246/139‑1 (MQ=255)
     ctctTATCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtc                                                                                                   >  1:1221509/1‑139 (MQ=255)
          aTCGCGAAAAACGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTa                                                                                              >  2:787622/1‑139 (MQ=255)
                aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatcat                                                                                                                      >  2:563199/1‑109 (MQ=255)
                aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatcat                                                                                                                      <  1:563199/109‑1 (MQ=255)
                aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCcgcg                                                                                                            <  2:983368/119‑1 (MQ=255)
                aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCcgcg                                                                                                            >  1:983368/1‑119 (MQ=255)
                aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtct                                                                                                  >  2:830428/1‑129 (MQ=255)
                aaaaaCGTAACGGTCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCAtct                                                                                                  <  1:830428/129‑1 (MQ=255)
                             tCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                      >  1:1059929/1‑112 (MQ=255)
                             tCCTCTCCTTTGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGcc                                                                                                      <  2:1059929/112‑1 (MQ=255)
                                      ttGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATc                                                                  <  1:787622/139‑1 (MQ=255)
                                      ttGTGAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATc                                                                  <  1:1513993/139‑1 (MQ=255)
                                       tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATcc                                                                 >  2:777064/1‑139 (MQ=255)
                                       tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATcc                                                                 >  1:774185/1‑139 (MQ=255)
                                       tgtgAGAAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATcc                                                                 >  2:1313122/1‑139 (MQ=255)
                                           agaAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGtt                                                             >  2:381449/1‑139 (MQ=255)
                                            gaAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCgg                                                                                                                   <  1:844264/84‑1 (MQ=255)
                                            gaAATCAAATTACCGCGTAGTGAAACCGAAGCGCCATCATGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCgg                                                                                                                   >  2:844264/1‑84 (MQ=255)
                                              aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCAt                                                                   <  1:457611/130‑1 (MQ=255)
                                              aaTCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCAt                                                                   >  2:457611/1‑130 (MQ=255)
                                                tCAAATTACCGCGTAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTaaa                                                        <  2:774185/139‑1 (MQ=255)
                                                tCAAATTACCGCGTAATGAAACCGTAACGCCTTCGGGCATGGTTTTTGCGAAATCAGCGGGC‑CTTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTaaa                                                        <  1:620800/139‑1 (MQ=255)
                                                             tAGTGAAACCGAAGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCAt                                           <  1:381449/139‑1 (MQ=255)
                                                                        aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatcat                                                                                                                      >  1:945982/1‑53 (MQ=255)
                                                                        aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGcatcat                                                                                                                      <  2:945982/53‑1 (MQ=255)
                                                                        aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTg                                                                            <  2:1340275/95‑1 (MQ=255)
                                                                        aaGCGCCATCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTg                                                                            >  1:1340275/1‑95 (MQ=255)
                                                                        aaGCGCCATCGTGCATGGATTTTGCGAAATCAACGGTCA‑TTTTTTTGCAACATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCa                                >  1:1450431/1‑139 (MQ=255)
                                                                          gcgcCATCGTGCATGGTTTTTGCG‑‑ATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAg                            >  2:1441909/1‑139 (MQ=255)
                                                                                tCGTGCATGGTTTTTGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCAt                        >  1:941999/1‑139 (MQ=255)
                                                                                            tttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGa                                                                                                           <  2:476695/44‑1 (MQ=38)
                                                                                            tttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGa                                                                                                           >  1:476695/1‑44 (MQ=38)
                                                                                            tttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaaa              <  1:137317/137‑1 (MQ=255)
                                                                                            tttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGaaaa              >  2:137317/1‑137 (MQ=255)
                                                                                            tttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACt            <  1:696195/139‑1 (MQ=255)
                                                                                             ttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAggcgg                                                                                 <  2:524811/69‑1 (MQ=255)
                                                                                             ttGCGAAATCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAggcgg                                                                                 >  1:524811/1‑69 (MQ=255)
                                                                                                  aaaTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAggc                                                                                   >  1:196476/1‑62 (MQ=255)
                                                                                                  aaaTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAggc                                                                                   <  2:196476/62‑1 (MQ=255)
                                                                                                  aaaTCAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTaa      >  2:874974/1‑139 (MQ=255)
                                                                                                      cAACGGTCA‑TTTTTTTGCATCATCGGTTCCGCGATAGCCATCTTCTATTGCATGAGGCGGCGGTGGCGCTGCATCCTGTTTTAAACCGCCCTGGTCATCTGCCAACGCATAAGGCATGACAAGAAAACTTGCTAATACa  <  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: ‑