Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA41.2% 17.5 / 17.2 34L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (14/6);  new base A (0/14);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 4.05e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCT  >  NC_000913/357742‑358032
                                                                                                                                                     |                                                                                                                                                
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                 <  1:635849/150‑1 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                               <  1:24643/151‑1 (MQ=255)
   aggCGTAGCGCATCAGGCAATTCCGGCCGCGGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                  <  2:362842/144‑1 (MQ=255)
       gTAGCGCATCAGGCAATTCCAGCGGCAGACCTGTGTCAGCGGCGACCGTAATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAGGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgcc                                                                                                                                        <  1:857634/151‑1 (MQ=255)
         aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑ccgccg                                                                                                                                       >  2:843443/1‑150 (MQ=255)
                cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTcc                                                                                                                               >  2:1146354/1‑151 (MQ=255)
                    cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCagt                                                                                                                            >  2:808767/1‑150 (MQ=255)
                        tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTgg                                                                                                                       >  2:755213/1‑151 (MQ=255)
                         ccgcccccggtcttgtcTCAGCGCCTCCCGTAATTCTTTCCCCCCAA‑CACCCGGTCTTTTACCCAACCCCATGTGAAAAAATGCCTGCGCGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTg                                                                                                                        <  1:298750/136‑1 (MQ=255)
                           aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                    >  1:252279/1‑151 (MQ=255)
                                ctgatcTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAA‑CAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTATTTTTTAGGCGCAACAAAA‑CCGCCGCTCCCTCCactagcgacatcg                                                                                                                  >  2:298741/6‑140 (MQ=11)
                                         gtCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑gaa                                                                                                        >  2:1158957/1‑148 (MQ=255)
                                            aGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAagcggaaga                                                                                                       >  2:604256/1‑140 (MQ=255)
                                                gCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑aaaaaaactt                                                                                                 >  1:59811/1‑141 (MQ=255)
                                                    cggtttttCTTTCCCCCCAA‑CACCCGCGTTTTTACCCACCCCAATGTAAAAAAAGGCCGCGGGGGGGGCGGTCCGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCg                                                                                          <  2:82945/145‑1 (MQ=255)
                                                       tGATTC‑TTTCCCCCAACCACCCCCGCTTTTCCCCACCGCCAGGTAAAAAAATCCTTTCGCGGCGGCTGCCGGCAACTTTTTTTCCCCCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAg                                                                                              <  1:223726/144‑1 (MQ=255)
                                                                  cGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGACAAG‑AAAAAACATTAAgtaccccctgtt                                                                                   >  2:801461/1‑133 (MQ=25)
                                                                         caacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCGCGAG‑AAAAAAAATTAAAAACCTCATCTTGTGTggggatggg                                                                      >  2:78091/1‑144 (MQ=25)
                                                                            ccGCGCTTTTCCCCAACCCATTGTAAAAAAATCCCTTGGCGGCGGCTGTCGGCCGGTTTTTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                                     <  2:927817/149‑1 (MQ=255)
                                                                                   tttACCCACCGCTATGTCAAAAAAGCCCTGCGCGGGGGCTGCTTGCGACTTTTTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                             <  2:356415/150‑1 (MQ=255)
                                                                                          aaCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCACAG‑AAATAAACTTAAACCCCCCCCGTTTTTTGGGCATTGCTGCCGGGAACAAta                                                        >  2:473687/1‑147 (MQ=21)
                                                                                                       aaacggctttCGCGGCGGTTTGCGGCCATGTTTTTTCCCCCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                        <  1:439165/143‑1 (MQ=255)
                                                                                                       aaaaTGCCTCGGGGGCGTCTGTCTGCAAGTTTTTTTCCCCCCGCAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                        <  2:43814/151‑1 (MQ=255)
                                                                                                              cTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACCTTTAAACCC‑CCCTTTTTTGTGGGAATAGCCGCTGGTGACAAAGTGGGAAGGGGAgggggggg                                     >  1:156939/1‑140 (MQ=255)
                                                                                                                    ggggctttccgcccAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAg                               <  2:181523/135‑1 (MQ=255)
                                                                                                                    ggcggcTTTCTGCAATTTTTTTTCCCGCCCAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                             <  1:160803/149‑1 (MQ=255)
                                                                                                                     ggggttGTCTGCCCGTTTTTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                                <  1:329074/140‑1 (MQ=255)
                                                                                                                     gcggtTGTTTGCCATTTATTTTGCCGCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc                           <  1:500030/150‑1 (MQ=255)
                                                                                                                                 atgtttttttcCCCCCACAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg                  <  2:982734/141‑1 (MQ=255)
                                                                                                                                   gtttttttcCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt                <  2:1121203/143‑1 (MQ=255)
                                                                                                                                     tATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTTAAACCCG‑CCTGTTTTGTGGGGATAGCCGTTGGGGACAAAGTGGGAAGGGGAGAGGGGCGGATGAGCTCCAGaaaccccccaatt            >  1:808056/1‑138 (MQ=255)
                                                                                                                                       tttttCCCCCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt             <  1:1055638/147‑1 (MQ=255)
                                                                                                                                        ttttCCCCCCCAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg         <  1:717254/150‑1 (MQ=255)
                                                                                                                                            ccgccaCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc      <  1:976489/148‑1 (MQ=255)
                                                                                                                                            cccccgaAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc   <  1:422268/148‑1 (MQ=255)
                                                                                                                                               ccgcAAAAACG‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt  <  2:252271/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCT  >  NC_000913/357742‑358032

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