Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273899,4380AC58.8% 7.8 / 16.2 17G162G (GGT→GGGybjPlipoprotein
Reads supporting (aligned to +/- strand):  ref base A (6/1);  new base C (0/10);  total (6/11)
Fisher's exact test for biased strand distribution p-value = 5.66e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCC  >  NZ_CP009273/899355‑899525
                                                                                    |                                                                                       
cATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACccgcc                                                                                   >  1:500057/1‑90 (MQ=255)
          aaaaTGCCGAACATTTTAGGGGTTTTATTGGTTGATAAAGGCTTTAGCGTTTCTCGATGGACTGAGGGAGCGTCCCCGCCGTGGCATGaa                                                                         <  1:370340/89‑1 (MQ=255)
                    aCATTTTAGGGGGTTTTTTGGTTGATAAAGGCGTTAGCGTTTCTCGAGGGACTGAGGGAGCGTCCCCGCCGTGGCATGAACGCTACCGcc                                                               <  1:424379/90‑1 (MQ=255)
                     aaaTTTAGGGGGTTTTTTTGTTGTAACAGGTCGTGGGGGATTCTCAATGGACGGCGGGGCCGTCCCCGCCGTGGCATGAACGCTACCGcc                                                               <  2:40734/89‑1 (MQ=255)
                                 tttATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGcaca                                                  >  2:377701/1‑90 (MQ=255)
                                         ttGATCAGGGGGTTAGCTTTTCTCGAGGGACTGACGGGGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGAt                                          <  2:56184/90‑1 (MQ=255)
                                         ttGATCAAGGCGTTAGCGATTCTCGACGGCCGGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGAt                                          <  1:194679/90‑1 (MQ=255)
                                         ttGATAAGGGCGTAAGCGATTCTCGATGGACTGAGGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGAt                                          <  1:275023/90‑1 (MQ=255)
                                         ttGACCAAGGCGTTAGCTTTTCTCGATGGACTGAGGGGGCTTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGAt                                          <  2:193474/90‑1 (MQ=255)
                                         ttGAACAGGGCGTTAGCGATCCTCGAGGGACTGCGGGAGGGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGAt                                          <  2:470489/90‑1 (MQ=255)
                                          tGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATg                                         >  1:360780/1‑90 (MQ=255)
                                                    ctttAGCTTTTCTCGATGGACGGGCGGGGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAacac                                <  2:333706/88‑1 (MQ=255)
                                                               ctcGATGGACGGAGGGAGCGCCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTg                     <  2:315513/90‑1 (MQ=255)
                                                                             ggAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACtt       >  2:277068/1‑90 (MQ=255)
                                                                                gCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGAt                                          >  2:265409/1‑52 (MQ=255)
                                                                                gCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGAt                                          <  1:265409/52‑1 (MQ=255)
                                                                                  gTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATcc  >  1:384922/1‑90 (MQ=255)
                                                                                    |                                                                                       
CATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCC  >  NZ_CP009273/899355‑899525

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