Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,357,7550AT52.2% 8.7 / 27.3 23P530P (CCA→CCTyhcDfimbrial biogenesis outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base A (6/5);  new base T (0/12);  total (6/17)
Fisher's exact test for biased strand distribution p-value = 4.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTC  >  NZ_CP009273/3357681‑3357836
                                                                          |                                                                                 
tGGCACAATAAACTTAGTGTTTAAAAAAAAAAGGAAGAGGATGGTGATCTGTTTTACGTCAACATCAGTTTTCTTTTTGGCAGAAAAAac                                                                    <  1:662356/90‑1 (MQ=255)
  ccagaaTCACCTTAGTGCTTCAAAAAATAAGGAAGACGATGGTGATTTGTTCTCCGTAAACATCAGTATCCCATTTGGAAGATAAAacac                                                                  <  1:815135/89‑1 (MQ=255)
     gcaTCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTTTTCTACGTCAACATCAGTATTCCTTTTGGCAGATAAAACACCGc                                                               <  1:802406/90‑1 (MQ=255)
     ccaTCAACTTAGTCCTTCAAAAAAAAATGAAGACGATGGTGATCTGTTCTACGTCAACATAAGTATTCCTTTTGGCAGATAAAACACCGc                                                               <  1:390197/89‑1 (MQ=255)
       aTCAACTTAGTGCTTCAGAAAATAATGAAGACGAGGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCca                                                             <  2:312505/90‑1 (MQ=255)
       aTCAACTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCca                                                             <  2:819361/90‑1 (MQ=255)
       aTCAACTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCca                                                             <  1:276233/90‑1 (MQ=255)
        tCAACTTAGTGCTTCAAAAAAAAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCcac                                                            <  1:28544/90‑1 (MQ=255)
          aaCTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCacac                                                          <  2:142708/90‑1 (MQ=255)
          aaCTTAGTGCTTCAAAAAAAAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCCacac                                                          <  1:397054/90‑1 (MQ=255)
           aCTTAGTGCTTAAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATAAAACACCGCCACACt                                                         <  1:78103/90‑1 (MQ=255)
                     tCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGc                                               >  2:488000/1‑90 (MQ=255)
                                             gATCTGTTCTCCGTCAACATCAGTTTTCCTTTTGGAAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTAt                       <  1:489967/90‑1 (MQ=255)
                                                  gTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCca                                                             >  1:25205/1‑47 (MQ=255)
                                                  gTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCca                                                             <  2:25205/47‑1 (MQ=255)
                                                     tttcGTCAAAAACATTTTTCCTTTTGGCAAAAAAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTgg               <  2:846237/87‑1 (MQ=255)
                                                     cTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTgg               >  2:920418/1‑90 (MQ=255)
                                                       aCGTCAACATCAGTATTCCATTTGGCAGATCAAac                                                                    >  1:843465/1‑35 (MQ=255)
                                                       aCGTCAACATCAGTATTCCATTTGGCAGATCAAac                                                                    <  2:843465/35‑1 (MQ=255)
                                                       aCGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGgtg             >  2:274446/1‑90 (MQ=255)
                                                            aaCATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACg                                        <  2:15727/58‑1 (MQ=255)
                                                            aaCATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACg                                        >  1:15727/1‑58 (MQ=255)
                                                                  aGTATTCCTTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgtc  <  1:274446/90‑1 (MQ=255)
                                                                          |                                                                                 
TGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTC  >  NZ_CP009273/3357681‑3357836

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