Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,6560TG33.3% 66.8 / 36.7 48A192A (GCA→GCC) yadLputative fimbrial protein YadL
Reads supporting (aligned to +/- strand):  ref base T (0/32);  new base G (16/0);  total (16/32)
Fisher's exact test for biased strand distribution p-value = 4.43e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCA  >  NC_000913/151591‑151723
                                                                 |                                                                   
gggATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGc                                                                   >  1:1724763/1‑68 (MQ=255)
   aTGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGgaga                                                                >  1:893038/1‑68 (MQ=255)
   aTGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGgaga                                                                >  1:53482/1‑68 (MQ=255)
   aTGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGgaga                                                                >  1:3477099/1‑68 (MQ=255)
    tGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCCGaa                                                               >  1:213746/1‑68 (MQ=255)
     gCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaa                                                              <  1:3505531/68‑1 (MQ=255)
     gCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaa                                                              <  1:2913795/68‑1 (MQ=255)
     gCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGACTGGTGCAGaaa                                                              <  1:2555472/68‑1 (MQ=255)
        tatataGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCCGAAAAca                                                           >  1:3155310/1‑67 (MQ=255)
            taGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTAcc                                                       <  1:1622900/68‑1 (MQ=255)
            taGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAACTTAcc                                                       >  1:3979340/1‑68 (MQ=255)
             aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGCGAAAACTTAccc                                                      >  1:3282807/1‑68 (MQ=255)
             aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGGAAAATTAacc                                                      >  1:3117425/1‑68 (MQ=255)
             aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAACTTAccc                                                      >  1:1612085/1‑68 (MQ=255)
             aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAAATTAccc                                                      >  1:302416/1‑68 (MQ=255)
             aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCGGAAAAATTAccc                                                      >  1:1667694/1‑68 (MQ=255)
             aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAACTTAccc                                                      >  1:760629/1‑68 (MQ=255)
             aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAATTAccc                                                      >  1:2097935/1‑68 (MQ=255)
             aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAATTAccc                                                      >  1:1046741/1‑68 (MQ=255)
                tCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAATTAcccccc                                                   >  1:1092262/1‑65 (MQ=255)
                  cTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGt                                                 <  1:2223429/68‑1 (MQ=255)
                  cTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGt                                                 <  1:2612414/68‑1 (MQ=255)
                    gAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGa                                               <  1:1869894/68‑1 (MQ=255)
                     aaTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGAc                                              <  1:1246440/68‑1 (MQ=255)
                     aaTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGAc                                              <  1:265181/68‑1 (MQ=255)
                       tCTGTAAGTCTGTTAGTTNTAGGTTACTGTAAATGTGGCTGGGGCAGGAAAAATACCCCCCGGGACAt                                            >  1:3091325/1‑68 (MQ=255)
                       tCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACAt                                            <  1:2871142/68‑1 (MQ=255)
                         tGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTg                                          <  1:432960/68‑1 (MQ=255)
                          gTAAGTCTGTTAGTTGTAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                                         <  1:800878/68‑1 (MQ=255)
                          gTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                                         <  1:1277726/68‑1 (MQ=255)
                          gTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                                         <  1:2491843/68‑1 (MQ=255)
                          gTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                                         <  1:840997/68‑1 (MQ=255)
                          gTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                                         <  1:2947900/68‑1 (MQ=255)
                             aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAACTTACCCGCAGTGACATTGtttt                                      >  1:178886/1‑67 (MQ=255)
                             aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAATTACCCGCAGTGACAATGtttt                                      >  1:961205/1‑67 (MQ=255)
                               tCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGGGAAAAATTAACCGCAGGGGCATATTTTatg                                    >  1:3386085/1‑66 (MQ=255)
                               tCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAAATTACCCCCAGGGACATTGtttttg                                    >  1:3310324/1‑65 (MQ=255)
                                 tGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACAGt                                  <  1:703338/68‑1 (MQ=255)
                                    tagttaTAGGTTACTGTAAATGTGGCTGGGGCAGAAAACTTACCCCCCGGGGCATTGTTTTCGGgttt                               >  1:2871607/1‑68 (MQ=255)
                                          taGGTTACTGTAAATGTGGCTGGGGCAGAAAAAATACGCGCCGGGACATTGTTTATGGGTtttttttt                         >  1:2178569/1‑68 (MQ=255)
                                            ggTTACTGTAAATGTGGCTGGGGGGGAAAACTTACCCGCCGTGGCATTGTTTTCGGGTtttttttttt                       >  1:3514455/1‑68 (MQ=255)
                                            ggTTACTGTAAATGTGGCTGGGGGAGGAAAATTACCCGCAGGGGCCTTTTTTTTGGGTtttttttttt                       >  1:519480/1‑68 (MQ=255)
                                            ggTTACTGTAAATGTGGCTGGGGCGGGAAAATTACCCCCCGGGGCCTTTTTTTCGGGTtttttttttt                       >  1:868300/1‑68 (MQ=255)
                                            ggTTACTGTAAATGTGGCTGGGGCAGAAAACTTACCCGCCGGGGCATTGTTTACGGTTTTAttttttt                       >  1:2373977/1‑68 (MQ=255)
                                            ggTTACTGTAAATGTGGCTGGGGCAGAAAACTTACCCCCCGGGGCATTGTTTTCGGtttttttttttt                       >  1:1886851/1‑68 (MQ=255)
                                            ggTTACTGTAAATGTGGCTGGGGCAGAAAAATTACCCCCAGCGGCATTTTTTTCGGGTtttttttttt                       >  1:298456/1‑68 (MQ=255)
                                                aCTGTAAATGTGGCTGGGGGAGAAAACTTACCCCCCGGGACATTGTTTTCGGGTTTTTTTTTTTTCAc                   >  1:3247393/1‑68 (MQ=255)
                                                aCTGTAAATGTGGCTGGGGGAGAAAACTTACCCCCAGGGGGATTGTTTTCGGTTTTATTTTTTTTCAc                   >  1:3855343/1‑68 (MQ=255)
                                                aCTGTAAATGTGGCTGGGGCAGAAAAATTACCCCCCGTGACATTGTTTTCGGGTTTTTTTTTTTTCAc                   >  1:3402741/1‑68 (MQ=255)
                                                  tGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa                 <  1:771838/68‑1 (MQ=255)
                                                  tGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa                 <  1:677104/68‑1 (MQ=255)
                                                  tGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa                 <  1:2859333/68‑1 (MQ=255)
                                                  tGTAAATGTGGCTGGGGCAGAAAACTTACCCGCAGTGAAATTGTTTACGGTTTTTTTTTTTTTCACGa                 >  1:1291788/1‑68 (MQ=255)
                                                  tGTAAATGTGGCTGGGGCAGAAAACTTACCCGCAGGGGCATTGTTTTCGGGTTTTTTTTTTTTCACGa                 >  1:2099172/1‑68 (MQ=255)
                                                   gTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGAc                <  1:3555486/68‑1 (MQ=255)
                                                    tAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGAcc               <  1:3237761/68‑1 (MQ=255)
                                                       aTGTGGCTGGTGCAGAAACCTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg            <  1:1376378/68‑1 (MQ=255)
                                                       aTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGTCCAgg            <  1:2409764/68‑1 (MQ=255)
                                                       aTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg            <  1:1888501/68‑1 (MQ=255)
                                                       aTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg            <  1:3729832/68‑1 (MQ=255)
                                                       aTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg            <  1:3160377/68‑1 (MQ=255)
                                                       aTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg            <  1:867774/68‑1 (MQ=255)
                                                       aTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg            <  1:2573792/68‑1 (MQ=255)
                                                       aTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg            <  1:640814/68‑1 (MQ=255)
                                                       aTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg            <  1:1732225/68‑1 (MQ=255)
                                                       aTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg            <  1:2607774/68‑1 (MQ=255)
                                                       aTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTGTCCACGGCCAgg            <  1:2543336/68‑1 (MQ=255)
                                                        tgtgGCTGGGGCAGAAAACTTACCCGCAGGGACATTGTTTTTGGGTTTTTTTTTTTTCACGACCAggg           >  1:1924792/1‑67 (MQ=255)
                                                               ggTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACt    <  1:1820552/68‑1 (MQ=255)
                                                                 tGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCa  <  1:1936965/68‑1 (MQ=255)
                                                                 tGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCa  <  1:3441196/68‑1 (MQ=255)
                                                                 |                                                                   
GGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCACTCA  >  NC_000913/151591‑151723

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: ‑