Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273737,5620AC55.6% ‑5.6 / 17.9 18S151A (TCT→GCT) dtpDdipeptide permease DtpD
Reads supporting (aligned to +/- strand):  ref base A (7/1);  new base C (0/10);  total (7/11)
Fisher's exact test for biased strand distribution p-value = 2.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGCCAACCGCCGCCAGGCCAAAGCCCATCGCCCAGCTGTACTCTTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGAC  >  NZ_CP009273/737481‑737645
                                                                                 |                                                                                   
tGCCAACCGCCGCCAGGCCAAAGCCCATCGCCCAGCTGTACTCTTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGt                                                                             >  2:306919/1‑90 (MQ=255)
       cgccgcCAGGCCAAAGCCCACCGCCCAGCTGTCCTCTTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGCCCCCACGTTACccgc                                                                      <  2:325568/90‑1 (MQ=255)
            ccAGGCCAAAGCCCATCGCCCGGCTGTCCTCTTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGCCCCCACGTTACCCGCCGCAt                                                                 <  2:590402/90‑1 (MQ=255)
               ggCCAAAGCCCATCGCCCAGCTGTACTCTTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACa                                                              >  1:47473/1‑90 (MQ=255)
                         caccGCCCGGCTGTCCTCTCCCGGGGCGAACCCACAGGCGATGGGGGCAAAAAAAGCCCCCCCGTTCCCCGCCGCATACATCAGCGAGaa                                                    <  2:661841/87‑1 (MQ=255)
                               ccAGCTGTCCTCTTCCGGGGCGTAACCACAGGCGATAGGTGCGATAATAGCCCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGcc                                              <  1:306919/90‑1 (MQ=255)
                                   cTGTCCTTTCCCGGGGCCAAACCCCAGGCGAAAGGGGCGATAAAAGCCCCCCCGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCa                                          <  2:312151/90‑1 (MQ=255)
                                         tctTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGCCCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGAcgc                                    <  1:36364/90‑1 (MQ=255)
                                          ctTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATc                                                            <  1:87980/65‑1 (MQ=255)
                                          ctTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATc                                                            >  2:87980/1‑65 (MQ=255)
                                          ctTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGAcgcg                                   >  1:629288/1‑90 (MQ=255)
                                               tGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCg                              >  2:29138/1‑90 (MQ=255)
                                                 ggCGTAACCACAGGCGATAGGTGCGATAATAGCCCCCCCGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGt                            <  1:149797/90‑1 (MQ=255)
                                                 ggCGAACCCACAGGCGATGGGGGCAAAAAAACCCCCCCCTTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGt                            <  2:629288/90‑1 (MQ=19)
                                                          acaGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATa                   >  1:748609/1‑90 (MQ=255)
                                                                       gCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGc      >  1:676159/1‑90 (MQ=255)
                                                                           aaataGCCCCCCCTTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGAc  <  2:676159/89‑1 (MQ=38)
                                                                           aaataGCCCCCATGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGAc  <  2:416247/89‑1 (MQ=38)
                                                                                 |                                                                                   
TGCCAACCGCCGCCAGGCCAAAGCCCATCGCCCAGCTGTACTCTTCCTGGGCGTAACCACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGAC  >  NZ_CP009273/737481‑737645

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