Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,556,2070TC34.9% 58.3 / 19.0 43S215G (AGT→GGT) ‡rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (25/3);  new base C (0/15);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 5.38e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556070‑3556345
                                                                                                                                         |                                                                                                                                          
cTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTcc                                                                                                                                          >  1:206456/1‑140 (MQ=255)
               cTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGa                                                                                                                           >  1:930841/1‑140 (MQ=255)
                  ttCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGAttt                                                                                                                        >  1:625429/1‑140 (MQ=255)
                    cACCGACGACAAAAAAGCGTTCGGTTATTTTTTCCCTTTCGACTTCAAGGGGGCCGGTATCCCCCGCCCCGTGGTCGCGCGGGAGGTTATGAATTTTTTGTTCATGGAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCa                                                                                                                      <  2:818601/140‑1 (MQ=255)
                    cACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCa                                                                                                                      >  1:1156856/1‑140 (MQ=255)
                        gacgacAAAAAGGCGTTGGGTAATTTTTTCTCTTTCGACTTCAAGCGAGACGTTTTTCCCCGGCCCCGGGTCGCGGGGCAGTTTATGAATATTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                                  <  2:510746/140‑1 (MQ=255)
                          cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGTGATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                                >  1:1138188/1‑140 (MQ=255)
                           aaaaaaaagGCGTTGGGTGTTTTTTTCATTTTCGACTTAAAGCGAGACGTTTTCCCCCGCCCCCGGGTCGCGCGGCAGGTTATAATTTTTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAg                                                                                                               <  2:625429/137‑1 (MQ=255)
                                        tCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                  >  2:520605/1‑140 (MQ=255)
                                             tattttttCCCTTTCGACTTCAAGCGAGACGTTATTCCCCGGCCCCGGGCCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTTACGCTCAGCAACATGGCGCGGCacac                                                                                             <  2:670174/136‑1 (MQ=255)
                                             tattttttCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGAAGTTTATAAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac                                                                                             <  2:930841/136‑1 (MQ=255)
                                                          tCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                <  2:1105577/140‑1 (MQ=255)
                                                          tCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGACACACCAGCTAATAGAAc                                                                                >  1:1312647/1‑140 (MQ=255)
                                                           cGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACt                                                                               >  1:89823/1‑140 (MQ=255)
                                                            gACTTCAAGCGAGACGTTATTCCCCGCCCCCTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                              <  2:63360/140‑1 (MQ=255)
                                                                 aaaGCGAGGCGTTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGACACACCAGCTAATAGAACTTCTCCa                                                                         <  2:1312647/139‑1 (MQ=255)
                                                                      gagggggTTTTCCCCGCCCCCGGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                    <  2:1156856/135‑1 (MQ=255)
                                                                      gagaGGGTATTCCCCGGCCCCGGGTCGCGCGGCAGTTTATAAATATTCTGTTCAGGAAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                    <  2:696587/140‑1 (MQ=255)
                                                                      gagaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                    <  2:653485/140‑1 (MQ=255)
                                                                      agaaggttATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATATTCTTTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                    <  2:755037/133‑1 (MQ=255)
                                                                           ggTTTTCCCCGCCCCCGGGTCGCGCGGCAGGTTATAAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGTGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca                                                               <  2:1364156/140‑1 (MQ=255)
                                                                           ggTATTCCCCGCCCCCGGGTCGCGCGGCAGGTTATAAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca                                                               <  2:844276/140‑1 (MQ=255)
                                                                                 aaccGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                         >  1:510208/3‑140 (MQ=255)
                                                                                  cccGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGTGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                        >  1:1364156/1‑140 (MQ=255)
                                                                                  cccGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                        >  2:864411/1‑140 (MQ=255)
                                                                                       ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtctgtct                                                   >  1:755037/1‑134 (MQ=255)
                                                                                          tggggcgcgtgGCAGGTTAGTATTTCTCTGTTCGGGCGGGGAAATCCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCACGGCACACCAGCTAATAGAACTTCTCCACGCGTCTGCACAATGTTCCCGCGCt                                                <  2:982104/135‑1 (MQ=255)
                                                                                                 cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                         >  1:1209295/1‑140 (MQ=255)
                                                                                                     gCAGGTTATAAATTTTCTGTTCATGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGc                                     <  1:1433899/140‑1 (MQ=255)
                                                                                                      gagtttAAGAATTTTCTGTTCATGGAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                    <  1:792849/136‑1 (MQ=255)
                                                                                                      cAGGTTATGAATATTCTGTTCATGCAGGGAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                    <  1:864411/140‑1 (MQ=255)
                                                                                                       agtttATAAATTTTCTGTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                   <  2:89823/137‑1 (MQ=255)
                                                                                                        gtttATAAATATTCTGTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAg                                  <  2:274801/138‑1 (MQ=255)
                                                                                                                attttCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaaa                          <  2:206456/137‑1 (MQ=255)
                                                                                                                atatTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaaa                          >  2:1176660/1‑140 (MQ=255)
                                                                                                                    tCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                      >  2:1349608/1‑140 (MQ=255)
                                                                                                                    tCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                      >  1:1210367/1‑140 (MQ=255)
                                                                                                                        ttCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  2:631155/1‑140 (MQ=255)
                                                                                                                        ttCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  2:1103743/1‑140 (MQ=255)
                                                                                                                         tCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCAc                 >  2:492508/1‑140 (MQ=255)
                                                                                                                           aTGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACTATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg               >  2:568864/1‑140 (MQ=255)
                                                                                                                            tGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg              >  1:754627/1‑140 (MQ=255)
                                                                                                                             gCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc             >  2:32967/1‑140 (MQ=255)
                                                                                                                             gCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc             >  1:1399700/1‑140 (MQ=255)
                                                                                                                             gCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc             >  2:678850/1‑140 (MQ=255)
                                                                                                                                 ggAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgaga         <  1:678850/140‑1 (MQ=255)
                                                                                                                                  gAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGTTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc        >  1:1376075/1‑140 (MQ=255)
                                                                                                                                        cTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt  <  2:510208/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
CTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556070‑3556345

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑