Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,6560TG26.1% 86.4 / 32.5 46A192A (GCA→GCCyadLputative fimbrial protein YadL
Reads supporting (aligned to +/- strand):  ref base T (0/34);  new base G (12/0);  total (12/34)
Fisher's exact test for biased strand distribution p-value = 2.57e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGG  >  NC_000913/151594‑151713
                                                              |                                                          
aTGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGc                                                          <  1:542322/65‑1 (MQ=255)
       tataGTTCCTGAATCTGTACGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGgaaaa                                                  >  1:831163/1‑64 (MQ=255)
        ataGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGGGaaaaac                                                 >  1:1894864/1‑63 (MQ=255)
        ataGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGaaaatg                                                 >  1:1985850/1‑63 (MQ=255)
        ataGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGgaaaat                                                 >  1:2795637/1‑63 (MQ=255)
         taGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAAAtt                                                >  1:2904914/1‑65 (MQ=255)
         taGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAAAtt                                                >  1:2178755/1‑65 (MQ=255)
         taGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAACtt                                                >  1:421518/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTCGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAAATTa                                               >  1:1291065/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGGGAAAACATa                                               >  1:3121905/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGGGAAAAATTa                                               >  1:1241187/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGCGAAAAcgtt                                               >  1:1454553/1‑62 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAACTTa                                               >  1:3148263/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAAATTa                                               >  1:389933/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAAATTa                                               >  1:891952/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCGGGAAAATTa                                               >  1:130613/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCGGAAAACTTa                                               >  1:2597095/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCCGAGAAATTa                                               >  1:3157829/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAACTTa                                               >  1:312791/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAACATa                                               >  1:223012/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAATTa                                               >  1:170690/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAATTa                                               >  1:3450602/1‑65 (MQ=255)
              ccTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCg                                           <  1:2684273/65‑1 (MQ=255)
                tGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCa                                         <  1:3207302/65‑1 (MQ=255)
                tGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCa                                         <  1:2957732/65‑1 (MQ=255)
                tGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCa                                         <  1:1636184/65‑1 (MQ=255)
                  aaTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGt                                       <  1:3447427/65‑1 (MQ=255)
                      tGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACa                                   <  1:1938527/65‑1 (MQ=255)
                         aaGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTg                                <  1:2535961/65‑1 (MQ=255)
                          aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                               <  1:2902925/65‑1 (MQ=255)
                          aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                               <  1:1935926/65‑1 (MQ=255)
                          aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                               <  1:1073693/65‑1 (MQ=255)
                          aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                               <  1:1855690/65‑1 (MQ=255)
                            tCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCGGAAAACTTACCCGCAGTGACATTGttt                             >  1:2921130/1‑65 (MQ=255)
                            tCTGTTAGTTATAGGCTACTGTAAATGTGGCTGGGGCGGCAAAACTAACCCCTATGGCATTGttt                             >  1:368797/1‑65 (MQ=255)
                                 tagttaTAGGTCACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGt                        <  1:1491952/65‑1 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGGAGGAAAATTACCCCCAGGGACATTGTTTTTGGGTttttttc                >  1:1054032/1‑64 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGGAGAAAACTTACCCGCCGGGACCTTTTTTTCGGGtttttttt                >  1:1733518/1‑65 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGGAGAAAACTTACCCCCAGGGACATTGTTTTTGGTtttttttt                >  1:2782509/1‑65 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGCGGAAAAGATAACACCCGGGACATTGTTTGCGGTTTTAtttt                >  1:2915764/1‑65 (MQ=255)
                                          gTTACTGTAAATGTGGCTGGGGGAGAAAACTTACCCCCAGCGAAATTGTTTTACGtttttttttt               >  1:3410956/1‑65 (MQ=255)
                                          gTTACTGTAAATGTGGCTGGGGGAGAAAAATTACCCGCAGGGACATTGTTTTCGGGTtttttttt               >  1:1064893/1‑65 (MQ=255)
                                          gTTACTGTAAATGTGGCTGGGGCAGAAAACATACCCCCAGTGGCATTGTTTTTGGGTtttttttt               >  1:826921/1‑65 (MQ=255)
                                               tGTAAATGTGGCTGGGGGAGGAAAATTACCCGCGGGGACATTTTTTTCGGTTTTTTTTTTTTTCa          >  1:3210227/1‑65 (MQ=255)
                                               tGTAAATGTGGCTGGGGGAGAAAACTTACCCGCAGGGGCATTGTTTTCGGTTTTTTTTTTTTTCa          >  1:932173/1‑65 (MQ=255)
                                                gTAAATGTGGCTGGGGCAGAAAACTTACCCGCCGTGACATTGTTTTCGGTTTTATTTTTTTCCAc         >  1:40576/1‑65 (MQ=255)
                                                  aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa       <  1:404808/65‑1 (MQ=255)
                                                  aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa       <  1:423376/65‑1 (MQ=255)
                                                  aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa       <  1:979576/65‑1 (MQ=255)
                                                  aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa       <  1:2299773/65‑1 (MQ=255)
                                                  aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa       <  1:2301079/65‑1 (MQ=255)
                                                  aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa       <  1:2421753/65‑1 (MQ=255)
                                                     tgtgGCTGGGGGAGAAAACTTACCCGCCGTGGAATTGTTTTCGGGTTTTTTTTTTTTCAccacca    >  1:1146172/1‑65 (MQ=255)
                                                     tgtgGCTGGGGCAGAAAACTTACCCCCAGCGGCATTTTTTTTGGGTTTTTTTTTTTTCACGACCa    >  1:1249111/1‑65 (MQ=255)
                                                       ttgCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:2334961/63‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAGCTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:1822113/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTCCCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:2161811/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:3195169/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:2740256/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:1187589/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:1490514/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:72194/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:63037/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:2029956/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:2199038/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:3113032/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:3005789/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:2331017/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:3468715/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:2366961/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg  <  1:3303848/65‑1 (MQ=255)
                                                              |                                                          
ATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGG  >  NC_000913/151594‑151713

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