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

CTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACC  >  NC_000913/357726‑358013
                                                                                                                                                      |                                                                                                                                           
cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAt                                                                                                                                            >  2:764309/1‑151 (MQ=255)
          aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                   <  2:648091/150‑1 (MQ=255)
             ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaa                                                                                                                               <  2:39649/151‑1 (MQ=255)
             ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                >  1:286530/1‑150 (MQ=255)
                   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                         <  1:565513/151‑1 (MQ=255)
                       gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                     >  1:798586/1‑151 (MQ=255)
                         aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                   >  2:774691/1‑151 (MQ=255)
                            gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATACCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                 >  1:25661/1‑150 (MQ=255)
                               tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                                                                                                             >  2:798994/1‑151 (MQ=255)
                                 aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                            >  2:71558/1‑150 (MQ=255)
                                         ccAGCCGCAGCCCTGTGTCAGCGGCTCCCGTGTTTCATTCCCCCCAACACCCGCGCATTTATCCACCGCCATGTGAAAAAATGCTTCCGCGGCGGCTTTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                   <  2:669285/151‑1 (MQ=255)
                                         ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                    >  1:653304/1‑150 (MQ=255)
                                           aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                 >  1:133839/1‑151 (MQ=255)
                                              cccagccctttccacgcgcCTACCGTAGTTCTTTCCCCCCAACACCCGCACATTTCCCCACCGCCTTGGAAAAAAATGCCTCCGCGGCGGTTGCCGCCCGGCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                              <  1:42657/136‑1 (MQ=255)
                                                       gtgtCAGCGCCTCCCGTGATTCATCCCCCCCAACACCCGCCCTTTTACCCAACGCCATGTCCAAAAATGCCTGCGGGGCGGGTGTCTGCCAGCTATTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                                      <  1:233055/150‑1 (MQ=255)
                                                         gtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggggggggaaaa                                                                                    >  1:14996/1‑139 (MQ=255)
                                                           cAGCGGCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCGTGTAAAAAATGGCTTGCGCGGGGCCTGTCTGCCATCTTTTTTCCCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                                                                    <  2:680066/148‑1 (MQ=255)
                                                            aGCGCCTCCCGTGATTCTTTCCCCCCAACCCCCGCGCTTTAACCCACCGCCTTGTGAAAAAAGGCCTGCGCGGCGGCTGTCGGCCAGCTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                   <  1:256580/148‑1 (MQ=255)
                                                              cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGCGaaaaaaa                                                                                  >  1:742869/1‑143 (MQ=255)
                                                                ccTACCGTGATTCTTCCCCGCAAACACCCGCCCTTTTCCCCAACGCCATGTAAAAAAATCCCTCCGCGGCGCTTTCCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                   <  1:84641/143‑1 (MQ=255)
                                                                    ccGTGATTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCATGTAAAAAAATCCCTGCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                       <  1:160788/151‑1 (MQ=255)
                                                                      gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGGaaaaaaaaataagcaaac                                                                       >  1:713579/1‑135 (MQ=25)
                                                                          ttCTTTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAAGGCCTTCGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                  <  2:151122/150‑1 (MQ=255)
                                                                          ttCTTTCCCGCCAACACCCGCCCTTTTACCCACCCCCATGTGAAAAAATGCCTTCGCGGCGGTTGTCCGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                  <  1:729441/150‑1 (MQ=255)
                                                                          ttCATTCCCGCAACCACCCGCGCTTTAACCCACCCCCATTTGCAAAAATCCTTCCGCGGCGGCTGTCTGCCACCTTTTTTGCCGCGGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                  <  2:443619/150‑1 (MQ=255)
                                                                           tCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAgaatttaaccc                                                                         >  1:573985/1‑133 (MQ=25)
                                                                                          cccGCGCTTTTACCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGTTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                 <  2:682027/150‑1 (MQ=38)
                                                                                                 aTTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAAAAAAAACATTTAAACCCCCATGTTTTGTGGGGATCGCCGCTgg                                             >  2:487542/1‑148 (MQ=25)
                                                                                                      cccAACCCCATGTCAAAAAATCCCTTCGCGGCGGCTGCCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa                                       <  1:562401/148‑1 (MQ=255)
                                                                                                         aaCCCCATGTACAAAAATCCCTCCGCGGGGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                     <  1:458175/148‑1 (MQ=255)
                                                                                                          aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAACTTAATCCCGTCATGGTGTGTGGGAATAGCTggtggtgaaaa                                       >  2:741978/1‑140 (MQ=37)
                                                                                                                 gtaaaaaaaTCCCTCCGCGGCGGTTGTCTGCCAGCTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa                                <  1:50233/141‑1 (MQ=255)
                                                                                                                 gggCAAAAAGGCTTGCGCGGGGGCTGTCGGCCAGCTATTTTCCCGCCACAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                           <  1:409403/148‑1 (MQ=255)
                                                                                                                  tgCAAAAATGCCTCCGGGGGGTCTGTCTGCCAGTTATTTTCCCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                            <  1:990987/148‑1 (MQ=255)
                                                                                                                     aaaaaTGCCTCGGCGGGGGCTGCCGGCAATTTTTTTTCCCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                           <  1:803090/146‑1 (MQ=255)
                                                                                                                         aTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAACATTTAACCCCCCCTGGTTTTGTGGGAACGCTGCTGGGAAAAAAGTGGGAAGGGGAgggggg                    >  1:478774/1‑144 (MQ=255)
                                                                                                                             cctgcccggcgtCTGTCGCCCGGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACATTGTGGAAAGGGAAGTGCGGCGaa                <  1:834864/144‑1 (MQ=255)
                                                                                                                                    ggcggtTGTCTGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc        <  2:356261/151‑1 (MQ=255)
                                                                                                                                    ggcgccTGTCGGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc         <  1:346635/150‑1 (MQ=255)
                                                                                                                                           gTCGGCCTGTTATTTTCCCGCCGAAACAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc  <  1:403072/150‑1 (MQ=255)
                                                                                                                                                ccAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg      <  1:235537/141‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                           
CTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACC  >  NC_000913/357726‑358013

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