Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,342,6160GA48.2% 16.3 / 187.5 114V518V (GTC→GTTyfaLputative autotransporter adhesin YfaL
Reads supporting (aligned to +/- strand):  ref base G (25/34);  new base A (31/24);  total (56/58)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 116 total aligned reads displayed.

CAAGATCTACCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTGGAGAAA  >  NC_000913/2342567‑2342664
                                                 |                                                
cAAGATCTACCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTg                                                  <  1:1180240/50‑1 (MQ=255)
 aaGATCTACCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGa                                                 >  1:4668491/1‑50 (MQ=255)
  aGATCTACCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTAAc                                                >  1:895442/1‑50 (MQ=255)
  aGATCTACCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTAAc                                                >  1:1268113/1‑50 (MQ=255)
   gATCTACCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACt                                               <  1:269228/50‑1 (MQ=255)
   gATCTACCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACt                                               <  1:327205/50‑1 (MQ=255)
   gATCTACCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTAACt                                               >  1:8391502/1‑50 (MQ=255)
    aTCTACCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACtt                                              >  1:2110568/1‑50 (MQ=255)
     tCTACCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACtt                                              >  1:2577547/1‑49 (MQ=255)
      cTACCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTc                                             >  1:1038747/1‑49 (MQ=255)
       tACCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTc                                             >  1:2326350/1‑48 (MQ=255)
        aCCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTAACTTCAAg                                          >  1:3861778/1‑50 (MQ=255)
         ccACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCaa                                           <  1:5269206/48‑1 (MQ=255)
         ccACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAgg                                         >  1:7951641/1‑50 (MQ=255)
         ccACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAgg                                         <  1:5963250/50‑1 (MQ=255)
         ccACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAgg                                         <  1:7844946/50‑1 (MQ=255)
         ccACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTAACTTCa                                            <  1:7224586/47‑1 (MQ=255)
         ccACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTAACTTCAAgg                                         >  1:5818961/1‑50 (MQ=255)
         ccACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTAACTTCAAg                                          <  1:1325813/49‑1 (MQ=255)
          cacaTCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGc                                        <  1:4605515/50‑1 (MQ=255)
          cacaTCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGc                                        <  1:6468514/50‑1 (MQ=255)
          cacaTCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGc                                        >  1:3589857/1‑50 (MQ=255)
          cacaTCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGc                                        >  1:5101277/1‑50 (MQ=255)
           acaTCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGcc                                       <  1:7590753/50‑1 (MQ=255)
           acaTCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGc                                        >  1:2326396/1‑49 (MQ=255)
            caTCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCa                                      >  1:1015071/1‑50 (MQ=255)
            caTCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCa                                      >  1:6522525/1‑50 (MQ=255)
            caTCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCa                                      >  1:2953336/1‑50 (MQ=255)
             aTCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCAt                                     >  1:7770196/1‑50 (MQ=255)
              tCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATc                                    <  1:6932986/50‑1 (MQ=255)
              tCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATc                                    <  1:3352450/50‑1 (MQ=255)
              tCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATc                                    >  1:7541077/1‑50 (MQ=255)
              tCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATc                                    <  1:757099/50‑1 (MQ=255)
              tCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATc                                    <  1:492038/50‑1 (MQ=255)
              tCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATc                                    >  1:4201913/1‑50 (MQ=255)
              tCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATc                                    <  1:1406395/50‑1 (MQ=255)
               ccTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCg                                   >  1:2175566/1‑50 (MQ=255)
                cTGACTGCCGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCg                                   <  1:7270812/49‑1 (MQ=255)
                cTGACTGCCGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGt                                  >  1:3614176/1‑50 (MQ=255)
                  gACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTg                                >  1:8432450/1‑50 (MQ=255)
                  gACTGCCGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTg                                >  1:2621946/1‑50 (MQ=255)
                   aCTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGtt                                 <  1:2661746/48‑1 (MQ=255)
                   aCTGCCGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGtt                                 <  1:1024535/48‑1 (MQ=255)
                    cTGCCGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGTa                              >  1:3971354/1‑50 (MQ=255)
                      gCCGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGTat                             >  1:6628238/1‑49 (MQ=255)
                       ccGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTatga                           <  1:1978548/50‑1 (MQ=255)
                       ccGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTat                             <  1:1039000/48‑1 (MQ=255)
                       ccGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGTatg                            <  1:2485243/49‑1 (MQ=255)
                       ccGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGTat                             <  1:1846811/48‑1 (MQ=255)
                        cGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTatgat                          >  1:7174434/1‑50 (MQ=255)
                        cGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTatgat                          >  1:3666468/1‑50 (MQ=255)
                        cGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGTatgat                          <  1:5249334/50‑1 (MQ=255)
                        cGATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGTatgat                          >  1:7870358/1‑50 (MQ=255)
                         gATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTatgatg                         <  1:2872855/50‑1 (MQ=255)
                         gATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTatgatg                         >  1:6860508/1‑50 (MQ=255)
                         gATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTatgatg                         <  1:4749239/50‑1 (MQ=255)
                         gATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTatgatg                         <  1:6798801/50‑1 (MQ=255)
                         gATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTatgatg                         <  1:6130036/50‑1 (MQ=255)
                         gATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGTatgatg                         >  1:4727602/1‑50 (MQ=255)
                         gATAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGAatgatg                         >  1:4068225/1‑50 (MQ=255)
                          aTAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGTATGAtgt                        <  1:1883934/50‑1 (MQ=255)
                           tAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGTATGAtgtg                       <  1:4874766/50‑1 (MQ=255)
                           tAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGTATGAtgtg                       >  1:463718/1‑50 (MQ=255)
                           tAACCGTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGTATGAtgtg                       <  1:7799225/50‑1 (MQ=255)
                             aCCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTATGATGTGcc                     <  1:8394517/50‑1 (MQ=255)
                              ccGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTATGATGTGCCg                    >  1:3514402/1‑50 (MQ=255)
                                gTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTATGATGTGCCGct                  <  1:7525360/50‑1 (MQ=255)
                                gTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGTATGATGTGCCGct                  <  1:1581659/50‑1 (MQ=255)
                                gTTCCGCTGCTGCCACTAACTTCAAGGCCATCGTTGTATGATGTGCCGct                  <  1:3253143/50‑1 (MQ=255)
                                 ttCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTATGATGTGCCGctc                 >  1:5421129/1‑50 (MQ=255)
                                  tCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTATGATGTGCCGctct                >  1:6937226/1‑50 (MQ=255)
                                   ccGCTGCTGCCACTGACTTCAAGGCCATCGTTGTATGATGTGCCGctct                <  1:6755840/49‑1 (MQ=255)
                                   ccGCTGCTGCCACTGACTTCAAGGCCATCGTTGTATGATGTGCCGctct                <  1:1832605/49‑1 (MQ=255)
                                   ccGCTGCTGCCACTGACTTCAAGGCCATCGTTGTATGATGTGCCGctct                <  1:3784916/49‑1 (MQ=255)
                                   ccGCTGCTGCCACTAACTTCAAGGCCATCGTTGTATGATGTGCCGctct                <  1:1029275/49‑1 (MQ=255)
                                   ccGCTGCTGCCACTAACTTCAAGGCCATCGTTGTATGATGTGCCGctct                <  1:6286635/49‑1 (MQ=255)
                                    cGCTGCTGCCACTGACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTg               <  1:5770706/49‑1 (MQ=255)
                                     gctgctGCCACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGcc             <  1:8467668/50‑1 (MQ=255)
                                     gctgctGCCACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGc              >  1:7683075/1‑49 (MQ=255)
                                      ctgctgCCACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGc              >  1:7578462/1‑48 (MQ=255)
                                       tgctgcCACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAg           >  1:2417881/1‑50 (MQ=255)
                                         ctgcCACTGACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAg           <  1:4572355/48‑1 (MQ=255)
                                         ctgcCACTGACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGt         <  1:7650637/50‑1 (MQ=255)
                                         ctgcCACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAgg          <  1:7659762/49‑1 (MQ=255)
                                         ctgcCACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGt         <  1:8389798/50‑1 (MQ=255)
                                          tgcCACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTg        >  1:27403/1‑50 (MQ=255)
                                          tgcCACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTg        >  1:7338429/1‑50 (MQ=255)
                                           gcCACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTg        >  1:5698790/1‑49 (MQ=255)
                                            ccACTGACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTg        <  1:757960/48‑1 (MQ=255)
                                            ccACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTg        >  1:3759219/1‑48 (MQ=255)
                                            ccACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTg        >  1:5119802/1‑48 (MQ=255)
                                             cACTGACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTg        >  1:2941892/1‑47 (MQ=255)
                                             cACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTg        >  1:4269530/1‑47 (MQ=255)
                                             cACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTg        >  1:6991520/1‑47 (MQ=255)
                                             cACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTGgag     <  1:5331728/50‑1 (MQ=255)
                                             cACTAACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTGgag     <  1:2674843/50‑1 (MQ=255)
                                              aCTGACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTGgag     >  1:4972485/1‑49 (MQ=255)
                                               cTGACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTGgaga    >  1:5754523/1‑49 (MQ=255)
                                               cTGACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTGGAGaa   <  1:4482293/50‑1 (MQ=255)
                                                tGACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTGGAGaaa  >  1:7515691/1‑50 (MQ=255)
                                                 |                                                
CAAGATCTACCACATCCTGACTGCCGATAACCGTTCCGCTGCTGCCACTGACTTCAAGGCCATCGTTGTATGATGTGCCGCTCTGCCAGGTGGAGAAA  >  NC_000913/2342567‑2342664

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