Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8750AT40.5% 43.9 / 13.8 37N272K (AAT→AAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (17/5);  new base T (0/15);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 3.12e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357727‑358019
                                                                                                                                                       |                                                                                                                                                 
tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTAtt                                                                                                                                                 >  1:118940/1‑151 (MQ=255)
          ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                        <  1:243324/150‑1 (MQ=255)
            ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                      >  1:920564/1‑150 (MQ=255)
             ctaaaaaGCGTAGCGCATCAGGCTATTCCACCCGCAGTCCTGGCTCCGCGCCTTGCTTGATCCATTTCCGCCAACACCCGCGCATTGTTCCACCGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCTGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                      <  1:1172270/144‑1 (MQ=17)
               gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                   <  2:659290/150‑1 (MQ=255)
                aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GAGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaa                                                                                                                                     >  1:477322/1‑147 (MQ=255)
                  aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                               >  1:1100999/1‑151 (MQ=255)
                       tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                          >  2:31148/1‑151 (MQ=255)
                        aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                         >  1:919355/1‑151 (MQ=255)
                            cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                     >  1:725895/1‑151 (MQ=255)
                            cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                     >  2:783481/1‑151 (MQ=255)
                             aTCAGGCAATGCCAGCCGCAGACCTGAGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGACGTtct                                                                                                                      >  1:355597/1‑149 (MQ=255)
                                          aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gg                                                                                                      >  1:190301/1‑151 (MQ=255)
                                            ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                    >  2:814014/1‑151 (MQ=255)
                                                  ccgtgtgCCGGCGTCTCCGGTGATTTTTTCCCCCAAACACCCGCCCTTTTACCCAACCCCATG‑GCAAAAAGTCCCTC‑C‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAag                                                                                                <  2:1219728/146‑1 (MQ=255)
                                                   ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGCGGC‑ggggg                                                                                                   >  1:5610/1‑144 (MQ=255)
                                                     tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggggggacagg                                                                                             >  1:693371/1‑142 (MQ=255)
                                                             ctgCTACGGTATTTCTTTCCCGCCAACACCCCCGCTTTTACCCAACCCCTGGTG‑AAAAA‑AGCCTTCC‑GCGGGGGCTTTCTCCAGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAAt                                                                                          <  2:731707/142‑1 (MQ=255)
                                                                ctccggtgtTTCTTTCCCGCCAACACCCGCCCTTTTACCCACCGCCTTGTGAAAAAA‑TCCCTC‑C‑GGGGCGCTTGTCTCCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAg                                                                                 <  2:836702/145‑1 (MQ=255)
                                                                  aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGaaaaaaaattaaaagc                                                                              >  2:585065/1‑138 (MQ=25)
                                                                      attttCATCCCCCCCACCACCCGCCCTTTCACCCACCGCATTGTCAAAAAC‑GGCCTC‑C‑GGGGCGCCTGTCTCCCAGTTTTTTTCCCGCCACAAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGc                                                                                <  1:247189/142‑1 (MQ=255)
                                                                         taaatcccccccAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAA‑T‑CCTT‑CGGGGGGGGCTGTCTGCCAGTTTTTTTTCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAt                                                                         <  2:431267/139‑1 (MQ=255)
                                                                           cATTCCCCCCACCACCCGCGCATTTCTCCAACGCAATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTTTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt                                                                      <  1:1141385/150‑1 (MQ=255)
                                                                                       caacCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGAGAAATAAATTTTAACCCCGCCTGTTTTGTgg                                                                >  2:502570/1‑144 (MQ=25)
                                                                                                aTTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAGAAAAAAGCTTTAACACCG‑CCTGTTTTGTgggggagtct                                                        >  2:39833/1‑135 (MQ=25)
                                                                                                      caacccccATGTAAAAAAA‑TCCCTC‑C‑GGGGGGCCTGCCTCCCATTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                               <  1:918972/140‑1 (MQ=255)
                                                                                                            ccATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGAAAAAAATTTTAAATCCCCATGTTGTGTGGGAATAGCTGCTGGGaagcaagagggaa                                     >  2:491663/1‑139 (MQ=25)
                                                                                                            ccATGTAAAAAAA‑TCCTTC‑C‑GGGGCGGTTGTCTCCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                    <  2:625501/151‑1 (MQ=255)
                                                                                                                        attcccTC‑C‑GGGGCGGCGTTCGGCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                         <  1:308755/147‑1 (MQ=255)
                                                                                                                         acgCCTT‑C‑GCGGCGCGTGTCGGCCATCTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                        <  2:318772/149‑1 (MQ=255)
                                                                                                                                   ggggcggcTTTCCGCCGGCTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                        <  2:775044/141‑1 (MQ=255)
                                                                                                                                     ggcgtcTGTCCGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGaa                      <  1:1070772/143‑1 (MQ=255)
                                                                                                                                     ggcggtTTTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc               <  2:875240/150‑1 (MQ=255)
                                                                                                                                     ggcggtTGTCGCCCATTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                <  1:977023/149‑1 (MQ=255)
                                                                                                                                     ggcggcTGTCTGCCATTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc              <  2:919422/151‑1 (MQ=255)
                                                                                                                                     ggcggcTGTCTCCCTGTTTTTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                <  1:473264/149‑1 (MQ=255)
                                                                                                                                                  cAGTTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa  <  1:498929/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
TGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357727‑358019

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