Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,134,4070AC57.1% 7.8 / 14.0 14G229G (GGT→GGGasmAouter membrane assembly protein AsmA
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base C (0/8);  total (4/10)
Fisher's exact test for biased strand distribution p-value = 1.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGAT  >  NZ_CP009273/2134319‑2134491
                                                                                        |                                                                                    
tttGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCAc                                                                                     <  2:93653/90‑1 (MQ=255)
    gTTAATTGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTg                                                                                 >  1:424301/1‑90 (MQ=255)
                  gTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTaa                                                                   >  1:536446/1‑90 (MQ=255)
                  gTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTaa                                                                   >  2:123038/1‑90 (MQ=255)
                                tcccACTGGGCCTGAAACCCCCCCTGCCCCTTAATCCCTTTTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAgcc                                                     <  2:242020/88‑1 (MQ=255)
                                                  cTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGa                                   >  1:22303/1‑90 (MQ=255)
                                                         gcccctgAATACCTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAg                            <  1:123038/90‑1 (MQ=255)
                                                                       tGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTTTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgta              <  1:12035/90‑1 (MQ=255)
                                                                        tttttGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgtac             <  1:135115/89‑1 (MQ=255)
                                                                        tttttGGAAAACGGGCCCCCTGCAACTGCCAATTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgtac             <  2:499738/89‑1 (MQ=255)
                                                                         ttttGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTAAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgtacc            <  1:485919/90‑1 (MQ=255)
                                                                         ttttGGTAAATCGCCCCCCTCCACCCGCAATTTATTTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgtacc            <  2:611100/90‑1 (MQ=255)
                                                                                 aaTCGGCCCCCTGCAACTGCCATTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAggg    <  1:302904/90‑1 (MQ=255)
                                                                                   tCGGCACCCTGCAACTGCCAGTTATTTTTTTAAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGat  <  2:536446/90‑1 (MQ=255)
                                                                                        |                                                                                    
TTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGAT  >  NZ_CP009273/2134319‑2134491

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