Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA41.9% 41.9 / 21.1 44L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (16/9);  new base A (0/18);  total (16/28)
Fisher's exact test for biased strand distribution p-value = 1.12e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357742‑358028
                                                                                                                                                       |                                                                                                                                             
gATAAGCGTAGCGCCTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                             <  1:803661/151‑1 (MQ=255)
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                              <  1:223995/150‑1 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                            <  2:41241/151‑1 (MQ=255)
   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                           >  1:669994/1‑150 (MQ=255)
   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                         <  2:135736/151‑1 (MQ=255)
      cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                      <  2:970472/151‑1 (MQ=255)
        tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                     <  2:833948/150‑1 (MQ=255)
                    cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                         >  2:1033389/1‑150 (MQ=255)
                         ccAGCCGCAGCCCTGTGTCAGCGGCTACCGTGTTTCATTCCCGCCAACACCCGCGCATTTACCCACCCCCATGTGAAAAAATGCTTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                    <  1:165693/134‑1 (MQ=18)
                          cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                  >  2:113506/1‑151 (MQ=255)
                               gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                >  2:771791/1‑148 (MQ=255)
                                   aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggaga                                                                                                          >  2:889914/1‑145 (MQ=255)
                                   aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                               >  2:244458/1‑145 (MQ=255)
                                    ccTGTGTCGGGGGCTCCCGTGTTTCTTTTCCCCCAACACCCGCGTTTTTACCCAACCCCATGGGAAAAAATCCCTCCGGGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                           <  2:888673/148‑1 (MQ=255)
                                        tgtCAGCGGCTCCGGTGTTTCATTCCCCCAAACACCCGCGCTTTTCCCCACCGCCATGTAAAAAAATGCCTTCGGGGCGTCTGTCGCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                          <  2:880820/145‑1 (MQ=255)
                                         gtCGGCGGCTCCGGTAATTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCAATGTGAAAAAATCCCTCCGCGGGGGCTTTCTGCCATTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                           <  2:447368/143‑1 (MQ=255)
                                             gCGGCTCCCGTTTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCTTGTAAAAAAATGCCTCCGGGGCGGTTTTCTGCCATTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCtt                                                                                              <  2:708350/150‑1 (MQ=255)
                                                   aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGgaaaaaaatttatacc                                                                                        >  1:1319091/1‑138 (MQ=25)
                                                   aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGaaaaaaattttaaacc                                                                                        >  2:63188/1‑138 (MQ=25)
                                                       tatttCATTCCCCCCAACACCCGCGATTTTACCCACCGCCATGTCAAAAAAGGCCTTCGCGGCGGCTGTCTGCCAGGTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                        <  2:151544/143‑1 (MQ=255)
                                                           tCATTCCCCCCAACACCCACGAATTTCTCCACCCCCATGTACAAAAAGGCGTTCGCGGGGGTTGTCTGCACGCTTTTTTCCCGCCACAACAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa                                                                                   <  2:583011/147‑1 (MQ=255)
                                                             ttcccccccAACACCCGCGCTTTTACCCACCGCCTTGTCAAAAAATCCCTTCGGGCGGGTTTTCCGCAACCTATTTTCCCGCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                              <  2:15508/143‑1 (MQ=255)
                                                              tcccccccAACACCCGCCCTTTTACCAACCCCATTGTAAAAAAACCCTTCGGGGGGGGTTTTCTCCCAGTTATTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt                                                                             <  1:445121/144‑1 (MQ=255)
                                                                 ccccaaacccccGCGCATTTACCCACCGCCAGGAGAAAAAATCCCTCGGGGGCGGCTGTCTGCCGTCTTTTTCCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                          <  2:1275957/140‑1 (MQ=255)
                                                                 ccccaaacaccCGCTCTTTTACCCACCCCATTGTAAAAAAAGGCCTTCGGGGCGGTTTTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtg                                                                            <  2:126057/143‑1 (MQ=255)
                                                                  ccccaccaCCCGCGCTTTCATCCCCCGCTGTGTAAAAAAATTCCCGCGGGGCGGCGTCCACCCAGCTTTTTTCCCCCGGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                          <  2:1307605/147‑1 (MQ=255)
                                                                      aacaccCCCGCATTTACCCACCGCCATGTGAAAAATTGCCTTCGCGGCGTCTGTCTCCCATTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                    <  2:185310/151‑1 (MQ=255)
                                                                               gcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAAACCCTCCTGTTTTGTGGGAAACgc                                                                  >  2:610625/1‑144 (MQ=25)
                                                                                  tttccccACCGCCATGTGAAAAAATCCCTGCGGGGGGGCTTTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                         <  2:1225706/145‑1 (MQ=255)
                                                                                  tttACCCACCGCCATGTAAAAAAATGCCTCCGCGGGGCCTGCCGCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                          <  2:908709/149‑1 (MQ=255)
                                                                                        cAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGAGGCAGAGAAAAAAATTTTAGAACCCCACGTTGTGTGGGGATCGGCGCTGGtgaaaaagt                                                   >  2:1018711/1‑142 (MQ=21)
                                                                                         aaCGCCATGTAAAAAAATCCCTGCGGGCCGGTTGTCGCCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                  <  1:1165073/150‑1 (MQ=255)
                                                                                             ccATGTGAAAAAATCCCTCCGCGGCGGCTGTCTCCCAGCTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                            <  2:1279255/151‑1 (MQ=255)
                                                                                                  tgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAAAAAAAAATTTAAACCCGCCACGGTTGGTGGGGAACCCTGCTGGGGACAAAGAGGGAAggg                                           >  2:486236/1‑148 (MQ=255)
                                                                                                     aaaaaGCCCTTCGGGGCGGCTTTCTGCCTGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                         <  2:106927/147‑1 (MQ=255)
                                                                                                     aaaaaCGCCTCCGTGCGGGCCTTCTGCCATTATTTTTCCCCCAACAATAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                        <  2:356179/148‑1 (MQ=255)
                                                                                                           gCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAAATTTTAACCCCTCCTGTTTTGTGGGGAACGCTGCTTGGGACAAAGTGGGAAGGGGAgggggggg                                 >  2:205143/1‑142 (MQ=255)
                                                                                                              ttCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAAATTTAAACCCGCCAGTTTTGTGGGGAAAGCTTCTGGGGACAAAGTGGGAAGGGGAggggggggg                                >  2:420308/1‑139 (MQ=255)
                                                                                                                     gcggtTGTCTCCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                      <  2:998501/151‑1 (MQ=255)
                                                                                                                     gcggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAAGAAAAAAAATTTAACCCCCCCTGTTGTGTGGGGATCCCTGCTGGTGAAAAAGTGGGAAGGGGAGGGGGGGGGATCAGATcc                       >  2:1186128/1‑150 (MQ=255)
                                                                                                                       ggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAAAAAAAACATTAAACCCGCCCTGTTTTGTGGGGATCACCGCCTGGAACAAATTGGGAAGGGGAGGgcgggggaagagctgc                       >  2:224169/1‑136 (MQ=255)
                                                                                                                           gTCTGCAAGTTTTTTTCCCCCCCAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                 <  2:549864/150‑1 (MQ=255)
                                                                                                                             cTGCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                     <  2:1233299/144‑1 (MQ=255)
                                                                                                                               tccAGCTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                  <  1:200379/145‑1 (MQ=255)
                                                                                                                                 cAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAACATTTAAACCCCCCCGTTTTGTGGGAATCGCTGCCGGTGAAAAAGTGGAAAAGGGAGGGCGGGGGAACAACTcc                       >  2:394392/1‑138 (MQ=255)
                                                                                                                                           tcccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc  <  2:401488/148‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                             
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357742‑358028

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