Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC50.0% 32.0 / 10.3 36K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (13/5);  new base C (0/18);  total (13/23)
Fisher's exact test for biased strand distribution p-value = 7.42e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357730‑358027
                                                                                                                                                      |                                                                                                                                                      
aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                        >  1:246111/1‑151 (MQ=255)
    gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                                    <  2:167985/151‑1 (MQ=255)
          cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                              >  2:422950/1‑151 (MQ=255)
               aaGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTATTTTTTCGCCGGAACAAAAc                                                                                                                                         >  2:142740/1‑151 (MQ=11)
               aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                          <  1:823128/150‑1 (MQ=255)
                     aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                                   >  2:611468/1‑151 (MQ=255)
                        gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCCGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                                >  2:301362/1‑151 (MQ=255)
                                 aaTTCCAGCCGCAGCCCTGTGTCAGCGGCTCCCGTGATTCATTTCCGCCAACACCCGCGCATTTACCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCagta                                                                                                                        <  1:78712/150‑1 (MQ=255)
                                   ttCCAGCCGCAGCCCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACACCCGCGCTTTTATCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCagta                                                                                                                        <  1:955276/148‑1 (MQ=21)
                                      cacccGCAGACCTGTGTAGGCGGCCCCCGTTATCCATTCCCCCCACCACCCGCCCTTTTACCCACCCCCTTGTAAAAAAATCCCTTCGCGGCGGCTTTCGCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagta                                                                                                                        <  2:713824/142‑1 (MQ=255)
                                         ccggggtcatgtgTCAGCGGCTCCGTTGATTCATTCCCCCCAACAACCGCGAATTTATCAACCTCATTGTGCAAAACTGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCCCCGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                               <  1:852602/142‑1 (MQ=255)
                                            cAGACCTTTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                               >  2:271621/1‑148 (MQ=255)
                                             aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                               >  2:357204/1‑147 (MQ=255)
                                               aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAACGCCGTTCTCTCCAGTAGTGGCGGGGGaa                                                                                                            >  1:311339/1‑148 (MQ=255)
                                                ccggtgtCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCTTTTATCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                          <  1:444797/146‑1 (MQ=21)
                                                       cAGCGGCCCCCGTGTTTCATTCCCCCCAACACCCGCGCTTTTACCCAACGCCATGTAAAAAAATGCCTTCGCGGCGTCTTTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAga‑g‑a                                                                                                      <  2:892836/144‑1 (MQ=255)
                                                           ggCTCCCGTGATTCATTCCCGCCAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTaa                                                                                           <  1:375060/151‑1 (MQ=25)
                                                               cccGTTATTCATCACCACCAACACCCGGCCTTTTAACAACCCCCAGGTGAAAAACGGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTCCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGccc                                                                                       <  2:633727/150‑1 (MQ=255)
                                                                   gtttcctttcccGCAAACACCCCCTCTTTTCTCCACCCCCGTGTAAAAAAGTGCCTGCGGGGCGGCTGTCGGCGAGTTTTTTTCCCCCGGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATa                                                                                                  <  2:472094/126‑1 (MQ=255)
                                                                             ccGCCAACAACCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                          <  1:197103/150‑1 (MQ=255)
                                                                                  aacaccCGCGCTTTTACCCACCCCAATGTGAAAAAATCCCTCCGCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt                                                                     <  2:441608/150‑1 (MQ=255)
                                                                                     acccGCGCTTTTACCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                                   <  2:92739/147‑1 (MQ=38)
                                                                                     aacCGCGCTTTTACCCAACCCGAGGAACAAAATGGCCGCGGCGGCGGTTGGCGCCGGGTTATTTTCCCCCCGGAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                                 <  2:691809/151‑1 (MQ=255)
                                                                                         gcgcATTTACCCACCGCAATGTGCAAAAAGGCTTGCGGGGCGGCTGCCTCCAGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                             <  1:290140/151‑1 (MQ=255)
                                                                                               ttATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAAGTTTTAAACCC‑CCCTGTTTTGTGGGGAACGCTGCTGGTAAc                                                        >  1:236948/1‑150 (MQ=25)
                                                                                               ttACCCACCGCCATGTGAAAAAATGCTTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                        <  1:788095/150‑1 (MQ=255)
                                                                                                  cccACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTCCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                     <  1:189747/149‑1 (MQ=255)
                                                                                                             gtgAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAg                                            <  1:142748/147‑1 (MQ=255)
                                                                                                             gtcaaaaaaTGCCTCCGGGGCGGCTGTCGCCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAgg                                           <  2:440742/144‑1 (MQ=255)
                                                                                                                cAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑A‑AAAAAAATTTTAACCCCCCACGTTGTGGGGGGATCGGTGCCGGTAAAAATGA‑GGAAAGGGGAggg                                     >  1:284186/1‑149 (MQ=255)
                                                                                                                 aaaaaTCCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTgc                                    <  1:632022/151‑1 (MQ=255)
                                                                                                                   aaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAAAAAATTTTAACCCGCCCTGTTGTGTGGGAAACGCTGCTGGTGAAAAAGT‑GGGGAGGGAAgggggg                                  >  1:980127/1‑146 (MQ=255)
                                                                                                                          ttCGCGCCGCCTGCCCCCCAGTTTTTTTCCCCCCACAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATc                            <  2:456804/150‑1 (MQ=255)
                                                                                                                                   gcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAAAAAATTTTAGCCCCGCCTGTTGTGTGGGGAACGCTGCCGGTGACAAAGG‑GGGAAGGGGAGGGGGGGGGATAAGCTCCAGa                   >  1:778621/1‑150 (MQ=255)
                                                                                                                                           ccAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGCGAGG‑GTAAAAAGATTTAA‑CCCCTCATGTTTTTTGGTGAATACTTCTGGTGAAAATGT‑GGGAAGAGGAGTgaggtggaataagtcccg                    >  1:482377/1‑127 (MQ=255)
                                                                                                                                           ccAACTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa          <  1:232567/151‑1 (MQ=255)
                                                                                                                                               cTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAAAAACTTTAAACCCCTCCTTTTGTTTGGCGAACGCTGCTGGGGAAAAAAT‑GGGAAAGGGAGGGGGGGGGATCAGCCCCAGAAACGCCCCAAtt       >  1:111689/1‑150 (MQ=255)
                                                                                                                                                 tttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt       <  1:109716/147‑1 (MQ=255)
                                                                                                                                                  ttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg    <  2:122914/150‑1 (MQ=255)
                                                                                                                                                     ttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGCAAAA‑GAAAAAAATTTTAAACCCG‑CCCGTTTTGTGGGGATAGCTGGTGGGGACAATGT‑GGGAAGGGGAGGGGGGGGGATCAACTCCCGAAACCCGCTAATTTAggt  >  2:839844/1‑147 (MQ=255)
                                                                                                                                                     ttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGG‑AAAAAAATTTAA‑CCCCCCCTGTTTTGTGGGGAACGCCGCCGGTGACAATGTGGGAAAGGGAGGGGGGGGGAA‑CAGCCCCCGAACCccgcccaattaggt  >  2:736833/1‑139 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357730‑358027

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