Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,6560TG20.0% 180.5 / 32.2 65A192A (GCA→GCCyadLputative fimbrial protein YadL
Reads supporting (aligned to +/- strand):  ref base T (0/52);  new base G (13/0);  total (13/52)
Fisher's exact test for biased strand distribution p-value = 6.09e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCAC  >  NC_000913/151594‑151720
                                                              |                                                                 
aTGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGc                                                                 <  1:160452/65‑1 (MQ=255)
  gCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAg                                                               <  1:1053793/65‑1 (MQ=255)
   cATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGa                                                              <  1:1308901/65‑1 (MQ=255)
     tatataGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaa                                                            <  1:2267642/65‑1 (MQ=255)
      atataGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGCGaaaa                                                           >  1:2883616/1‑65 (MQ=255)
       tataGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAAc                                                          <  1:645291/65‑1 (MQ=255)
         taGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGGGAAAAAtt                                                        >  1:1404375/1‑65 (MQ=255)
         taGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCGGGAAAAtt                                                        >  1:1764022/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGGAGGAAAATTa                                                       >  1:2880612/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTa                                                       <  1:2904923/65‑1 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGGAAAAttt                                                       >  1:2684889/1‑64 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAACttt                                                       >  1:1384197/1‑64 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAACTTa                                                       >  1:1629744/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAACTTa                                                       >  1:478604/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAACTTa                                                       >  1:2311160/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAAATTa                                                       >  1:2273345/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCGGAAAACATa                                                       >  1:1181911/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAGACCTa                                                       >  1:1698580/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAACTTa                                                       >  1:394666/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAACTTa                                                       >  1:94756/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAAttt                                                       >  1:1432704/1‑64 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAATTa                                                       >  1:1036183/1‑65 (MQ=255)
          aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAATTa                                                       >  1:2881020/1‑65 (MQ=255)
            ttCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTAcc                                                     <  1:182072/65‑1 (MQ=255)
              ccTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGCGAGAACATACCCg                                                   >  1:458298/1‑65 (MQ=255)
               cTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTTCCCGc                                                  <  1:2879646/65‑1 (MQ=255)
                tGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCa                                                 <  1:2657848/65‑1 (MQ=255)
                tGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCa                                                 <  1:870974/65‑1 (MQ=255)
                  aaTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGt                                               <  1:2738245/65‑1 (MQ=255)
                  aaTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGt                                               <  1:2017915/65‑1 (MQ=255)
                  aaTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGt                                               <  1:2258272/65‑1 (MQ=255)
                  aaTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGt                                               <  1:1354351/65‑1 (MQ=255)
                       gTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACAt                                          <  1:353480/65‑1 (MQ=255)
                          aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                                      <  1:588117/65‑1 (MQ=255)
                          aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                                      <  1:1161665/65‑1 (MQ=255)
                          aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                                      <  1:278013/65‑1 (MQ=255)
                          aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                                      <  1:189618/65‑1 (MQ=255)
                          aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                                      <  1:135297/65‑1 (MQ=255)
                          aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGt                                      <  1:2611155/65‑1 (MQ=255)
                          aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTAGTGCAGAAAACTTACCCGCAGTGACATTGt                                      <  1:1632953/65‑1 (MQ=255)
                            tCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGGGAAAAATTACCCGCCGGGAAATTGttt                                    >  1:743746/1‑65 (MQ=255)
                            tCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAAATTACCCGCAGGGGCATTGttt                                    >  1:2257415/1‑65 (MQ=255)
                            tCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAATTACCCCCAGGGACATTGttt                                    >  1:853692/1‑65 (MQ=255)
                                ttagttaTAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACgg                                <  1:1570355/65‑1 (MQ=255)
                                 tagttaTAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGt                               <  1:2674415/65‑1 (MQ=255)
                                 tagttaTAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGt                               <  1:2677743/65‑1 (MQ=255)
                                 tagttaTAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGt                               <  1:2529636/65‑1 (MQ=255)
                                  agttaTAGGTTACTGTAAATGTGGCTGGGGGAGGAAAAATACCCCCCGGGGCATTCTTTTCgggt                              >  1:1694604/1‑63 (MQ=255)
                                      ataGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTAt                          <  1:2867281/65‑1 (MQ=255)
                                      ataGGTTACTGTAAATGTGGCTGGGGGAGAAAAAATACCCCCAGGGGCATTGTTTTAGGGttttt                           >  1:2349743/1‑64 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGGGGAAAAACTACCCCCCGGGGCCTTTTTTACGGGtttttttt                       >  1:1881768/1‑65 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGGAGGAAACTTACCCGCCGGGGCATTTTTTTTGGTtttttttt                       >  1:2090483/1‑65 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGGAGGAAAATTACCCGCAGCGGCATTGTTTTCGGTtttttttt                       >  1:1876625/1‑65 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGGAGAAAAATTATCCGCAGTGACGATGTTTACGGCtttttttt                       >  1:1709500/1‑65 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGGAGAAAAATTACCCGCCGGGGCATTTTTTTTGGTtttttttt                       >  1:1956279/1‑65 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGGAGAAAAATTACCCCCGGTGTCATTTTTTTTGGGtttttttt                       >  1:1385516/1‑65 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGGAGAAAAATTACCCCCCGGGGCATTGTTTTCGGGtttttttt                       >  1:808212/1‑65 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGGAGAAAAATTAACCCCCGTGGCATTGTTTTCGGGtttttttt                       >  1:1060275/1‑65 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGCAGGAAACATTCCCCCCGTGACAATGTTTTCGGGTTTCtttt                       >  1:950365/1‑65 (MQ=255)
                                         ggTTACTGTAAATGTGGCTGGGGCAGAAAAATTAACCCCAATGCCACTGTTTACGGTTTTAtttt                       >  1:78952/1‑65 (MQ=255)
                                             aCTGTAAATGTGGCTGGGGGGGAAAACTTACCCGCAGGGACACTTTTTTCGGGTttttttttttt                   >  1:934185/1‑64 (MQ=255)
                                             aCTGTAAATGTGGCTGGGGGGGAAAAATTACCCGCAGTGGCATTGTTTTTGGttttttttttttt                   >  1:2303432/1‑64 (MQ=255)
                                               tGTAAATGTGGCTGGGGGCGAAAAATTACCCGCCGCGACATTTTTTTTGGtttttttttttttta                 >  1:2608919/1‑62 (MQ=255)
                                               tGTAAATGTGGCTGGGGCCGAAAAATTACCCGCCGGGACAATGCTTTAGGGTTTTTTTTTTTTcc                  >  1:1341050/1‑65 (MQ=255)
                                                gTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCAc                <  1:656322/65‑1 (MQ=255)
                                                gTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCAc                <  1:2898430/65‑1 (MQ=255)
                                                  aaaTGTGGCTTGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa              <  1:207934/65‑1 (MQ=255)
                                                  aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa              <  1:1403673/65‑1 (MQ=255)
                                                  aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa              <  1:1277807/65‑1 (MQ=255)
                                                  aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa              <  1:2336555/65‑1 (MQ=255)
                                                  aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGa              <  1:2353854/65‑1 (MQ=255)
                                                  aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTCTACGGTTTTATTTTTTTCCACGa              <  1:1152761/65‑1 (MQ=255)
                                                   aaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGAc             <  1:593742/65‑1 (MQ=255)
                                                    aTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGAcc            <  1:1669419/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:1442333/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:1891688/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:846417/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:2495524/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:2516976/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:2172477/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:1109092/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:399933/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:1219812/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:2054612/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:2932237/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:1226025/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:1019299/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:1249561/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:131265/65‑1 (MQ=255)
                                                       tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:1920542/65‑1 (MQ=255)
                                                       tgGCTGGGTCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAgg         <  1:715278/65‑1 (MQ=255)
                                                            ggTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGc    <  1:2910120/65‑1 (MQ=255)
                                                            ggTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGc    <  1:1089248/65‑1 (MQ=255)
                                                              tGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCAc  <  1:1806454/65‑1 (MQ=255)
                                                              |                                                                 
ATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTGTTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCAC  >  NC_000913/151594‑151720

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