Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,816,5660TC69.2% 31.2 / 17.3 26L125L (TTA→TTGyqaADedA family protein YqaA
Reads supporting (aligned to +/- strand):  ref base T (2/6);  new base C (5/13);  total (7/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATC  >  NC_000913/2816487‑2816646
                                                                               |                                                                                
ggTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCTCTGAACGGTCGCTGCTGCAACTGCAACATAGCGt                                                                                  >  2:484275/1‑80 (MQ=255)
  tGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGCAAc                                                                               <  2:484234/81‑1 (MQ=255)
  tGACCATATGTGCAAGGCTAAAATTAATGCCACCACAACATTGCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGCAAc                                                                               <  2:484235/81‑1 (MQ=255)
      cATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGCAACGCtt                                                                           >  1:484267/1‑81 (MQ=255)
      cATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGCAACGCtt                                                                           <  2:484236/81‑1 (MQ=255)
          tgtgCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGCAACGCTTTAc                                                                        <  2:484237/80‑1 (MQ=255)
                ggCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGCAACGCTTTACCAAGGca                                                                 >  1:484270/1‑81 (MQ=255)
                      cAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGCAACGCTTTACCAAGGCACaaaaa                                                           <  2:484238/81‑1 (MQ=255)
                               ccaccaCATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACAAAAAAAAGATTAc                                                  <  2:484239/81‑1 (MQ=255)
                               ccaccaCATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGCAACGCTTTACCAAGGCACAAAAAAAAGATTAc                                                  <  2:484240/81‑1 (MQ=255)
                                   cacaTCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGt                                              <  1:484241/81‑1 (MQ=255)
                                   cacaTCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGCAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGt                                              >  1:484273/1‑81 (MQ=255)
                                    acaTCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTc                                             <  1:484242/81‑1 (MQ=255)
                                    acaTCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTc                                             <  2:484243/81‑1 (MQ=255)
                                            cccTGAACGGTCGCTGCTGCAACTGCAACATAGCGCAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACgag                                     <  2:484244/81‑1 (MQ=255)
                                                 aaCGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATgcg                                >  2:484280/1‑81 (MQ=255)
                                                          tgctgcAACTGCAACATAGCGCAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCatccatcc                       >  1:484276/1‑81 (MQ=255)
                                                            ctgcaactgcaaCATAGCGTAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATccc                      <  1:484245/80‑1 (MQ=255)
                                                            ctgcaactgcaaCATAGCGTAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCg                     >  2:484278/1‑81 (MQ=255)
                                                            ctgcaactgcaaCATAGCGCAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATccc                      <  1:484246/80‑1 (MQ=255)
                                                               caactgcaacATAGCGCAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGc                    <  1:484248/79‑1 (MQ=255)
                                                               caactgcaacATAGCGCAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCt                   >  1:484277/1‑80 (MQ=255)
                                                               caactgcaacATAGCGCAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCt                   <  2:484247/80‑1 (MQ=255)
                                                                 actgcaacATAGCGTAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAAc                <  1:484249/81‑1 (MQ=255)
                                                                 actgcaacATAGCGCAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAAc                <  2:484250/81‑1 (MQ=255)
                                                                 actgcaacATAGCGCAACGCTTTACAAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAAc                <  2:484251/81‑1 (MQ=255)
                                                                         aTAGCGCAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAg        <  1:484252/81‑1 (MQ=255)
                                                                               tAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATc  <  1:484253/81‑1 (MQ=255)
                                                                               |                                                                                
GGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACAAAAAAAAGATTACCGGTCCCCACGAGATGCGCATCCATCCCGCTAACAGACACAGTAAATC  >  NC_000913/2816487‑2816646

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