Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,775,8680AC53.3% 22.6 / 39.9 30V73G (GTA→GGA) ypjAadhesin‑like autotransporter YpjA/EhaD
Reads supporting (aligned to +/- strand):  ref base A (10/4);  new base C (0/16);  total (10/20)
Fisher's exact test for biased strand distribution p-value = 3.33e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-05
Rejected as consensus: Frequency below/above cutoff threshold.

TAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCT  >  NZ_CP009273/2775779‑2775954
                                                                                         |                                                                                      
tAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATa                                                                                        >  2:1234689/1‑90 (MQ=255)
            tattaCCCTGCCCCAAAAAATTTTTATGCCTCCCTATATCTTAATACCCCCCAGATTAAAAAAGACTTCCATTAGATCCCCCGCCATGaa                                                                            <  1:1177154/90‑1 (MQ=255)
               taGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACAt                                                                         >  2:443507/1‑90 (MQ=255)
                    tcccccAAAAAATTTTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGTTTTCAATTAGATCCCCCGCCATGAACATTctgc                                                                    <  1:15002/88‑1 (MQ=255)
                     gCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTctgct                                                                   >  2:42245/1‑90 (MQ=255)
                     gCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTctgct                                                                   >  2:1177154/1‑90 (MQ=255)
                          acaAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCa                                                              >  2:777677/1‑90 (MQ=255)
                              aaTTGTTTTGCCTTCCTATATCTTAATACCCACCAAATTAAATAAGACTTCAATTAGACCCCCCGCCATGAACATTCTGCTGGCCATGGt                                                          <  1:39138/90‑1 (MQ=255)
                                  gTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGat                                                      <  2:820228/90‑1 (MQ=255)
                                       gccTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAatat                                                 <  2:864149/88‑1 (MQ=255)
                                        aCTTCCTATATTTTGATACCCCCCAGATTAAATAAGTTTTCAATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATg                                                <  1:1107614/90‑1 (MQ=255)
                                         cTTCCTATATCTTGATATCCACCAGATTAAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGa                                               <  2:714414/90‑1 (MQ=255)
                                                 ctCTTGATACCCCCCAGATTAAATAAGTTTTCAATTATATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTAtt                                       <  1:927276/89‑1 (MQ=255)
                                                 atCTTGATATCCACCAGATTCAATAAGACTTCCAATAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTAtt                                       <  2:177738/90‑1 (MQ=255)
                                                      gATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATg                                                            >  1:275721/1‑64 (MQ=255)
                                                      gATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATg                                                            <  2:275721/64‑1 (MQ=255)
                                                      gATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATttgttg                                  >  1:944689/1‑90 (MQ=255)
                                                      gATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATttgttg                                  >  2:305456/1‑90 (MQ=255)
                                                          tCCACCAAATTAAATAATTCTTCAATTAAATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTAcc                              <  1:682276/90‑1 (MQ=255)
                                                                 gATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCa                       >  2:1024814/1‑90 (MQ=255)
                                                                   ttCAATAAGACTTCCATTAGATACCCCGCCATGAACATTctgc                                                                    <  2:1066118/43‑1 (MQ=255)
                                                                   ttCAATAAGACTTCCATTAGATACCCCGCCATGAACATTctgc                                                                    >  1:1066118/1‑43 (MQ=255)
                                                                    tCAATAAGACTTCTATTATATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATcc                    <  1:777677/90‑1 (MQ=255)
                                                                              ttCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGt          <  2:1087551/90‑1 (MQ=255)
                                                                              ttCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGt          <  1:1024814/90‑1 (MQ=255)
                                                                               tCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGtt         >  2:718479/1‑90 (MQ=255)
                                                                                ccATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTg        <  2:621871/90‑1 (MQ=255)
                                                                                caaTTATATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTg        <  2:889608/88‑1 (MQ=255)
                                                                                  aTTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGcc      <  1:396426/90‑1 (MQ=255)
                                                                                   ttAGACACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCa     <  1:533915/90‑1 (MQ=255)
                                                                                      gATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCt  <  2:534269/90‑1 (MQ=255)
                                                                                         |                                                                                      
TAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCT  >  NZ_CP009273/2775779‑2775954

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