Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,336,7850AC60.0% 15.4 / 27.5 25G288G (GGT→GGGyfaLAIDA‑I family autotransporter adhesin YfaL/EhaC
Reads supporting (aligned to +/- strand):  ref base A (9/1);  new base C (0/15);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 4.89e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-05
Rejected as consensus: Frequency below/above cutoff threshold.

GTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC  >  NZ_CP009273/2336711‑2336874
                                                                          |                                                                                         
ttATTCCAAATAACCAGCTTTTTCCCGTCGGAAATAAAAAAGGAAACTCCACTAAAGCCGAGATAAAAAAACCCCCCCGCCGCAGaggaa                                                                            <  1:684830/89‑1 (MQ=255)
gTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGaggaa                                                                            >  1:244524/1‑90 (MQ=255)
     gCCAATACCCAGCTTTTTTCCGTCGGAAATACAAAAGGTAACTTCCCTAAACCCGAGATAAAAAAAGCCCCCCGCCGCAGAGGAAGGAcc                                                                       <  2:878378/90‑1 (MQ=255)
     gCCAATACCCAGCTTTTTCCCGCCGGAAATATAAAAGGTAACTTCACTTAAGCCGAGATACAAAAAGCCCCCCGCCGCAGAGGAAGGAcc                                                                       <  2:808122/90‑1 (MQ=255)
      cagaTACCCAGCTTTTTTCCGTCGCCAATACCAAAGGTACCTTCACTTAAGCCGAAATACAAAAACCCCCCCGCCGCAGAGGAAGGACCa                                                                      <  1:71112/87‑1 (MQ=255)
        aaTAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCAtc                                                                    >  1:194106/1‑90 (MQ=255)
        aaTAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCAtc                                                                    <  2:322589/90‑1 (MQ=255)
               aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg                                                             >  2:526468/1‑90 (MQ=255)
               aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg                                                             >  1:839304/1‑90 (MQ=255)
               aGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAg                                                             >  1:545532/1‑90 (MQ=255)
                   tttttCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCctg                                                         >  2:71112/1‑90 (MQ=255)
                    ttttCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCctgc                                                        <  2:347750/90‑1 (MQ=255)
                       tCCGTCGGAAATATCAAAGGTACCTTCACTTAAGCCGAGATACATAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCtgctgc                                                     <  2:320125/90‑1 (MQ=255)
                            cGGAAATAAAAAAGGTAACTTCCCTAAAGCCGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTAt                                                <  2:393589/90‑1 (MQ=255)
                                         ggTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGAc                                   >  1:225097/1‑90 (MQ=255)
                                                  aCTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGcc                          >  1:878516/1‑90 (MQ=255)
                                                   ttaaaGCCGAAAAAAATAACCCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCt                         <  1:481988/87‑1 (MQ=255)
                                                    taaaGCCGAAATACAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTc                        <  2:545532/88‑1 (MQ=255)
                                                         ccGAGATACATAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTc                   <  2:10738/90‑1 (MQ=255)
                                                         ccGAAATAAATAAGGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTc                   <  2:294595/90‑1 (MQ=255)
                                                         ccGAAATAAAAAAACCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTc                   <  1:70894/90‑1 (MQ=255)
                                                           gagaTACATAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTg                 <  2:225097/90‑1 (MQ=255)
                                                           gagaTACAAAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTg                 <  2:495035/90‑1 (MQ=255)
                                                                 cATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTa           >  1:64104/1‑90 (MQ=255)
                                                                   tAAAGCCCCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGc         <  2:653041/90‑1 (MQ=255)
                                                                          aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTAt                                                >  1:137621/1‑44 (MQ=255)
                                                                          aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTAt                                                <  2:137621/44‑1 (MQ=255)
                                                                          aCCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGtcatc  >  2:879586/1‑90 (MQ=255)
                                                                          |                                                                                         
GTATTGCCAATAACCAGCGTTTTTCCGTCGGCAATATCAAAGGTAACTTCACTTAAGCCGAGATACATAAAGCCACCCGCCGCAGAGGAAGGACCATCTCCATAGCCTGCTGCGCTATTGTTCTCATCGACTAACACGCCTCCGTTCTGGCTGTAGCTGTCATC  >  NZ_CP009273/2336711‑2336874

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: