Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA47.6% 1.5 / 11.9 21L267F (TTG→TTTcynRDNA‑binding transcriptional dual regulator CynR
Reads supporting (aligned to +/- strand):  ref base C (10/1);  new base A (0/10);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 3.40e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCTTTCTCGCAGTAGTGGTCAATTTGCTCTCTGGTGGCAAATTCCGCGCTGAGCAGGACCAGTTTTTCAT  >  NC_000913/357645‑358128
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
aaaTGCCTTCGCAGCGGCTGGCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCCTCTCCGGCCTACGTAGAGCCCTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCCACGCCATGTGCCAAAATGCCCTCGCGGCGGCTGTCTGCCCGCTCTTTTTCCCCCGCGACAAAAc                                                                                                                                                                                                                                               >  2:15106/1‑251 (MQ=18)
      cTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                                                                                                                                        >  1:500813/1‑251 (MQ=255)
                gCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCCTTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                                                                                                                              >  1:45132/1‑251 (MQ=255)
                              gTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCCGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGCTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTCGTGGCgggggggg                                                                                                                                                                                                                >  1:456694/1‑247 (MQ=255)
                                aGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCCGTAGTGGCggggggtagg                                                                                                                                                                                                              >  2:172958/1‑245 (MQ=255)
                                    tATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggaggagaa                                                                                                                                                                                                           >  1:126954/1‑241 (MQ=255)
                                          ggATGCGGGGTGACCGCTTTACCCGGCCTACGTAGAGCACGAACCCGTTGGGCATAAAAGGCGTGGCCCATCGGCAAATCCCCGCCGCAGACTTGTGTCGGCGGCTCCGTTGTTTATTTCCCGCAAACACCCGCGCATTTACCCAACGCAATTTGAAAAATGCCCTCCGCGGCGGCTGTCTGCCGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                                                                                                                                                   <  1:265889/251‑1 (MQ=255)
                                                   gTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaattttaacccctc                                                                                                                                                                                           >  1:84037/1‑235 (MQ=255)
                                                   gTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAGGAGaaaaagattttaatccagc                                                                                                                                                                                           >  1:154796/1‑235 (MQ=255)
                                                                                                                     gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAACTTTAAAACCCTCCTGGTTTGTGGGGATCGCTGCTGGTGACCAAGTGGGAAGGGGAGAgcgggggaagagcgccaca                                                                                                                        >  1:592308/1‑235 (MQ=255)
                                                                                                                                 cgaagtccTGTGTCGGCGGCTACCTGATTTCATTCCCACCAACAACCGCGCATTTACCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt                                                                                                            <  1:383058/245‑1 (MQ=255)
                                                                                                                                                               tatttcccGCCAACACCCGCGCATTTATCCACCGCCATGTGCAAAATCCCCTGCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGAt                                                                                <  1:15106/244‑1 (MQ=255)
                                                                                                                                                                                                 cccttgtgAAAAAATCCCTCCGCGGCGGCTGTCGGCCACTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCTTTCTCGCAGTAGTGGTCAATTTg                                              <  1:381084/247‑1 (MQ=255)
                                                                                                                                                                                                    aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAGAAAAAAATTTTAACCCCGCCTTTTGTTTGGGGAAACCTGCTGGTGACAATGTGGGAAAGGGAGGGCGGCGGATCAGCTCCCGAACCCCGCCCATTTAGTTTGCCCCCAAGACCACCCGTGGGTGTAGCCCCCCCTtctccccgcagtggggaattttttt                                           >  1:168620/1‑230 (MQ=255)
                                                                                                                                                                                                        acaaaactcccTTCGCGGCGGTTTTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCTTTCTCGCAGTAGTGGTCAATTTGctctct                                        <  1:36370/241‑1 (MQ=255)
                                                                                                                                                                                                                     cgggggcgttttcctccacttttttttcCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCTTTCTCGCAGTAGTGGTCAATTTGCTCTCTGGTGGCAAATTcc                           <  2:125422/228‑1 (MQ=255)
                                                                                                                                                                                                                                  cgccatttttttttcCCGCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCTTTCTCGCAGTAGTGGTCAATTTGCTCTCTGGTGGCAAATTCCGCGCTGAGCAgg               <  2:168620/240‑1 (MQ=255)
                                                                                                                                                                                                                                   cgCCAGTTATTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCTTTCTCGCAGTAGTGGTCAATTTGCTCTCTGGTGGCAAATTCCGCGCTGAGCAGGAc             <  1:461306/250‑1 (MQ=255)
                                                                                                                                                                                                                                           ttttttCCGGCGGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCTTTCTCGCAGTAGTGGTCAATTTGCTCTCTGGTGGCAAATTCCGCGCTGAGCAGGACCAGttt       <  2:375209/249‑1 (MQ=255)
                                                                                                                                                                                                                                            tttttccccccgaaacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTTGATGTAGCCCCGCTTTCTCGCAGTAGTGGTCAATTTGCTCTCTGGTGGCAAATTCCGCGCTGAGCAGGACCAGtt        <  2:537488/235‑1 (MQ=255)
                                                                                                                                                                                                                                             ttttccccccgcAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCTTTCTCGCAGTAGTGGTCAATTTGCTCTCTGGTGGCAAATTCCGCGCTGAGCAGGACCAGttttt     <  1:560608/242‑1 (MQ=255)
                                                                                                                                                                                                                                               tcccgccgCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCTTTCTCGCAGTAGTGGTCAATTTGCTCTCTGGTGGCAAATTCCGCGCTGAGCAGGACCAGTTTTtcat  <  1:172958/248‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
AAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCTTTCTCGCAGTAGTGGTCAATTTGCTCTCTGGTGGCAAATTCCGCGCTGAGCAGGACCAGTTTTTCAT  >  NC_000913/357645‑358128

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: