Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,458,3860AG73.7% 21.3 / 10.4 19E369G (GAA→GGA) gspLtype II secretion system protein GspL
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base G (6/8);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCC  >  NZ_CP009273/3458307‑3458470
                                                                               |                                                                                    
cTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGaa                                                                            >  1:422879/1‑90 (MQ=255)
            aTGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTc                                                                <  2:199723/90‑1 (MQ=255)
              gAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAg                                                              <  2:337201/90‑1 (MQ=255)
                        agTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATgg                                                    <  2:416623/90‑1 (MQ=255)
                          tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                  <  2:464329/90‑1 (MQ=255)
                          tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                  <  2:271106/90‑1 (MQ=255)
                          tGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGAc                                                  >  2:162722/1‑90 (MQ=255)
                                       gatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAg                                                >  2:62234/1‑79 (MQ=255)
                                       gatcgTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAg                                                <  1:62234/79‑1 (MQ=255)
                                             ttCTGTGAGTTAACACAGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGaaa                               <  2:458992/90‑1 (MQ=255)
                                             ttCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGaaa                               <  1:90191/90‑1 (MQ=255)
                                                             aGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTc                                                                >  2:34812/1‑41 (MQ=255)
                                                             aGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTc                                                                <  1:34812/41‑1 (MQ=255)
                                                             aGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCAtggtg               >  2:334861/1‑90 (MQ=255)
                                                             aGTCCTGGCTTCCGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCAtggtg               >  2:423897/1‑90 (MQ=255)
                                                                        ccGATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAt    >  1:45769/1‑90 (MQ=255)
                                                                        ccGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAt    >  1:385798/1‑90 (MQ=255)
                                                                        ccGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAAt    >  1:100585/1‑90 (MQ=255)
                                                                          gATGGGAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATcc  <  1:423897/90‑1 (MQ=255)
                                                                               |                                                                                    
CTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCC  >  NZ_CP009273/3458307‑3458470

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