Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,549,5590GA44.7% 27.8 / 18.2 38S217F (TCC→TTC) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base G (14/7);  new base A (0/17);  total (14/24)
Fisher's exact test for biased strand distribution p-value = 1.73e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGA  >  NZ_CP009273/3549439‑3549697
                                                                                                                        |                                                                                                                                          
aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAgtg                                                                                                                          >  1:888243/1‑139 (MQ=255)
aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAgtg                                                                                                                          >  1:107035/1‑139 (MQ=255)
 cgcgCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAgtg                                                                                                                          >  1:43816/1‑138 (MQ=255)
          cACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGAttt                                                                                                                 >  2:383253/1‑138 (MQ=255)
          cACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGAttt                                                                                                                 >  2:487302/1‑138 (MQ=255)
             cgacgaAAAAAAAGCGTCCGGTGATTTTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCGGTTTATGAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACg                                                                                                             <  1:175910/139‑1 (MQ=255)
                             ttgggTGATTTTTCCCCTTCCGACTACAGGCGGGACGTTATTCCCGGCCCCCGGGTCGCGCGGCAGTTTATAAATTTTCTGTTCTTGAAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                             <  1:63975/136‑1 (MQ=255)
                             ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                             >  1:486102/1‑139 (MQ=255)
                                 tttttttttcCACTTTCGACTCCAAGCGAGCCGTTTTTCCCCGGCCCTTGGTCGCGCGGCAGTTTATGAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                         <  2:384000/133‑1 (MQ=255)
                                 ttGATATTTTCCCTTTCGACTCCAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                         <  2:888243/138‑1 (MQ=255)
                                    atatTTTCCCTTCCGACTTAAGGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGAAGGTTATGAATATTCTGTCCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca                                                                                       <  1:487302/138‑1 (MQ=255)
                                    atatTTCCACTTTCGACTTAAAGCGAGACGGTATTCCCCGGCCCCTGGCCGCGCGGCAGGTTATAAATATTCTGTTCATGCAGGGAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca                                                                                       <  1:19676/138‑1 (MQ=255)
                                                 cGACTTCAAGCGAGCCGTTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATGAATATTCTGTCCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa                                                                          <  2:107035/138‑1 (MQ=255)
                                                  gACTTAAAGCGAGAGGGTATTCCCCGCCCCGTGGTCGCGCGGCGGTTTATGAATTTTCTGTTCATGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                         <  1:700636/138‑1 (MQ=255)
                                                   ccTTAAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCAGGTTATGAATTTTCTGTTCATGCAGGAAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                       <  1:400538/138‑1 (MQ=255)
                                                        aggCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                  <  2:43816/137‑1 (MQ=255)
                                                        aaGCGAGGCGGTTTTCCCCGGCCCCTGCTCGCGCGGCAGGTTATAAATATTCTGTTCAGGCAGGGAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                  <  2:638518/139‑1 (MQ=255)
                                                        aaGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCAGGTTATAAATATTCTGTTCAGGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                  <  1:450337/139‑1 (MQ=255)
                                                                  gTATTCCCCGGCCCCTGTGCGCGCGGCAGGTTATAAATATTCTTTTCATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca                                                        <  2:434971/139‑1 (MQ=255)
                                                                        cccGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                  >  1:814693/1‑139 (MQ=255)
                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                             >  1:494016/1‑139 (MQ=255)
                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                             >  2:444645/1‑139 (MQ=255)
                                                                             ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTACCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                             >  2:543732/1‑139 (MQ=255)
                                                                                  ggTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                        >  1:273733/1‑139 (MQ=255)
                                                                                     cgcgcgGCAGGTTATAAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                     <  2:791836/139‑1 (MQ=255)
                                                                                       cgcgGCAGGTTATGAATATTCTGTTCAAGGAGGGAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCaa                                   <  1:383253/139‑1 (MQ=255)
                                                                                            cAGGTTATAAATATTCTGTTCATGCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGc                              <  1:220115/139‑1 (MQ=255)
                                                                                             agtttAGGAAAATTCTGTAAAGGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCGCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                             <  1:644545/136‑1 (MQ=255)
                                                                                             aGGTTATGAATATTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                             <  1:510256/139‑1 (MQ=255)
                                                                                             aGGTTATAAATATTCTGTTCATGCGGGAAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                             <  1:572817/139‑1 (MQ=255)
                                                                                              gtttATGAATATTCTGTTCATGCAGGGAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                             <  1:658602/136‑1 (MQ=255)
                                                                                                     aaTATTCTGTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                     <  1:147746/139‑1 (MQ=255)
                                                                                                     aaTATCTTGGTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                     <  2:722445/139‑1 (MQ=255)
                                                                                                          tCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                             <  2:391320/94‑1 (MQ=255)
                                                                                                          tCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                             >  1:391320/1‑94 (MQ=255)
                                                                                                          tCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACACTGTTCGCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                >  2:690782/1‑139 (MQ=255)
                                                                                                                cATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa           <  1:641893/138‑1 (MQ=255)
                                                                                                                   gCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg       >  2:628456/1‑139 (MQ=255)
                                                                                                                      gggAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                            >  2:196206/1‑67 (MQ=255)
                                                                                                                      gggAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                            <  1:196206/67‑1 (MQ=255)
                                                                                                                        gAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                  >  2:447786/1‑75 (MQ=255)
                                                                                                                        gAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                  <  1:447786/75‑1 (MQ=255)
                                                                                                                        gAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgaga  >  1:121934/1‑139 (MQ=255)
                                                                                                                        |                                                                                                                                          
ACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGA  >  NZ_CP009273/3549439‑3549697

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