Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,474,4240TG71.4% 40.1 / 11.8 28K43N (AAA→AAC) rpsL30S ribosomal subunit protein S12
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base G (8/12);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGAT  >  NC_000913/3474181‑3474639
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                        
agctttgggACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCCCCACGTACGGTGTGGTAATGAACACCCGGGAGTTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAGCACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAg                                                                                                                                                                                                                   <  1:265735/244‑1 (MQ=255)
     tAGGACGCTTCACGCCATACTTGGAACGAGCCGGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATCTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGgt                                                                                                                                                                                                               <  1:2647/250‑1 (MQ=255)
                       aCTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTCACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATCTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACAt                                                                                                                                                                                            >  1:265819/1‑251 (MQ=255)
                                            ggTCTTTAACGCCGGAGCAGTCAAGCTCCCCACGATCGGTGTGGTAATGAACCCCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCAt                                                                                                                                                                      <  1:158252/251‑1 (MQ=255)
                                                     cGCCGGAGCAGCCAAGCGCACCCCGTACGGTGTGGTAAGGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAgc                                                                                                                                                             <  1:216614/251‑1 (MQ=255)
                                                        cGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAATGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCa                                                                                                                                                          <  1:23243/251‑1 (MQ=255)
                                                           aGCAGTCAAGCGCACCACGTACGGTGTGGTAATGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCCCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGCAACCGTTAGTCAGTAGAACACGGCCTACTTTACGCAGCGCGGAGTTCGGTTTGTTCGGAGAGGTAGTATATACACGAGTACATACGC‑ACAGTTTTTGCGGGAATGCTTCCAGCGCAGGc                                                                                                                                                       >  2:105698/1‑251 (MQ=255)
                                                                               aCGGTGTGGTAATGAACTCCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACttt                                                                                                                                   <  1:245217/251‑1 (MQ=255)
                                                                                    gtgGTAATGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCg                                                                                                                                >  1:288331/1‑249 (MQ=255)
                                                                                                ccccGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGc                                                                                                                  <  1:44127/250‑1 (MQ=255)
                                                                                                             ttAACACGCCCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGt                                                                                                      <  1:136943/250‑1 (MQ=255)
                                                                                                                          cACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCAtt                                                                                        >  1:263437/1‑251 (MQ=255)
                                                                                                                                  aggatcaCGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATCTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCt                                                                                <  1:126025/251‑1 (MQ=255)
                                                                                                                                    gatcaCGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATCTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGCTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTcc                                                                              >  2:37759/1‑251 (MQ=255)
                                                                                                                                             gTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAg                                                                     <  1:221170/251‑1 (MQ=255)
                                                                                                                                                     gCAGGTTGTGTCCTTCACCACCGATGTATGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTg                                                               <  1:247823/249‑1 (MQ=255)
                                                                                                                                                       aGGTTGTGACCTTCACCACCGATCTATGAAGTCGCTTCGAAGCCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATTCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTc                                                           <  2:83392/251‑1 (MQ=255)
                                                                                                                                                           tgtgACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTa                                                        >  1:137706/1‑250 (MQ=255)
                                                                                                                                                                cctcacactccGAGGTAGGTAGTCACTTCGAAATCGTTAGTCAGACTAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGGATATACACGGGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCATCGCAGGCTCGGTGCTTTTCGCAACTTTGCGAGCACTTGGTTTGCGTTCCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACg                                                  <  2:76522/242‑1 (MQ=255)
                                                                                                                                                                       caccGATGTAGGAAGTCACTACGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTGCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATaaa                                           <  1:227427/251‑1 (MQ=255)
                                                                                                                                                                          cGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACg                                        >  1:279861/1‑251 (MQ=255)
                                                                                                                                                                                        aCTTCGAAACCGTTATTCAGACGAACACGTCATAATTTACGCCGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATTCGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTGAGCTTTTGCTTCGTAAACACGGAATAAAACGTCCTCACACaacat                          <  2:265819/251‑4 (MQ=255)
                                                                                                                                                                                         cTTCAAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTAGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCATCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAACTGGTTAACTGTTGCCATTAAATCGCTCCTGGTTTTAGCTTTTGCTTCGTACACACGTACTAAAACGTCCTCACACAATATg                         >  1:177768/1‑251 (MQ=255)
                                                                                                                                                                                                     taggcAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCGC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGACATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAAt             >  2:199711/5‑251 (MQ=255)
                                                                                                                                                                                                             cGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAAGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAgg                      >  1:121299/1‑234 (MQ=255)
                                                                                                                                                                                                             cGAACACGGCATACTTTACGCAGCGCGGAGTGCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTGTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAgg                      <  2:121299/234‑1 (MQ=255)
                                                                                                                                                                                                              gAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGAGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCg    >  1:169542/1‑251 (MQ=255)
                                                                                                                                                                                                                acacGGCATACATTACGCAGCGCGGAGTTCGGTTTGTTAGGAGTGGTAGTATATACACGATTACATACTCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGAt  <  2:230679/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                        
AGCCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCAC‑GTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGAT  >  NC_000913/3474181‑3474639

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