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

GTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAT  >  NZ_CP009273/3549472‑3549647
                                                                                       |                                                                                        
gTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGaa                                                                                        >  1:316927/1‑90 (MQ=255)
             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:371455/1‑90 (MQ=255)
                                           gCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                             >  2:572270/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)
                                                                                      ggAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACaac  >  1:140598/1‑89 (MQ=255)
                                                                                       |                                                                                        
GTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAT  >  NZ_CP009273/3549472‑3549647

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