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

TTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCA  >  NZ_CP009273/3357671‑3357842
                                                                                    |                                                                                       
ttCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTgg                                                                                    <  2:278249/90‑1 (MQ=255)
          gggaagcaTCAATTTAGTGCTTAAAAAAATAATGAAGACGAGGGTGTTCTTTTCTACGTCAACATCAGTATCCTTTTTGGAAGAAAAAac                                                                          <  2:328536/86‑1 (MQ=255)
          gggCAGCATCACCTTAGTCCTTAAAAAAATAAAAAAGACGATGGTGATCTTTTCTACGTAAACATCAGTATTCCTTTTGGAAGACAAAac                                                                          <  2:345431/89‑1 (MQ=255)
             cagcaTCAACTTAGTGCTTCAGAAAAAAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACAcc                                                                       <  1:180249/90‑1 (MQ=255)
               gcaTCACCTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTCCGTCAACATCATTATTCCATTTGGCAGATCAAACCCCGc                                                                     <  2:380888/90‑1 (MQ=255)
               gcaTCAACTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTAAACATCAGTATTCCATTTGGCAGATAAAACACCGc                                                                     <  2:171030/90‑1 (MQ=255)
               gcaTCAACTTAGTGCTTAAAAAAATAATGAAGCCGATGGTGATCTGTTCTACGCCAACATCAGTATTCCTTTTGGCAGATCAACCACCGc                                                                     <  2:169972/90‑1 (MQ=255)
                 aTCAACTAAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTCCGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCca                                                                   <  2:63748/90‑1 (MQ=255)
                  tCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCcac                                                                  <  1:63947/90‑1 (MQ=255)
                  tCAACTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACCCCGCcac                                                                  <  2:427303/90‑1 (MQ=255)
                  tCAACTTAGTGCTTAAAAAAAAAATGAAGACGATGGTGATCTTTTCTACGTCAACATCAGTATTCCTTTTGGCAGATAAAACACCGCcac                                                                  <  2:67899/90‑1 (MQ=255)
                    aaCTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCCacac                                                                <  1:36565/90‑1 (MQ=255)
                    aaCTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCacac                                                                >  2:68858/1‑90 (MQ=255)
                    aaCTTAGTGCTTCAGAAAATAATGAAGACGAGGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATAAAACACCGCCacac                                                                <  1:557385/90‑1 (MQ=255)
                    aaCTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTTTTCTACGTCAACATCAGTATTCCTTTTGGCAGATAAAACACCGCCacac                                                                <  2:365648/90‑1 (MQ=255)
                    aaCTTAGTGCTTAAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCacac                                                                <  2:151269/90‑1 (MQ=255)
                    aaCTTAGTCCTTAAAAAAAAAATGAAGACGATGGTGATTTGTTCTCCGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCCacac                                                                <  2:30551/90‑1 (MQ=255)
                     aCTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCCACACt                                                               <  2:30498/90‑1 (MQ=255)
                      cTTAGTGCTTAAGAAAAAAATGAAGACGATGGTGCTCTTTTCTCCGTCAACACCATTATTCATTTTGGCAAACCAACCACCGCCACACTg                                                              <  2:344433/90‑1 (MQ=255)
                          gTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGtata                                                          >  2:269658/1‑90 (MQ=255)
                          gTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGtata                                                          >  1:412747/1‑90 (MQ=255)
                           tgtttCAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTAAACATCAGTTTTCCTTTTGGCAGATAAAACCCCGCCACACTGTATAc                                                         <  2:111307/87‑1 (MQ=255)
                                    aaaTAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGa                                                >  2:97128/1‑90 (MQ=255)
                                      ataataAAGACGAGGGTGATTTGTTCTCCGTCAACATCATTATTCCTTTTGGCAAATCAAACACCGCCACACTGTATACTCGCCATGACg                                              <  2:86911/90‑1 (MQ=255)
                                      aaaatAAAGAGGATGGTGTTTTTTTCTCCGCAAACACCAGTATTCTTTTTGGAAAAAAAAACCCCGCCACACTGTATACTCGCCATGACg                                              <  2:564266/88‑1 (MQ=255)
                                                  aTGGTGTTCTGTTCTACGTAAACATCAGTATTCCTTTTGGAAAACCAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAAccc                                  <  2:295095/90‑1 (MQ=255)
                                                   gggTGATCTTTTCTACGTCAACATCAGTTTTCCTTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCa                                 <  1:557900/89‑1 (MQ=255)
                                                       gATCTTTTCTACGCCAACATCAGTATCCCTTTTGGCAGATAAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTAt                             <  1:550865/90‑1 (MQ=255)
                                                            tttCTACGTCAAAATCATTATTCCTTTTGGCAGACAAACCACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTAc                        <  1:381658/89‑1 (MQ=255)
                                                            gTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTAc                        >  2:212622/1‑90 (MQ=255)
                                                            gTTCTACGCCAACATCAGTATTCCTTTTGGCAAATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTAc                        <  2:239251/90‑1 (MQ=255)
                                                              tCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTg                      >  2:6457/1‑90 (MQ=255)
                                                                 ccGTAAACATCAGTTTTCTTTTTGGAAGAAAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGgtg                   <  1:323722/89‑1 (MQ=255)
                                                                 acatCAACAATAGTGTCCCTTTTGACAGAGCAACCTCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGgtg                   <  2:207496/87‑1 (MQ=255)
                                                                 aCGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGgtg                   >  2:400270/1‑90 (MQ=255)
                                                                 aCGCCAAAATCGGTTTTCCTTTTGGCAACTAAATCACCGCCCCACTGTATACTCGCCATGACGATCATAAAACCCCCTATGGTACTGgtg                   <  2:388484/90‑1 (MQ=255)
                                                                      aaCATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATg              >  2:485255/1‑90 (MQ=255)
                                                                             gTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgtcg       >  2:72886/1‑90 (MQ=255)
                                                                             gTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgtcg       >  1:154264/1‑90 (MQ=255)
                                                                                  ccATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCa  >  1:215982/1‑90 (MQ=255)
                                                                                    |                                                                                       
TTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCA  >  NZ_CP009273/3357671‑3357842

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