Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,775,8680AC52.9% 10.6 / 18.9 17V73G (GTA→GGA) ypjAadhesin‑like autotransporter YpjA/EhaD
Reads supporting (aligned to +/- strand):  ref base A (7/1);  new base C (0/9);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 4.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-04
Rejected as consensus: Frequency below/above cutoff threshold.

TTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCT  >  NZ_CP009273/2775785‑2775954
                                                                                   |                                                                                      
ttGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGc                                                                                  >  2:418333/1‑90 (MQ=255)
              tcccccACAAAATTGTTATGACTTCCTATATCTTGATACCCCCCAGATTCAATAAGACTCCCATTAGATCCCCCGCCATGAACATTctgc                                                                    <  1:418333/88‑1 (MQ=255)
               gCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTctgc                                                                    <  2:171776/89‑1 (MQ=255)
               gCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTctgct                                                                   >  1:261298/1‑90 (MQ=255)
               gCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTctgc                                                                    >  1:171776/1‑89 (MQ=255)
                   cacaAAATTTTTATGACTTCCTATATCTTGATATCCCCCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGcc                                                               <  1:259066/90‑1 (MQ=255)
                      aaaaTTTTTATCTCTTCCTATCTCTTTATATCCAACAAATTCAAAAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATg                                                            <  1:491517/90‑1 (MQ=255)
                        aaTTTTTATGCCTCCCTATATCTTAATACCCACCAGATTAAAAAAGACTTCCATTAAATCCCCCGCCATGAACATTCTGCTGGCCATGGt                                                          <  1:310541/90‑1 (MQ=255)
                            gTTATGACTTCCTATATCTTGATATCCCCCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGat                                                      <  2:261298/90‑1 (MQ=255)
                                  aCTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATg                                                <  2:41523/90‑1 (MQ=255)
                                         atatCTTGATATCCCCCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTa                                         <  2:429939/90‑1 (MQ=255)
                                               tGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATttgtt                                   >  2:279383/1‑90 (MQ=255)
                                                gATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATttgttg                                  <  2:183606/90‑1 (MQ=255)
                                                               cAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCa                   >  2:279725/1‑90 (MQ=255)
                                                                         tCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGtt         >  2:179081/1‑90 (MQ=255)
                                                                               aGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATc   >  2:497150/1‑90 (MQ=255)
                                                                                gATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCt  <  1:279383/90‑1 (MQ=255)
                                                                                   |                                                                                      
TTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCT  >  NZ_CP009273/2775785‑2775954

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