Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,134,4070AC55.0% 17.8 / 27.9 20G229G (GGT→GGGasmAouter membrane assembly protein AsmA
Reads supporting (aligned to +/- strand):  ref base A (5/4);  new base C (0/11);  total (5/15)
Fisher's exact test for biased strand distribution p-value = 8.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-05
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCG  >  NZ_CP009273/2134319‑2134495
                                                                                        |                                                                                        
tttGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCAc                                                                                         <  2:336974/90‑1 (MQ=255)
    gTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTg                                                                                     >  2:75899/1‑90 (MQ=255)
               ggCGTTTATGTGCCTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGt                                                                          <  1:144830/90‑1 (MQ=255)
                                tcccACTGGGCCTGAAACCCCCCCTGCCCCTGAATCCCTTTTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAgcc                                                         <  1:336954/88‑1 (MQ=255)
                                tcccACTGGGCCGGGAACCCCCCCCCCCCCTAAATCCCTTTTTTTGGAAAATCGGCCCCCTCCAACTGCCAGTTAATTTTTTCAATAgcc                                                         <  1:401325/88‑1 (MQ=255)
                                     cTGGGCCGGAAACCCCCCCTCCCCCTGAATCCTTTTTTTTGGAAAATCGGCCCCCTGCACCTGCCAGTTAATTTGTTCAATAGCCGCCGt                                                    <  1:546641/90‑1 (MQ=255)
                                                 aCTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAgcc                                                         <  2:752117/73‑1 (MQ=255)
                                                 aCTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAgcc                                                         >  1:752117/1‑73 (MQ=255)
                                                        tgcccctgAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGTGGATAATCa                                 >  2:218605/1‑90 (MQ=255)
                                                              tgAATCCCTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCa                           <  2:324315/90‑1 (MQ=255)
                                                                   cctttGTTTTGGAAAATCGGCCCCCTGCAACTGCCATTTATTTTGTTAAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtac                      <  1:86321/87‑1 (MQ=255)
                                                                    ccTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGt                                                                          >  1:592445/1‑37 (MQ=255)
                                                                    ccTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGt                                                                          <  2:592445/37‑1 (MQ=255)
                                                                       ttttttGGAAAACCGCCCCCCTGCAACTGCCAGTTATTTTGTTAAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgta                  <  1:275520/88‑1 (MQ=255)
                                                                         ttttGGAAAATCGCCCCCCTGCAACTGCCAGTTAATTTTTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgtacc                <  1:72709/90‑1 (MQ=255)
                                                                             ggTAAATCGGCCCCCTGAAACTGCCAGTTAATTTGTTAAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTa            <  2:493641/90‑1 (MQ=255)
                                                                             ggTAAACCGGCCCCCTGCAACTGCCAGTTAATTTTTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTa            <  2:667606/90‑1 (MQ=255)
                                                                               aaaaTCGCCCCCCTGCAACTGCAAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAg          <  1:223847/89‑1 (MQ=255)
                                                                                   ccGGCCCCCTCCAACTGCCAGTTATTTTTTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGat      <  1:179148/89‑1 (MQ=255)
                                                                                       cACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCg  >  1:892192/1‑90 (MQ=255)
                                                                                        |                                                                                        
TTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCG  >  NZ_CP009273/2134319‑2134495

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