Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,7350AT45.7% 33.5 / 21.6 35Y2* (TAT→TAAydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (2/17);  new base T (16/0);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 4.19e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATG  >  NC_000913/1413593‑1413884
                                                                                                                                              |                                                                                                                                                       
ttACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATg                                                                                                                                                   >  1:223230/1‑149 (MQ=255)
        tGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCCTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGATTTTTTTACAGGAAAAcccc                                                                                                                                           >  1:79228/1‑147 (MQ=255)
         gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACAGCTc                                                                                                                                           <  2:619254/148‑1 (MQ=255)
          ctctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAAAACCTCCctt                                                                                                                                       >  1:332365/1‑151 (MQ=255)
            ctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAAAACcccc                                                                                                                                          >  2:2504/1‑146 (MQ=255)
              ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattat                                                                                                                                     <  2:189022/149‑1 (MQ=255)
                ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATTAACACCTCCcttattttt                                                                                                                                 >  2:344352/1‑151 (MQ=255)
                ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattattat                                                                                                                                  <  2:352007/150‑1 (MQ=255)
                ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattatt                                                                                                                                    <  2:568539/148‑1 (MQ=255)
                    aCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAACACCTCCCTTAtttttt                                                                                                                                >  2:164843/1‑148 (MQ=255)
                             tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTc                                                                                                                     <  2:171974/149‑1 (MQ=255)
                             tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCg                                                                                                                    <  2:372212/150‑1 (MQ=255)
                                       cTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTTTATGAAAACCTCCCTTATTATTTTTCGTGGTTTGtttttttc                                                                                                            >  2:209728/1‑146 (MQ=255)
                                         gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGAAAACCCCCATTTTTTTTTTCAGTGGGTCGtttttt                                                                                                              >  2:32943/1‑144 (MQ=255)
                                             gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGa                                                                                                   <  1:499041/151‑1 (MQ=255)
                                                        cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAAAACCCCCATTATTTTTTCTCGTGGGTTGTTTTTTTCATCTTTTTTGTGCTtttttt                                                                                        >  2:567807/1‑151 (MQ=255)
                                                        cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCtttt                                                                                        <  1:24966/151‑1 (MQ=255)
                                                         aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTAATTAAAAACCCCCTTTTTTTTTACAATGGTGTCGTTTTTTCATCttttaaacgccttttttt                                                                                       >  2:232894/1‑136 (MQ=255)
                                                                  tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTCATTAAAACCCCCCCTATTTTTTTCCTTGGGTCGTTTTTTTCATTTTTTGAGTGCTTCtttttttctttccc                                                                              >  2:635736/1‑142 (MQ=255)
                                                                      tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGAACACCCCCCTTTTTTTTTTCCGTGGGTTGTTTTTTTCCTCTTTTGAGTGCTTTTTTTTTTCTTccacc                                                                            >  1:17784/1‑148 (MQ=255)
                                                                          ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGCGCCACCTTCCTTTTCAATAGTGGCGGTGAGTTTATTAATGGAGACCTCCCGTAAAATTTCCACCTGGTCGTTTTTTTCACCTTTCCTGTGCTTTTTTTTCACTTTCACCCTaaatc                                                                     >  2:474383/1‑148 (MQ=255)
                                                                           tAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCg                                                                     <  2:332271/150‑1 (MQ=255)
                                                                           tAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAATTTTTTTAATGAAAACCTCCCTTTTTTTTTTTAATGGGTCGTTTATTTCTCCTTTTTAGTGCTTTTTTTTTTCTTCCTCCATAACCgg                                                                    >  2:309095/1‑151 (MQ=255)
                                                                                    aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGAtt                                                           <  1:473563/151‑1 (MQ=255)
                                                                                          ttGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAACACCCCCATTATTTTTTTCAGTGGGTCGTTTTTTTCCTCTTTTCAGTGGTTCTTTTTTTCTTTCCCCATTACCGGGTCCGGGTTCGCCTc                                                     >  1:433664/1‑151 (MQ=255)
                                                                                          ttGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAATAAACACCTCCATTATTTTTTTCAGTGGGTTGTTTTTTTCATCTTTTTAGTGGTTCTTTTTTACTTTCCCCCAAAACCGGTTGGGGtttgcat                                                       >  1:429758/1‑146 (MQ=255)
                                                                                             tCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAAAACCCCCCTTTTTTTTTTCAATGGGTTGTTTTTTTCCTCTTTTGAGTGCTTTTTTTTTTCTTTCCCCCTAAACCGTTTGGGGTTCGCATcc                                                    >  1:7370/1‑149 (MQ=255)
                                                                                                   tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTacac                                             <  2:498574/150‑1 (MQ=255)
                                                                                                    agtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACtt                                           <  1:690690/151‑1 (MQ=255)
                                                                                                     gtatGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACCTCCACCATAGCCGGTTCGGGATTCGCATCCGTTACACTTc                                          <  2:238756/151‑1 (MQ=255)
                                                                                                             gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTa                                  <  2:493477/151‑1 (MQ=255)
                                                                                                                          tCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAAt                      <  1:543558/150‑1 (MQ=255)
                                                                                                                                         aTTTTATACATGAACACCGCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTATCGAGTGCTTCTTCTTCGCTTCCTTCATAACCGGTTCGGGATTCGCGTCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCTCCGTTa        <  2:496746/149‑1 (MQ=255)
                                                                                                                                           tttATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTg     <  1:124662/150‑1 (MQ=255)
                                                                                                                                            ttATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt   >  2:67019/1‑151 (MQ=255)
                                                                                                                                             tataCATGACCACCTCCATTATCATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATCACCGGTTCGTGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGa    >  2:444601/1‑149 (MQ=255)
                                                                                                                                              ataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg  <  1:606995/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                       
TTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATG  >  NC_000913/1413593‑1413884

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