Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913345,3280CA55.0% 23.0 / 66.4 40intergenic (+337/‑76)yahL/yahMuncharacterized protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (15/3);  new base A (0/22);  total (15/25)
Fisher's exact test for biased strand distribution p-value = 2.03e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-06

AATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACAT  >  NC_000913/345183‑345473
                                                                                                                                                    |                                                                                                                                                 
aaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCaaaaa                                                                                                                                              <  1:104949/151‑1 (MQ=255)
    aTCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCg                                                                                                                                           >  2:52784/1‑150 (MQ=255)
            aagttGTTTTGTGAATGTGTTGGTTGACATTCATATAAAAAAAATCATAATCCCATCATGTTTGTGTGGGGTTTTTTCTGTATCTTACGCACCGCACTCAACCCTGAAAGAAAATATGCTGTAAGGCTCATATAAAAAACCGCCATTAGCt                                                                                                                                  <  1:367349/147‑1 (MQ=255)
              gATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtca                                                                                                                                >  1:191993/1‑151 (MQ=255)
               atttttgGTGAATGTGTTGGTTGACATTCATATAAAAAAAATAATAATCCCATAATTTTTGTGTGGGGTCTTTTTTGTATCTTACCCACCGCACCAAACCCTGACAAAAAATATGCTTTAAGCCTCATATAAAAAACCGCCATTAGCt                                                                                                                                  <  1:377464/145‑1 (MQ=255)
                tGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtcat                                                                                                                               >  1:10363/1‑150 (MQ=255)
                    ttGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGa                                                                                                                          >  1:167292/1‑151 (MQ=255)
                       taaaTGTGTTGGTGGACTTTCATATGAAAAAAATAATAATCCCATCATTTTTGTGTGGGTTTTTTTCTGTATCTTACGCACCCCACTCAACCCTGACAAAAAATATGCTGAAAGGCTCATAAAAAAAACCGCCATTAGCTCATCAGGAag                                                                                                                        <  2:347195/148‑1 (MQ=255)
                          aTGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCa                                                                                                                    >  2:142940/1‑151 (MQ=255)
                              gttggttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGGCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacg                                                                                                                >  2:273992/1‑151 (MQ=255)
                                   tttaCATTAATATAAAAAAAACCAAAATCCCATCATGTTTGGGGGGGTTTTTTTTGTTATCTTCCGCATCCCACCCAGCCCTGAAAAAAAATATGCTTTAAGCCTCAAATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc                                                                                                           <  2:290831/148‑1 (MQ=255)
                                              ataaaaaaaaTCTAAATAACATTTTTTTGGTGTGTGTTTTTTTTTTTATCTTACTCATCCCACTCAACCCTAACAGAAAATTGGTTGGAAGGCTCAATTCAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACACCCAGGAttgttg                                                                                                <  2:187550/148‑1 (MQ=255)
                                              atGAAAAAAATCATAATCCCATCATGTTTGTGTGGGTTTTTTTCTGTATCTTCCGCACCGCACTCACGCCTACAAAAAAATATGCTGAAGGGCTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgttg                                                                                                <  1:340630/151‑1 (MQ=255)
                                              aaaaaaaaaaCCATAATCCCATCATGTTTGGGTGGGTTTTTTTCTTTATCTTACGCATCGCCCTCAAGCCTGAAAAAAAATATGCTGAAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgttg                                                                                                <  1:119056/148‑1 (MQ=255)
                                                 aaaaaaaTCAAATTCCCATCATGTTTGTGTGGGGTTTTTCCTGTATTTTACGCATCGAACACAACCTTAAAAAAAAATATGCTGAAGGCCTCATAAAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt                                                                                             <  1:16566/151‑1 (MQ=255)
                                                   aaaaataaaaatCCCATAATGTTTGTGGGGGGTTTTTTCTTTATTTTCCCCACCGCCCCCAAGCCTAAAAAAAAATATGCGGTAAGGCTCATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATgg                                                                                           <  2:227698/151‑1 (MQ=255)
                                                   aaaaaaaaaactcccATCTTTTTTGTGGGGGGTTTTTTCTGAATCTTCCCCACCCCCCCCAACCCAAAAAAAAAATAGGCGGTAAGACCAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATgg                                                                                           <  2:16683/138‑1 (MQ=255)
                                                   aaaaaTCATATTCCCAAGAGTTTTTTTGGGGTTCTTTCTGGGATCTTCCCCATCGCATCAAAGCCGAAAAAAAAATATGTTGTAGGGTTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATgg                                                                                           <  1:252942/151‑1 (MQ=255)
                                                   aaaaaTCAAAATCCCATCTTTTTTGTGGGGGTTTTTTTTTGTATCTTCCGCACCGCCCACAGGCCTAAAAAAAAAAAGGCTGAAAGGCTAATTTAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATgg                                                                                           <  2:260822/151‑1 (MQ=255)
                                                        tCAAAATCCCATATTTTTTGTTTGGGGTTTTTTCTTTATCTAACCCCTCGCCCTCAACCCGGAAAGAAAATATACTTTAAGCCTAATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgg                                                                                      <  2:333754/151‑1 (MQ=255)
                                                           aaaTCCCATATTGTTGTTGTGGTGTTTTTTTTGTTTCTCACGAATCACACTCAACCTTTAAAAAAAATTTGCTGAAAGCCTAATCAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGt                                                                                   <  2:354184/150‑1 (MQ=255)
                                                            aaTTCCAATCATTTTTGTGGGGGTTTTTTTTTTTTTTTTCCCCACCGCACCCAACCTTGACAAAAAATATGTTAAAGGCCTCAAATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGt                                                                                   <  1:174680/151‑1 (MQ=255)
                                                                ccATCATGTTTGTGTGGGTTTTTTTTTGTATCTTCCGCATCGCACTAAAGCCTGAAAAAAAATATGCTGAAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAg                                                                              <  1:190612/151‑1 (MQ=255)
                                                                      atGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCt                                                                          >  1:241815/1‑150 (MQ=255)
                                                                            gggtgGGTTTTTTTTTTTATCTTACGCACCGCCCCCAACCCGGAAAAAAAATATGCTGAAGGCCTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATgg                                                                    <  1:393098/148‑1 (MQ=255)
                                                                             tgtgGGGTTTTTTCTGTATCTTACGAACCGCACCAAGGCCTGACAAAAAATATGCTGAAGGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGcc                                                                  <  1:156281/151‑2 (MQ=255)
                                                                                 gggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc                                                               >  2:387782/1‑150 (MQ=255)
                                                                                   gTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCaa                                                            >  2:10924/1‑151 (MQ=255)
                                                                                       tttCTGTATCTTACTCATCGCACTCAAGCCTGACAGCAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGCTGCATGGAACGGGGTTCGAGGCCTCGATGGCGGTCTCACttt                                                        >  1:55437/1‑151 (MQ=255)
                                                                                       tttCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACttt                                                        >  1:160886/1‑151 (MQ=255)
                                                                                             gTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaa                                                   >  2:74160/1‑151 (MQ=255)
                                                                                                        caccccACTAAACCCTTAAAAAAAATATTCGGTAGGGCTAAAACAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                         <  2:286359/146‑1 (MQ=255)
                                                                                                          cccccctaaaGCCTGAAAGAAAATATGCTATAAGCCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAAt                                       <  2:326327/143‑1 (MQ=255)
                                                                                                          ccGCACACAAGCCTAACAAAAAATATGCAGAAAGGCTAAAATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAAt                                       <  1:182581/150‑1 (MQ=255)
                                                                                                             cACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAg                                    >  2:45716/1‑151 (MQ=255)
                                                                                                                          aaaaaaaTATCCTAAAAGCCTAATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTtattat                       <  2:179088/148‑1 (MQ=255)
                                                                                                                              aaaTATGCTGAAAGCCTAAAATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTatcatc                   <  2:84773/151‑1 (MQ=255)
                                                                                                                                     cTGAAAGGCCAAAATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCaa            <  2:78878/151‑1 (MQ=255)
                                                                                                                                           ggCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGa      >  2:124995/1‑151 (MQ=255)
                                                                                                                                               cATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACAt  >  2:145198/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACAT  >  NC_000913/345183‑345473

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