Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,877,8160TC72.7% 21.5 / 13.6 22pseudogene (408/1482 nt)yeaVBCCT family transporter YeaV
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base C (9/7);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCT  >  NZ_CP009273/1877728‑1877893
                                                                                        |                                                                             
tGGCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTc                                                                              >  1:182131/1‑90 (MQ=255)
 ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                             >  2:545846/1‑90 (MQ=255)
 ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                             >  1:16172/1‑90 (MQ=255)
 ggCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCg                                                                             >  1:66528/1‑90 (MQ=255)
           aTCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGc                                                                   >  2:764989/1‑90 (MQ=255)
              aTGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCa                                                                <  2:182131/90‑1 (MQ=255)
               tGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCAt                                                               >  2:320258/1‑90 (MQ=255)
                ggCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATg                                                              >  2:144869/1‑90 (MQ=255)
                    tATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGa                                                          >  1:711599/1‑90 (MQ=255)
                                 tgCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGAtt                                             <  1:522319/90‑1 (MQ=255)
                                  gCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGAttt                                            <  2:711599/90‑1 (MQ=255)
                                               ggTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATc                               >  1:12645/1‑90 (MQ=255)
                                                    gAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCAc                          >  2:660163/1‑90 (MQ=255)
                                                     aGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACt                         <  2:125766/90‑1 (MQ=255)
                                                         tttCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCg                     >  2:107233/1‑90 (MQ=255)
                                                          ttCCGGCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCgg                    >  1:337476/1‑90 (MQ=255)
                                                               gCATTATTGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCAc               <  1:790226/90‑1 (MQ=255)
                                                                 attattGCTGCTATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTg             <  1:833752/90‑1 (MQ=255)
                                                                       gctgctATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCAtt       <  2:330918/90‑1 (MQ=255)
                                                                         tgctATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTc     >  1:705173/1‑90 (MQ=255)
                                                                           ctATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTccc   <  1:707033/90‑1 (MQ=255)
                                                                            tATTACCGGCGTCCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCt  <  1:141193/90‑1 (MQ=255)
                                                                                        |                                                                             
TGGCCTCATTAATCATGGCTTATCACTTTCATGTGCGGAAAAACAAAGGTCTGAGCCTTTCCGGCATTATTGCTGCTATTACCGGCGTTCGCCCGCAAGGCCCATGGGGAAAACTGGTCGATTTGATGTTCCTGATCGCCACTGTCGGCGCACTGACCATTTCCCT  >  NZ_CP009273/1877728‑1877893

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: