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

GAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357728‑358019
                                                                                                                                                    |                                                                                                                                                    
gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttt                                                                                                                                                   >  2:708987/1‑151 (MQ=255)
gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttt                                                                                                                                                   <  2:306798/151‑1 (MQ=255)
     cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                              <  2:815440/151‑1 (MQ=255)
      gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                             <  1:1081563/151‑1 (MQ=255)
         ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                                                                          >  2:1027438/1‑151 (MQ=255)
                tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                     <  1:832153/149‑1 (MQ=255)
                     gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                             <  2:262286/151‑1 (MQ=255)
                      tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                            >  1:1045077/1‑151 (MQ=255)
                       aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGAGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                           >  1:985238/1‑151 (MQ=255)
                            aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct                                                                                                                      >  1:1102146/1‑151 (MQ=255)
                               aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                    >  1:71640/1‑150 (MQ=255)
                                ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                    >  1:244591/1‑149 (MQ=255)
                                  cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                 >  1:200823/1‑150 (MQ=255)
                                       ccAGCCGCAGACTTGTGTCAGCGGCTCCCGTATTTCATTCCCGCAAACACCCGCGCATTTACCCACCGCAGTGTGCAAAAATGCTTCGGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagt                                                                                                                 <  2:529796/145‑1 (MQ=255)
                                           ccGTAGACCTGGGTCGGGGCCTCCCGTGTTTCTTTCCCCCCAACCCCCGCGTTTTTCCCCACCCCCATGTAAAAAATTCCCTCCGGGGCGGCTGTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                        <  1:714356/150‑1 (MQ=255)
                                           ccGCGGCCCTGTGTCAGCGGCTCCCGTATTTCATCCCCCCCACCACCCGCGCATTTACCCACCGCCTTGGCCAAAATTCCCTCCGGGGGGGTTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                         <  1:830316/149‑1 (MQ=255)
                                           ccGCAGCCCTGGTCCGGCGGCTACCGTATTTCTTTCCCGCCAACACCCGCGCTTTTACCCAACCCATTGTGAAAAAATCCCTTCGGGGTGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                        <  2:1062208/150‑1 (MQ=255)
                                           ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggtg                                                                                                       >  1:384474/1‑149 (MQ=255)
                                                    tgtgTCAGCGCTTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCAATGTGAAAAAAGGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                               <  1:761859/149‑1 (MQ=255)
                                                      tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggccaaaaaaa                                                                                             >  2:1063213/1‑140 (MQ=255)
                                                      tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaa                                                                                             >  1:401778/1‑149 (MQ=255)
                                                            cGGCTCCGGTTATTCATTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTCAAAAAATCCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                 <  1:576557/140‑1 (MQ=255)
                                                             gcctcccGTGTTTCATTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTCAAAAAATCCTTCCGCGGGGGCTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                                             <  2:1097357/137‑1 (MQ=255)
                                                                  ccgttttccATTCCCCCCACCACCCCCGCATTTACCAACCCCCATGTAAAAAAGTGCCTCCGCGGCGGCTGTCTCCCACTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                                  <  1:809948/140‑1 (MQ=255)
                                                                     tgtttaatccccGCCAACACCCGGCCTTTTACCCAACGCCTTGTGCAAAATTGCTTCCGCGGGGGCTGTCGGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                                                                              <  2:659308/139‑1 (MQ=255)
                                                                        ttCTTTCCCGCCAACACCCGCGCTTTTATCCAACGCCATGTGAAAAAATGCTTTCGCGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                          <  2:881029/150‑1 (MQ=255)
                                                                        ttCATTCCCCCCAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                          <  2:184412/150‑1 (MQ=255)
                                                                            ttccccccaccaCCCGCCCTTTTACCCACCGCAATGTAAAAAATTGCTTCCGGGGCGGTTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                         <  1:281628/141‑1 (MQ=255)
                                                                              ccccccaccaCCCGCCCATTTCCCCAACGCAATGTACAAAATCCCTTCCGCGGCGGCTGTCTCCCACTTCTTTTCCCGCCGCAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                          <  1:397683/140‑1 (MQ=255)
                                                                                       acccGCCCTTTTACCCACCGCCTTGTCAAAAAAGGCCTCCGCGGCGGTTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                            <  2:968768/147‑1 (MQ=255)
                                                                                                tttccccacccccATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                  <  2:181612/138‑1 (MQ=255)
                                                                                                               gtgAAAAAAGCCCTCCGCGGGGGCTTTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                 <  2:1045220/150‑1 (MQ=255)
                                                                                                               gtcaaaaaaTGCCTCGGGGGCGGTTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                  <  1:434130/145‑1 (MQ=255)
                                                                                                                 gCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTAAAACCCCCATGTTGTGTGGGCATAGCTGCTGGGAaaaaaatgggaa                                      >  2:773690/1‑136 (MQ=255)
                                                                                                                 acAAAAATCCCTGCGCGGCGTTTCTCCCCGAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                <  1:708247/148‑1 (MQ=255)
                                                                                                                   aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAGGCTTTTAACCCCCCCCTTTTTTTGGCAATTACTGTGGGGAAAAATGTGGggaagggaagg                              >  2:57288/1‑149 (MQ=255)
                                                                                                                   aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAAATTTAAACCCCCCTCTTTTGTGGGGATCGCTGCTGGGTAAAAAGTGGGAAAGGGagagg                             >  2:734337/1‑147 (MQ=255)
                                                                                                                   aaaaaTGCCTCGGCGGCGCCTTTCTGCCAGTTTTTTTCCCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                              <  2:33307/149‑1 (MQ=255)
                                                                                                                                  ggcggcTTTCTCCAATTTTTTTTCCCGCACAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                 <  1:1032896/148‑1 (MQ=255)
                                                                                                                                  ggcggcTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                <  2:1004887/149‑1 (MQ=255)
                                                                                                                                        tGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGCGAGAGAAAAAACTTTTAACCCCTCATTTTTTGTGGGGATCGCTGCTGGGGACAAAGTGGGAAGGGGAGGGGGGGGGAACAGCCCCCGaa          >  1:536154/1‑149 (MQ=255)
                                                                                                                                           cTGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa             <  1:204621/143‑1 (MQ=255)
                                                                                                                                              ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa  <  2:1098103/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357728‑358019

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