Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,298,5140A.37.9% 64.2 / 40.6 29coding (1/54 nt)ychTprotein YchT
Reads supporting (aligned to +/- strand):  ref base A (6/12);  new base . (7/4);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACATTAGAGTTATTATCTCTAATGCGTCACTTCCAGGTGGCGTAAGCAAGATTACTCACTTCTGGGTACTGATTACGTGATCCAAATCAAATTTTTGCAAAGCTGACACCTTTCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTACTTCAAATTTTTCATCGGGTTATTTGCGCTGGTCAACCCGGTAGGGATTATTCCCGTCTTTATCAGCATGACCAGTTATCAGACAGCGGCAGCGCGAAACAAAACTAACCTTACAGCCAACCTGTCTGTGGCCAT  >  NC_000913/1298268‑1298764
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
aaCATTAGAGTTATTATCTCTAATGCGTCACTTCCAGGTGGCGTAAGCAAGATTACTCACTTCTGGGTACTGATTACGTGATCCAAATCAAATTTTTGCAAAGCTGACACCTTTCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACGAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAA‑TGCTg                                                                                                                                                                                                                                                       <  2:375123/251‑1 (MQ=255)
      agagTTATTATCTCTAATGCGTCACTTCCAGGTGGCGTAAGCAAGATTACTCACTTCTGGGTACTGATTACGTGATCCAAATCAAATTTTTGCAAAGCTGACACCTTTCAGCATCGCTTTTCGCCATTATAGCGAACAGTTAATAAATTGTAGTATGATGTGGTGGCTACATTAGCATGTTTTGCGCAACTAGATAACAATAACGAAGGATAGCAATTTTAAGGAGTTAGGAGGTGAAAAATGCTGTCaa                                                                                                                                                                                                                                                   <  1:171801/250‑1 (MQ=255)
                   ctAATGCGTCACTTCCAGGTGGCGTAAGCAAGATTACTCACTTCTGGGTACTGATTACGTGATCCAAGTCAAATTTTTGCAAAGCTGACACCTTTCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTAtt                                                                                                                                                                                                                                         <  1:395324/247‑1 (MQ=255)
                      aTGCGTCACTTCCAGGTGGCGTAAGCAAGCTTACTCACTTCTCGGTAATGATTACGTGATCCAAATCAAATTTTTGCACAGCTGACACCTTTCAGCATCCCTTTTCGCCATTATAGCTAACAGTTAATAACTTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTCTTGTCAgcg                                                                                                                                                                                                                                  >  2:323860/1‑251 (MQ=255)
                                tCCAGGTGGCGTAAGCAAGATTACTCACTTCTGGGTACTGATTACGTGATCCAAATCAAATTTTTGCAAAGCTGACACCTTTCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTc                                                                                                                                                                                                                         >  1:1621/1‑250 (MQ=255)
                                  cAGGTGGCGTAAGCAAGATTACTCACTTCTGGGTACTGATTACGTGATCCAAATCAAATTTTTGCAAAGCTGACACCTTTCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAAc                                                                                                                                                                                                                      <  1:126299/251‑1 (MQ=255)
                                        gCGTAAGCAAGATTACTCACTTCTGGGTACTGATTACGTGATCCAAATCAAATTTTTGCAAAGCTGACACCTTTCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAA‑TGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCttattt                                                                                                                                                                                                               <  1:245305/251‑1 (MQ=255)
                                                             tCTGGGTACTGATTACGTGATCCAAATCAAATTTTTGCAAAGCTGACACCTTTCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCgg                                                                                                                                                                                             <  1:120950/249‑1 (MQ=255)
                                                               tGGGTACTGATTACGTGATCCAAATCAAATTTTTGCAAAGCTGACACCTTTCAGCATCGCTTTTCGCCATTCTAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAA‑TGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCtttt                                                                                                                                                                                        >  2:465926/1‑251 (MQ=255)
                                                                        aTTACGTGATCCAAATCAAATTTTTGCAAAGCTGACACCTTTCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCgg                                                                                                                                                                                 >  1:24906/1‑250 (MQ=255)
                                                                                        caaatTTTTGCAAAGCTGACACCTTTCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCGGTCAAAAGGCGTAGTGTCAGCGCGTCTTTTCAGCCTTATTGATGGCTAACAGTATCCGGCGTTGGCTACGGGGCTaaacacgattcag                                                                                                                                                                <  2:57289/251‑12 (MQ=255)
                                                                                         aaatTTTTGCAAAGCTGACACCTTTCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAA‑TGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGAc                                                                                                                                                              <  1:465926/251‑1 (MQ=255)
                                                                                                   aaaGCTGACACCTTTCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGa                                                                                                                                                     <  2:339590/251‑1 (MQ=255)
                                                                                                              cTTTCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGAAATTTTAAGTAGTTAGGAGGTGAAAA‑TGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCACCCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGttt                                                                                                                                          >  1:26680/1‑250 (MQ=255)
                                                                                                               tttCAGCATCGCTTTTAGCCCTTATAGCTAACAGTTAATAAATTGTAGTATGCTTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGttt                                                                                                                                          >  1:399758/1‑250 (MQ=255)
                                                                                                                 tCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTAc                                                                                                                                        <  1:328755/250‑1 (MQ=255)
                                                                                                                         gCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTACTTCAAAttt                                                                                                                               >  2:337601/1‑251 (MQ=255)
                                                                                                                           ttttCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGCGGTGAAAA‑TGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTACTTCAAATTTTTc                                                                                                                            >  2:346124/1‑251 (MQ=255)
                                                                                                                                     tataGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTACTTCAAATTTTTCATCGGGtt                                                                                                                    <  1:337601/250‑1 (MQ=255)
                                                                                                                                                 ttAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAA‑TGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTACTTCAAATTTTTCATCGGGTTATTTGCGCTGGTc                                                                                                       >  1:449870/1‑250 (MQ=255)
                                                                                                                                                    ataaATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTGAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCGTGTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTGGATTTTCCCGTTTACTTCAAATTTTTCATCGGGTTATTTGCGCTGGTCAAc                                                                                                    <  1:234215/251‑1 (MQ=255)
                                                                                                                                                             agtaTGATTTGGTGGCTACATTAGCATGTTTTGCACACCTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAA‑TGCTGTCAAAAGGCGTATTGTCAGCGCGACTTTTCAACCTTATTTATGGCTACCATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTACTTCAAATTTTTCATCGGGTTATTTGCGCTGGTCAACCCGGTAggg                                                                                           >  1:245445/1‑250 (MQ=255)
                                                                                                                                                                      tggtggCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAA‑TGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTACTTCAAATTTTTCATCGGGTTATTTGCGCTGGTCAACCCGGTAGGGATTATTccc                                                                                  <  1:177085/250‑1 (MQ=255)
                                                                                                                                                                                                     aGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCATACCTTTTTTGATTTTCCCGTTTACTTCAAATTTTTCATCGGGTTATTTGCGCTGGTCAACCCGGTAGGGATTATTCCCGTCTTTATCAGCATGACCAGTTATCAGACAg                                                   <  1:248118/251‑1 (MQ=255)
                                                                                                                                                                                                                            aaTTTTAAGTCGTTAGGAGGTGAAAA‑TGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTACTTCAAATTTTTCATCGGGTTATTTGCGCTGGTCAACCCGGTAGGGATTATTCCCGTCTTTATCAGCATGACCAGTTATCAGACAGCGGCAGCGCGAAACAAAACTAAcc                           >  2:392041/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTACTTCAAATTTTTCATCGGGTTATTTGCGCTGGTCAACCCGGTAGGGATTATTCCCGTCTTTATCAGCATGACCAGTTATCAGACAGCGGCAGCGCGAAACAAAACTAACCTTAc                       >  1:265283/1‑250 (MQ=255)
                                                                                                                                                                                                                                    gtagtTAGGAGGTGAAAA‑TGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTACTTCAAATTTTTCATCGGGTTATTTGCGCTGGTCAACCCGGTAGGGATTATTCCCGTCTTTATCAGCATGACCAGTTATCAGACAGCGGCAGCGCGAAACAAAACTAACCTTACAGcc                   >  1:446431/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aggaggTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTACTTCAAATTTTTCATCGGGTTATTTGCGCTGGTCAACCCGGTAGGGATTATTCCCGTCTTTATCAGCATGACCAGTTATCAGACAGCGGCAGCGCGAAACAAAACTAACCTTACAGCCAACct              <  1:407275/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gaggTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTACTTCAAATTTTTCATCGGGTTATTTGCGCTGGTCAACCCGGTAGGGATTATTCCCGTCTTTATCAGCATGACCAGTTATCAGACAGCGGCAGCGCGAAACAAAACTAACCTTACAGCCAACctgt            <  2:24906/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTGTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTGTCCCGTTTACTTCAAATTTTTCATCGGGTTATTTGCGCTGGTCAACCCGGGAGGGATTATTCCCGTCTTTATCAGCATGACCAGTTATCAGACAGCGGCAGCGCGGAACAAAACGAACCTTACAGCCAACCTGTCTGGGGCCAt  <  2:449870/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
AACATTAGAGTTATTATCTCTAATGCGTCACTTCCAGGTGGCGTAAGCAAGATTACTCACTTCTGGGTACTGATTACGTGATCCAAATCAAATTTTTGCAAAGCTGACACCTTTCAGCATCGCTTTTCGCCATTATAGCTAACAGTTAATAAATTGTAGTATGATTTGGTGGCTACATTAGCATGTTTTGCACAACTAGATAACAATAACGAATGATAGCAATTTTAAGTAGTTAGGAGGTGAAAAATGCTGTCAAAAGGCGTATTGTCAGCGCGTCTTTTCAACCTTATTTATGGCTAACATTATCCGGCTTTTGCTTCGGAGCTAACCGTGATTCAGACCTTTTTTGATTTTCCCGTTTACTTCAAATTTTTCATCGGGTTATTTGCGCTGGTCAACCCGGTAGGGATTATTCCCGTCTTTATCAGCATGACCAGTTATCAGACAGCGGCAGCGCGAAACAAAACTAACCTTACAGCCAACCTGTCTGTGGCCAT  >  NC_000913/1298268‑1298764

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