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

CTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCG  >  NC_000913/357726‑358014
                                                                                                                                                     |                                                                                                                                           
cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAt                                                                                                                                            >  2:216725/1‑151 (MQ=255)
   aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAtttt                                                                                                                                         >  2:197877/1‑151 (MQ=255)
     cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                                       >  1:182952/1‑151 (MQ=255)
        gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                    >  2:760505/1‑151 (MQ=255)
             ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                >  2:940721/1‑150 (MQ=255)
               tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                              >  2:928073/1‑150 (MQ=255)
                gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                             <  1:824979/150‑1 (MQ=255)
                gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                              <  1:345863/149‑1 (MQ=255)
                   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                         <  1:18606/151‑1 (MQ=255)
                    aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                        <  2:244464/151‑1 (MQ=255)
                     gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCCGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                        >  1:942168/1‑150 (MQ=255)
                       gTAGCGCATCAGGCAATTCCTGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                     >  1:889348/1‑151 (MQ=255)
                        tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                    >  2:55688/1‑151 (MQ=255)
                                     aaTTCCAGCCGCAGAGCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAACGCCGTTCTCTCCagt                                                                                                         >  1:119241/1‑149 (MQ=255)
                                        tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                    >  1:833336/1‑151 (MQ=255)
                                         ccAGCCGCAGCCCTGTGTCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCATTTACCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTTTCTCCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                     <  2:12790/149‑1 (MQ=21)
                                           aGCCGCAGACCTGTGTCAGCGGCTACTGTGGTTCATTTCCGCCAACAACCGCGCATTTATCAAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGGCGCAACAAAACCGCCGTGCTCTCGAGTAGTGGCgg                                                                                                 >  1:808498/1‑151 (MQ=21)
                                            gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                 >  1:475324/1‑150 (MQ=255)
                                                   cccTGTGCCGGCGCCTACCGTGTTTCATTCCCCCCAACACCCGCCCTTTTATCCACCCCCATGTGCAAAAATCCCTCCGCGGCGGCTCTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                         <  1:427372/150‑1 (MQ=255)
                                                        tgccagggacTACGGTGATACTTTCCCGCAAAAACCCGCGCTTTTATCCAACGGCATGTCCAAAATTGCCTTCGGGGCGGCTGTCCGCCGGTTATTTTCCCGCCGAAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                     <  2:363467/141‑1 (MQ=255)
                                                             gcgcCTACGGTGATTCTTTCCCGCCAACACCCGCGCATTTACCCACCGCCTTGTGAAAAAATGCTTCCGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                               <  2:584613/151‑1 (MQ=255)
                                                               gtcTACCGTGATTCTTTCCCCCCAACACCCGCGCATTTACCCAACGCCTTGTGCAAAAATCCCTTCGCGGCGGCTTTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                                      <  2:713037/140‑1 (MQ=255)
                                                                cctcccttaaTTCTTTCCCCCCAACACCCGCGCTTTTCCCCACCGCCATGTGAAAAAATCCCTCCGCGGCGCCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                               <  1:726809/139‑1 (MQ=255)
                                                                 cTACGGTTATTCTTTCCCACAAACACCCGCACTTTTACCAAACCCCTGGTAAAAAAATCCTTCCGCGGCGGTTGTCGCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                            <  1:850144/150‑1 (MQ=255)
                                                                  tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGgaaaaaattttaaac                                                                          >  2:963721/1‑139 (MQ=25)
                                                                    ccgttttcaaTTCCCCCCACCACCCGCACATTTACCCACCGCAATGTCCAAAATTCCCTCCGCGGCGGCTTTCTCCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                            <  2:470988/138‑1 (MQ=255)
                                                                                    ccAACACCCGCGCTTTTACCCACCGCCTTGTACAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                         <  2:891928/150‑1 (MQ=255)
                                                                                           ccGCGCTTTTACCCACCCCCGTTGCCAAAAAGGCTTGCGCGGGGGCTTCCTCCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCATCATGTTGTGTGGCAATc                                                      <  1:749897/146‑1 (MQ=255)
                                                                                            cgcgCATTTACCCACCGCCATGTGCAAAAAGGCCTCCGCGGCGGTTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                <  2:343914/151‑1 (MQ=255)
                                                                                             gcgcTTTTACCCACCGCAATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                               <  2:736152/151‑1 (MQ=255)
                                                                                              cgcTTTTACCCAACGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                               <  2:430683/150‑1 (MQ=255)
                                                                                                       ccacccccATGTGAAAAAATCCCTCCGCGGCGGCTGTCGGCAAGTTATTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                       <  2:973080/143‑1 (MQ=255)
                                                                                                         accGCCATGTCAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                    <  2:252399/148‑1 (MQ=255)
                                                                                                                     aaaaaTGCTTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                       <  2:554164/151‑1 (MQ=255)
                                                                                                                      aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTTAAACCCTCCTGTTGTTTGGGGATCGCCGCTGGTAAAAATATGGGAAGGGGAggggg                      >  2:1002766/1‑147 (MQ=255)
                                                                                                                       aaaTGCCTCCGCGGCGGCTTTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                     <  2:593834/151‑1 (MQ=255)
                                                                                                                       aaaGGCTTTCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                     <  2:102866/151‑1 (MQ=255)
                                                                                                                        aaTGCCTCCGGGGCGGCTTTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                    <  2:132330/151‑1 (MQ=255)
                                                                                                                                   gggggcTGTCTCCAATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt           <  1:835340/146‑1 (MQ=255)
                                                                                                                                   ggcggtTGTCTCCCGGCTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc            <  2:455541/148‑1 (MQ=255)
                                                                                                                                   ggcggcTGTCTGCCGGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc         <  2:980493/151‑1 (MQ=255)
                                                                                                                                               ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg  <  1:934882/146‑1 (MQ=255)
                                                                                                                                                cAGTTATTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc   <  2:511472/144‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
CTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCG  >  NC_000913/357726‑358014

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