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

ACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556076‑3556345
                                                                                                                                    |                                                                                                                                          
aCTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCa                                                                                                                                    >  1:1028370/1‑140 (MQ=255)
 cTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAg                                                                                                                                   >  2:1644154/1‑140 (MQ=255)
    aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAgtgt                                                                                                                                >  2:1360835/1‑140 (MQ=255)
    aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAgtgt                                                                                                                                >  1:583453/1‑140 (MQ=255)
              cACGGACGACAAAAAAGCGTTGGGTGATTTTTTCACTTGCGACTTAAGGCGGGACGTTATCCCCCGGCCCTTGGGCGCGCGGCAGGTTATAAATTTTCTGTTCA‑GGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCa                                                                                                                      <  2:1028370/140‑1 (MQ=255)
                      acAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                              >  2:401286/1‑140 (MQ=255)
                                       tattttttCTCTTCCGACTTAAAGCGAGACGGTTTTCCCCGGCCCCTGGTCGCGGGGCAGTTTATAAATATTCTGTTCA‑GGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac                                                                                             <  2:444684/136‑1 (MQ=255)
                                         ttttttCACTTTCCACTAGAAGCGAGCCGTTATTCCCCGGCCCTTGGCCGCGGGGCAGTTTATGAATTTTCTGTTCA‑GGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                           <  2:1590235/138‑1 (MQ=255)
                                                aCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa                                                                                    >  1:287263/1‑140 (MQ=255)
                                                aCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa                                                                                    >  1:317850/1‑140 (MQ=255)
                                                aCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa                                                                                    >  1:581412/1‑140 (MQ=255)
                                                aCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa                                                                                    >  1:454375/1‑140 (MQ=255)
                                                  tttCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                  >  2:74613/1‑140 (MQ=255)
                                                  tttCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                  >  1:399038/1‑140 (MQ=255)
                                                   ttCCACTTCAAGCGAGACGTTATTCCCCGGCCGCTGGTCGCGGGGAAGGTTAAGAATTTTCT‑TTCAGTGCAGAAAAAACCCCGCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa                                                                                 <  2:563516/140‑1 (MQ=255)
                                                    tCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                >  1:180817/1‑140 (MQ=255)
                                                       aCTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                             >  1:972043/1‑140 (MQ=255)
                                                           cAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCCCCa                                                                         <  1:1360835/140‑1 (MQ=255)
                                                             aGCGAGCCGTTATTCCCCGCCCCCTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCA‑TGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                       <  2:399038/140‑1 (MQ=255)
                                                              gCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                      >  1:1186683/1‑140 (MQ=255)
                                                                     gttATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATATTCTGTTAA‑TGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca                                                               <  1:401286/138‑1 (MQ=255)
                                                                     ggTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCA‑TGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca                                                               <  2:227273/140‑1 (MQ=255)
                                                                     acagttccccGCCCCCTGGCCGCGGGGCAGTTTATAAATATTTTGTTCA‑TGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca                                                               <  2:155437/133‑1 (MQ=255)
                                                                         tcccccGGCCCCTGGTCGCGGGGCAGGTTATAAATTTTCTGTTCA‑GGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAt                                                           <  2:1319872/137‑1 (MQ=255)
                                                                            cccGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                        >  2:1284536/1‑140 (MQ=255)
                                                                              cgcccccTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCA‑TGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTcc                                                      <  2:239776/137‑1 (MQ=255)
                                                                                 ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  2:966952/1‑140 (MQ=255)
                                                                                 ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  1:1427759/1‑140 (MQ=255)
                                                                                 ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  2:1311968/1‑140 (MQ=255)
                                                                                  cccTGGCCGCGCGGCAGTTTATAAATTTTTTGTTCA‑GGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgcg                                                  <  2:583453/140‑1 (MQ=255)
                                                                                    cTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgctgt                                                >  1:155437/1‑137 (MQ=255)
                                                                                     tGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                               >  1:1623127/1‑140 (MQ=255)
                                                                                       gTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGc                                             >  2:606913/1‑140 (MQ=255)
                                                                                        tCGGGCGGCAGTTTATGAATTTCTTGTACA‑GGGAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                            <  1:1644154/140‑1 (MQ=255)
                                                                                           cgcgGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                         >  1:813852/1‑140 (MQ=255)
                                                                                             cgGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                       >  1:1205017/1‑140 (MQ=255)
                                                                                                gagtttATAAATATTCTGTTCA‑GGGAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                    <  1:219065/136‑1 (MQ=255)
                                                                                                 ggtttATAAATTTTCTGCTCA‑TGGAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                   <  1:476375/137‑1 (MQ=255)
                                                                                                 aGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                   >  2:1008782/1‑140 (MQ=255)
                                                                                                     tATGAATTTTCTGTTCA‑TGCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtg                               <  1:543150/140‑1 (MQ=255)
                                                                                                      aTGAATTTTCTGTTCA‑GGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtgt                              <  2:1205017/140‑1 (MQ=255)
                                                                                                             ttCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       >  1:913217/1‑140 (MQ=255)
                                                                                                             ttCTGTTCA‑TGCAGGGAAAAACTCCCCACCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       >  2:393618/1‑140 (MQ=255)
                                                                                                              tCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                      >  2:1585509/1‑140 (MQ=255)
                                                                                                                  ttCA‑TGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  <  2:972043/140‑1 (MQ=255)
                                                                                                                   tCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAGGGTGTTAAACGATGCCAc                 >  2:1543317/1‑140 (MQ=255)
                                                                                                                    cA‑TGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg                <  2:581412/140‑1 (MQ=255)
                                                                                                                       tGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg              >  1:789832/1‑140 (MQ=255)
                                                                                                                           gggAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag          >  2:1164771/1‑140 (MQ=255)
                                                                                                                                 aaCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc    >  2:438537/1‑140 (MQ=255)
                                                                                                                                   cTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt  <  1:1008782/140‑1 (MQ=255)
                                                                                                                                   cTCCCCACCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt  <  1:393618/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                          
ACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCA‑TGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556076‑3556345

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