Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,700,461 (T)5→4 intergenic (‑264/‑59) shoB ← / → ohsC toxic peptide ShoB/small regulatory RNA OhsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,700,4570T.97.5% 157.2 / ‑1.5 40intergenic (‑260/‑63)shoB/ohsCtoxic peptide ShoB/small regulatory RNA OhsC
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (12/27);  total (13/27)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

AATGGGGGAATTCTGTGGCATGGTGTAAAGCACAGCAAAATCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGTTCCTTACCAGTTC  >  NC_000913/2700326‑2700582
                                                                                                                                   |                                                                                                                             
aaTGGGGGAATTCTGTGGCATGGTGTAAAGCACAGCAAAATCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGtt                                                                                                                              >  2:35118/1‑133 (MQ=255)
  tGGGGGAATTCTGTGGCATGGTGTAAAGCACAGCAAAATCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGtttt                                                                                                                            <  2:13034/133‑1 (MQ=255)
      ggAATTCTGTGGCATGGTGTAAAGCACAGCAAAATCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATTTTGCTTTCTTCGTAAAGGCTCAAGttttatgc                                                                                                                        >  2:173719/1‑129 (MQ=255)
          ttCTGTGGCATGGTGTAAAGCACAGCAAAATCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATAAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTaaa                                                                                                                   >  1:513913/1‑133 (MQ=255)
                   aTGGTGTAAAGCACAGCAAAATCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAg                                                                                                          <  2:375030/133‑1 (MQ=255)
                       tgtAAAGCACAGCAAAATCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGGTCGAGTTGTATGCTGAAGATTCCAAGTTGc                                                                                                      >  1:159446/1‑133 (MQ=255)
                          aaaGCACATCAAAATGTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTg                                                                                                   <  1:220881/133‑1 (MQ=255)
                             gCACAGCAAAATCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTaa                                                                                                <  2:291558/133‑1 (MQ=255)
                               acaGCAAACTCTTCAATAACGACGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTaaaa                                                                                              <  1:173974/133‑1 (MQ=255)
                               acaGCAAAATTTTCATTAACGAAGCCAAATTTAATGTACTTTCGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTaaaa                                                                                              <  2:172456/133‑1 (MQ=255)
                               acaGCAAAATCTTCATTAACGAAGCCAAATTTCATGTACTTCCGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTaaaa                                                                                              <  2:159446/133‑1 (MQ=255)
                               acaGCAAAATCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTaaaa                                                                                              <  2:416461/133‑1 (MQ=255)
                               acaGCAAAATCTTCAATAACGAAACCAAATTTTATGTACTTTCGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTaaaa                                                                                              <  2:253555/133‑1 (MQ=255)
                                 aGCAAAATCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGTCTCGAGTTTTATGCTACAGATTGCACGTTGGTTGTAAAAGa                                                                                            >  2:186598/1‑133 (MQ=255)
                                 aGCAAAATCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTTTGCTAAAGATTGCATGTTGCTCGTAAAAGa                                                                                            >  1:180021/1‑133 (MQ=255)
                                 aGCAAAATCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGa                                                                                            <  2:355595/133‑1 (MQ=255)
                                  gCAAAATCTTAAATATCGAAGCCACTTTTAATTTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGAt                                                                                           <  2:208845/133‑1 (MQ=255)
                                    aaaaTCTTGAATAACAAACCCATTTTTAATCTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATaa                                                                                         <  1:249571/133‑1 (MQ=255)
                                    aaaaTCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTTGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATaa                                                                                         <  1:361286/133‑1 (MQ=255)
                                    aaaaTCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATaa                                                                                         <  1:35118/133‑1 (MQ=255)
                                      aaTCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGt                                                                                       >  1:512103/1‑133 (MQ=255)
                                         cTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTaca                                                                                    <  2:512103/133‑1 (MQ=255)
                                          ttCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTacac                                                                                   <  1:439219/133‑1 (MQ=255)
                                           tCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACt                                                                                  <  1:487733/133‑1 (MQ=255)
                                                    aaGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCCTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAACAGACAACTACACTGATCCATaa                                                                         >  1:239222/1‑133 (MQ=255)
                                                          aTTTTAATGTACTTATGACTTGACGGTCACTTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCtg                                                                   <  1:171061/133‑1 (MQ=255)
                                                                aTGTACTTACTAATTAGCGGTCACTTTATTCTTTCCATCAACATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGttgttg                                                             <  2:154235/133‑1 (MQ=255)
                                                                aTGTACTTACGAATTTGCGGTCACTTTAGTCTCTCCATCAAGATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGttgttg                                                             <  2:173212/133‑1 (MQ=255)
                                                                aTGTACTTACGAATTGGCGGTCTCGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGttgttg                                                             <  1:173719/133‑1 (MQ=255)
                                                                aTGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATAAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGttgttg                                                             <  2:513913/133‑1 (MQ=255)
                                                                     cTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCCCGAGTTTTATGCTAAAGATTGCAAGTTGCTTTTAAAAAATAAGATAACTGATACATTATAGCTGTTCTTGAGGGt                                                        >  1:173212/1‑133 (MQ=255)
                                                                                       cGTTAATCTTGCCATGAATATTGCTTTCGTCGAAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTa                                      <  2:181749/133‑1 (MQ=255)
                                                                                                  ccATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGATGGTGCATGCTGCACAAAATTAAAGTTAAAAAg                           >  2:24087/1‑133 (MQ=255)
                                                                                                  ccATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAg                           >  1:435518/1‑133 (MQ=255)
                                                                                                        atatTGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGAGGGTGCATACTGCACAAAATAAAAGTTAAAGAGTAAAAc                     >  1:157949/1‑133 (MQ=255)
                                                                                                           ttGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGCGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAAcccc                  <  1:18717/133‑1 (MQ=255)
                                                                                                            tGCTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAAccccc                 >  2:369313/1‑133 (MQ=255)
                                                                                                             gCTTTCTTGGTAAAGTCTTGAGGTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCg                <  1:186598/133‑1 (MQ=255)
                                                                                                              cTTTCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGt               >  2:301447/1‑133 (MQ=255)
                                                                                                               tttCTTCGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGtt              <  1:301447/133‑1 (MQ=255)
                                                                                                                 tcttcGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGTTcc            >  2:462270/1‑133 (MQ=255)
                                                                                                                   ttcGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGTTCCtt          >  1:285865/1‑133 (MQ=255)
                                                                                                                    tcGTAAAGGCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAGTCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGTTCCTTa         <  2:497298/133‑1 (MQ=255)
                                                                                                                           ggCTCGAGTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGTTCCTTACCAGTTc  <  1:462270/133‑1 (MQ=255)
                                                                                                                                   |                                                                                                                             
AATGGGGGAATTCTGTGGCATGGTGTAAAGCACAGCAAAATCTTCAATAACGAAGCCAATTTTAATGTACTTACGAATTGGCGGTCACGTTAATCTTTCCATCAATATTGCTTTCTTCGTAAAGGCTCGAGTTTTTATGCTAAAGATTGCAAGTTGCTTGTAAAAGATAAGTACACTGATCCATAATCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGTTCCTTACCAGTTC  >  NC_000913/2700326‑2700582

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