Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273899,4380AC73.7% ‑2.5 / 10.8 19G162G (GGT→GGGybjPlipoprotein
Reads supporting (aligned to +/- strand):  ref base A (5/0);  new base C (0/14);  total (5/14)
Fisher's exact test for biased strand distribution p-value = 8.60e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTG  >  NZ_CP009273/899356‑899527
                                                                                  |                                                                                         
aTTCCTCACAAAATGCCGGAAATTTTGGGGGTTTTTTTGTTTGATAAGGGCGTTAGCGTTTCTCGATGGACTGAGGGAGCGTCCccgccg                                                                                    <  2:8500/90‑1 (MQ=255)
     tCACAAAATGCCGGAAATTTTAGGGGGTTTATTGGTTGATCAAGGCGTTAGCGTTTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCa                                                                               <  1:410769/90‑1 (MQ=255)
             tGCCGGACATTTTAGGGGTTTTATTGGTTGATCAAGGGGTTAGCGTTTCTCGATGGACGGAGGGGGCGTCCCCGCCGTGGCATGAACGCt                                                                       <  2:606354/90‑1 (MQ=255)
              gCCGGACATTTTCGGGGTTTTATTGGTTGATCAAGGCTTTAGCGATTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTa                                                                      <  1:486651/90‑1 (MQ=255)
                                 ttATTGGTTAATAAGGGCGTTAGCTATTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACAt                                                   <  2:232367/90‑1 (MQ=255)
                                        ttGATAAAGGCGTTAGCGATTCTCGAGGGACTGAGGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGAt                                            <  1:488928/90‑1 (MQ=255)
                                        ttGACCAGGGCGTTAGCGTTTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGAt                                            <  2:159717/90‑1 (MQ=255)
                                         tGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATg                                           >  1:550318/1‑90 (MQ=255)
                                             aaaGGCGTTAGCGTTTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGAccc                                       <  1:269752/89‑1 (MQ=255)
                                                    ttAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTg                                >  2:785426/1‑90 (MQ=255)
                                                    ttAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTg                                >  2:173410/1‑90 (MQ=255)
                                                       gCGTTTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGcc                             <  1:292861/90‑1 (MQ=255)
                                                        cttttCCCGAGGGACGGAGGGGGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCt                            <  2:159511/87‑1 (MQ=255)
                                                          aTTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTc                          >  1:278144/1‑90 (MQ=255)
                                                             ctcGATGGACTGCCGGAGCCCCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTg                       <  1:142874/90‑1 (MQ=255)
                                                               cGATGGACGGACGGAGCGCCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGaa                     <  2:581737/90‑1 (MQ=255)
                                                                        gACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAaca            >  2:406461/1‑90 (MQ=255)
                                                                           ggAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACtt         <  1:24458/90‑1 (MQ=255)
                                                                               cGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATc     <  1:596944/90‑1 (MQ=255)
                                                                                  aCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTg  >  1:470588/1‑90 (MQ=255)
                                                                                  |                                                                                         
ATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTG  >  NZ_CP009273/899356‑899527

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