Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,549,5660TC44.4% 31.3 / 10.5 27S215G (AGT→GGT) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (15/0);  new base C (0/12);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 5.75e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAT  >  NZ_CP009273/3549485‑3549647
                                                                                 |                                                                                 
tttCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCa                                                                           >  1:289961/1‑90 (MQ=255)
       ttCAAGCGAGACGTTATTCCCCGGCCCCGGGCCGCGGGGCGGGTTATGAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGc                                                                    <  2:379095/90‑1 (MQ=255)
                gaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACg                                                           >  2:617457/1‑90 (MQ=255)
                 aCGTTATTCCCCGCCCCCTGGTCGCGCGGCAGGTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGc                                                          <  1:386017/90‑1 (MQ=255)
                           ccGCCCCCGGGGCGCGCGGGGGGTTATAAATTTTTTGTTCAGGGGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                <  2:467119/90‑1 (MQ=255)
                              gCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGGAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                             >  1:123271/1‑90 (MQ=255)
                              gCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                             >  1:232643/1‑90 (MQ=255)
                              gCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                             >  2:572270/1‑90 (MQ=255)
                              gCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                             >  2:371455/1‑90 (MQ=255)
                              gCCCCTGGTCGCGCGGCAGGT‑‑TGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                           >  1:599607/1‑90 (MQ=255)
                                cccTGGCCGCGCGGGAGGTTATAAATTTTTTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                           <  2:289961/90‑1 (MQ=255)
                                   tGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                        >  2:699498/1‑90 (MQ=255)
                                                    ataaaTTTTCTGTCCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                       <  2:695880/87‑1 (MQ=255)
                                                       aaTTTTCTGTTCATGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                    <  1:47822/90‑1 (MQ=255)
                                                          atTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                 >  1:695880/1‑90 (MQ=255)
                                                           ttCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                >  1:650308/1‑90 (MQ=255)
                                                           ttCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                >  1:542403/1‑90 (MQ=255)
                                                            tCTGTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc               <  2:46091/90‑1 (MQ=255)
                                                              tGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc             >  1:339842/1‑90 (MQ=255)
                                                              ggTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc             <  2:140598/89‑1 (MQ=255)
                                                               gTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa            >  1:50556/1‑90 (MQ=255)
                                                                      ggggggAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac     <  2:123271/88‑1 (MQ=255)
                                                                      cAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac     <  1:404913/90‑1 (MQ=255)
                                                                       aGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc               >  2:550386/1‑79 (MQ=255)
                                                                       aGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc               <  1:550386/79‑1 (MQ=255)
                                                                        ggaaaaacCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACaa   <  2:497132/87‑1 (MQ=255)
                                                                         ggAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACaac  >  1:140598/1‑89 (MQ=255)
                                                                                 |                                                                                 
TTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAT  >  NZ_CP009273/3549485‑3549647

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: