Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,7650CT37.0% 20.8 / 27.8 27intergenic (+134/‑58)kefC/folAglutathione‑regulated potassium‑efflux system protein KefC/type 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (10/7);  new base T (6/4);  total (16/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGT  >  NZ_CP009273/49681‑49851
                                                                                    |                                                                                      
cATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTcgacga                                                                                   <  1:556563/90‑1 (MQ=255)
cATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTcgacga                                                                                   <  1:375254/90‑1 (MQ=255)
           ggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCt                                                                        >  2:124933/1‑90 (MQ=255)
           ggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCt                                                                        <  1:424245/90‑1 (MQ=255)
                     gggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAg                                                              >  2:166306/1‑90 (MQ=255)
                          agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCg                                                         >  2:269713/1‑90 (MQ=255)
                          agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCg                                                         >  2:146481/1‑90 (MQ=255)
                          agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCg                                                         >  2:119121/1‑90 (MQ=255)
                          agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCg                                                         >  1:66675/1‑90 (MQ=255)
                          agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTa                                                                     >  2:247774/1‑78 (MQ=255)
                          agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTa                                                                     <  1:247774/78‑1 (MQ=255)
                          agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCg                                                         >  1:57573/1‑90 (MQ=255)
                          agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCg                                                         <  2:57573/90‑1 (MQ=255)
                           gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAATGGCGa                                                        >  1:222591/1‑90 (MQ=255)
                             gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgacga                                                                                   >  1:365338/1‑61 (MQ=255)
                             gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgacga                                                                                   <  2:365338/61‑1 (MQ=255)
                                              gAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTCTATCGGGaa                                     <  1:122911/90‑1 (MQ=255)
                                              gAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaa                                     <  1:146481/90‑1 (MQ=255)
                                              gAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTa                                                                     >  1:335756/1‑58 (MQ=255)
                                              gAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTa                                                                     <  2:335756/58‑1 (MQ=255)
                                                atatAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTa                                                                     <  1:151831/56‑1 (MQ=255)
                                                atatAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTa                                                                     >  2:151831/1‑56 (MQ=255)
                                                atatAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAt                                   >  2:100141/1‑90 (MQ=255)
                                                   taAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtctc                                <  1:20751/90‑1 (MQ=255)
                                                          ttCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc                         >  1:604987/1‑90 (MQ=255)
                                                                           cGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt        <  1:119121/90‑1 (MQ=255)
                                                                           cGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGtt        <  1:335014/90‑1 (MQ=255)
                                                                                 aGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGt  >  2:489632/1‑90 (MQ=255)
                                                                                 aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcgg            >  2:221622/1‑80 (MQ=255)
                                                                                 aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcgg            <  1:221622/80‑1 (MQ=255)
                                                                                 aGTCGACGACGGTTTACGCTTTACGTATAATGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGt  >  2:609556/1‑90 (MQ=255)
                                                                                    |                                                                                      
CATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGT  >  NZ_CP009273/49681‑49851

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