Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,385,853 (T)8→9 intergenic (+52/+4) yhcN → / ← yhcO DUF1471 domain‑containing stress‑induced protein YhcN/putative barnase inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,385,8451.T88.9% 108.7 / 6.2 36intergenic (+44/+12)yhcN/yhcODUF1471 domain‑containing stress‑induced protein YhcN/putative barnase inhibitor
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base T (20/12);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAA  >  NC_000913/3385725‑3385972
                                                                                                                         |                                                                                                                                  
aaGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCtttttttttgcg                                                                                                                     <  1:530422/133‑5 (MQ=255)
   gcgcAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTa                                                                                                                   <  2:6335/133‑1 (MQ=255)
    cgcAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTaa                                                                                                                  >  2:609428/1‑133 (MQ=255)
         cGGCCTACCAGATTACTGTAGCTCGTAGCGCTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACg                                                                                                             <  1:187128/133‑1 (MQ=255)
                      tACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGc                                                                                                >  2:388349/1‑133 (MQ=255)
                      tACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTCTTTTAGCTAATGGCGAAGATTAAAAccc                                                                                                >  2:299560/1‑131 (MQ=255)
                        cTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCaa                                                                                              >  2:435495/1‑133 (MQ=255)
                              cTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCTCTAATGAGGAAAATCAAAACGAAAATGccc                                                                                        >  1:202663/1‑133 (MQ=255)
                              cTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGccc                                                                                        >  1:124318/1‑133 (MQ=255)
                                  tAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTCAAACGCAAAAGCCCTTcc                                                                                    >  1:212987/1‑133 (MQ=255)
                                  tAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                    >  1:437026/1‑133 (MQ=255)
                                  tAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTcc                                                                                    >  1:613400/1‑133 (MQ=255)
                                      ggTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGct                                                                                <  2:437026/133‑1 (MQ=255)
                                      ggTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGct                                                                                <  2:473961/133‑1 (MQ=255)
                                      ggTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGct                                                                                <  1:351176/133‑1 (MQ=255)
                                        tGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGctct                                                                              <  1:388349/133‑1 (MQ=255)
                                        tGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGctct                                                                              >  1:535725/1‑133 (MQ=255)
                                        tGACACCTAGCACGCGACTGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGctct                                                                              <  1:299560/133‑1 (MQ=255)
                                                ggCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGc                                                                      >  2:558953/1‑133 (MQ=255)
                                                        aCGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCaa                                                              <  1:435495/133‑1 (MQ=255)
                                                                   gTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATAtt                                                   >  1:477657/1‑133 (MQ=255)
                                                                   gTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCACTCTTCCTCTGCCTCATCAAACAGCAATATTa                                                  >  2:427459/1‑133 (MQ=255)
                                                                    tACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTa                                                  >  2:566560/1‑133 (MQ=255)
                                                                              cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGGCTCATCAAACAGCATTATTAATGCGCCaaa                                        >  1:217893/1‑133 (MQ=255)
                                                                              cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAGACAGCATTATTAAGGCGCCaaa                                        >  1:187036/1‑133 (MQ=255)
                                                                              cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCaaa                                        >  2:333660/1‑133 (MQ=255)
                                                                              cccTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCaaa                                        >  2:357663/1‑133 (MQ=255)
                                                                                  cATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGc                                   <  1:107637/133‑1 (MQ=255)
                                                                                     cGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGa                                 >  2:330388/1‑133 (MQ=255)
                                                                                           gTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGTAATATTAGCGCTCCGAAACGGGGACGCGtt                           >  1:192766/1‑133 (MQ=255)
                                                                                              cgacgaTATTTCCCCCGGTTCGGGTGCTCTTTCTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTtc                        <  2:166001/133‑1 (MQ=255)
                                                                                                      ttGCCCCCGGGGCGGTGGCTTTTTTTTTTCGCTAATGACGATCATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAg                <  2:202663/133‑1 (MQ=255)
                                                                                                            ccGGTTCGGGTGCTTATTTTTTGCCCTATTGATGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAAc          <  2:187036/133‑1 (MQ=255)
                                                                                                              ggTTCGGGGGCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGCCAACCCGTTTTCTCTCCCAGATGAAcca        >  2:305787/1‑131 (MQ=255)
                                                                                                                    ggggCTTTTTTTTTGGGCTAATGACGAACATTTAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCaa  <  2:212987/133‑1 (MQ=255)
                                                                                                                    ggggCTTTTTTTTTGCGCTAATGAGGAACATTCAAACGCAAGTGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCaa  <  2:217893/133‑1 (MQ=255)
                                                                                                                    ggggCTTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCaa  <  2:613400/133‑1 (MQ=255)
                                                                                                                         |                                                                                                                                  
AAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCTAATGACGAACATTAAAACGCAAATGCCCTTCCAGCTCTTCCTCTGCCTCATCAAACAGCAATATTAACGCGCCAAAACGGCGACGCGTTTTCTCTCCCAGATGAACAAATTCAA  >  NC_000913/3385725‑3385972

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