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

TCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556071‑3556345
                                                                                                                                         |                                                                                                                                          
tCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTccc                                                                                                                                         >  2:20829/1‑140 (MQ=255)
ccggccttGACGCGCTTTTCACCGACGAAAAAAAAGCGTTCGGTGATTTTTTCACTTTCGCCTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCA‑GGAAAAAcccccc                                                                                                                                        <  2:514597/133‑1 (MQ=255)
   gCACTGCCGCGTTTTTCACCGACGACAAAAAAGCGTTGGGTGATTTTTTCCCTTTCGACTTCAAGGGAGCCGGTATTCCCCGCCCCCGGGCCGCGGGGCAGTTTATAAATATTCTGTT‑CATGAGGGAAAAAACCCCCCGc                                                                                                                                      <  2:114835/140‑1 (MQ=255)
         aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAgtgt                                                                                                                                >  1:1263190/1‑140 (MQ=255)
         aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAgtgt                                                                                                                                >  1:646704/1‑140 (MQ=255)
              cTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCATGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGa                                                                                                                           >  2:120979/1‑140 (MQ=255)
                   cACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCa                                                                                                                      >  1:1730289/1‑140 (MQ=255)
                       gacgacAAAAAGGCGTTCGTTGATTTTTTCTCTTCCGACTTCAAGGGAGACGTTATCCCCCGCCCCCTGGTCGCGGGGCAGTTTATAAATATTTTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                                  <  2:1314092/140‑1 (MQ=255)
                            aaaaaaaGCGTTCGGTGATTTTTTCTCTTTCGACTTAAGGGGAGAGGGTATTCCCCGCCCCTTGGTCGCGCGGCGGTTTATGAATTTTTTTTT‑CAGGCAGGAAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCa                                                                                                             <  2:80641/139‑1 (MQ=255)
                                      ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                   >  1:863855/1‑140 (MQ=255)
                                      ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                   >  1:1515443/1‑140 (MQ=255)
                                       tCGTTGTTATTTTCTTTTTCGACTTCAAGAGAGCCGGTGTCACCCGGCCCCCCGGCGCGGGGAAGTTTAGTAATTTTTTGCT‑GATGCGGGGAAAATCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                  <  2:764534/140‑1 (MQ=255)
                                         ggTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTT‑CATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                                <  1:531798/140‑1 (MQ=255)
                                            tattttttCTTTTTCGATTTCAAGGGAGCCGTTTTTCCCCGCCCCCGGGTCGCGCGGCAGGTTATAAATTTTCTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac                                                                                             <  2:1587717/136‑1 (MQ=255)
                                              ttttttCTCTTTCGACTAAAAGCG‑GGCGGTTTTTCCCCGCCCCCGGGCCGGGGGGCGGTTTTAAAATTTTTGTTAAATGGGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                           <  2:423306/138‑1 (MQ=255)
                                              tatTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                           >  1:1693716/1‑140 (MQ=255)
                                                      ttttCGACTTCAAGCGAGACGGTATTCCCCGCCCCCTGGCCGCGCGGCAGTTTATGAATATTCTGTT‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                                   <  2:528881/139‑1 (MQ=255)
                                                        ttCGACTTCAAGGGAGACGGTATTCCCCGCCCCCTGGTCGCGGGGCAGTTTATAAATTTTCTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa                                                                                 <  2:1468693/140‑1 (MQ=255)
                                                        ttCGACTTAAAGCGGGACGTTTTCCCCGGGCCCCTGGTCGCGGGGCAGTTTATAAATTTTCTGTT‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa                                                                                 <  2:1391384/140‑1 (MQ=255)
                                                            aCTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                             >  1:1065179/1‑140 (MQ=255)
                                                                cAAGCGAGCGGTTATCCCCCGCCCCCGGGCCGCGCGGCAGTTAATAATTTTTCTGTC‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                         <  2:863855/140‑1 (MQ=255)
                                                                aaaGCGAGACGGTATTCCCCGGCCCCTGGTCGCGGGGCTGTTTATGAATTTTTTTTT‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                         <  2:1543575/139‑1 (MQ=255)
                                                                  agggagaCGGTATCCCCCGGCCCCGGGTCGCGCGGAAGTTTATAATTTTTTTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                       <  2:381976/137‑1 (MQ=255)
                                                                     gagcgggTATCCCCCGGCCCCGGGCCGGGGGGAGGTTTATAAATTTTCTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                    <  2:1730289/135‑1 (MQ=255)
                                                                       gaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                  >  2:1307005/1‑140 (MQ=255)
                                                                          ggTTTTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca                                                               <  1:1594844/140‑1 (MQ=255)
                                                                               tACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATg                                                          >  1:1423621/1‑140 (MQ=255)
                                                                                      ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGGGATTTCACGCTCAGCAACATGGCTCGGCCCCCGACCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  2:1083691/1‑140 (MQ=255)
                                                                                      ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  2:213816/1‑140 (MQ=255)
                                                                                      ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  2:1412072/1‑140 (MQ=255)
                                                                                      ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCCCCACGCATCTGCACAATGTTCCcgc                                                   >  1:1715265/1‑140 (MQ=255)
                                                                                      ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGAGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  1:1360271/1‑140 (MQ=255)
                                                                                      ccccGGGCCGCGGGGCAGGTTATAAATTTTCTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   <  2:793549/140‑1 (MQ=255)
                                                                                           ggTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg                                              >  1:149148/1‑140 (MQ=255)
                                                                                                cgcgGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                         >  1:1080043/1‑140 (MQ=255)
                                                                                                    gCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGc                                     >  1:614657/1‑140 (MQ=255)
                                                                                                     gaGGTTATGAATTTTCTGTT‑CAGGGGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                    <  2:973954/139‑1 (MQ=255)
                                                                                                     aagtttATAAATTTTCTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                    <  1:1695622/136‑1 (MQ=255)
                                                                                                     aaGGTTATAAATATTCTTTT‑CATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                    <  1:1144456/139‑1 (MQ=255)
                                                                                                      aGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                   >  1:580202/1‑140 (MQ=255)
                                                                                                                  ttCTGTT‑CATGGAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       <  1:1748869/140‑1 (MQ=255)
                                                                                                                   tCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACACTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                      >  2:1696652/1‑140 (MQ=255)
                                                                                                                    ttGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg                     <  1:1402342/139‑1 (MQ=255)
                                                                                                                      ttt‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                   <  2:1423621/139‑1 (MQ=255)
                                                                                                                       tt‑CATGGAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  <  2:1515443/140‑1 (MQ=255)
                                                                                                                       tt‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  2:1187468/1‑140 (MQ=255)
                                                                                                                            tGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg              >  2:78232/1‑140 (MQ=255)
                                                                                                                            ggaaGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg              <  2:549597/137‑1 (MQ=255)
                                                                                                                            ggCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg              <  2:1080043/139‑1 (MQ=255)
                                                                                                                             gCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCGACGGGc             >  1:863508/1‑140 (MQ=255)
                                                                                                                                  gAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc        >  2:574192/1‑140 (MQ=255)
                                                                                                                                       aCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCg   >  2:785628/1‑140 (MQ=255)
                                                                                                                                        cTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCCCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt  <  2:1715265/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
TCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556071‑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: ‑