Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913897,3550TG26.8% 142.8 / 10.2 56V91G (GTG→GGG) potIputrescine ABC transporter membrane subunit PotI
Reads supporting (aligned to +/- strand):  ref base T (11/30);  new base G (15/0);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 4.75e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACC  >  NC_000913/897307‑897404
                                                |                                                  
cggcgACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGtggtgctggtg                                                   <  1:1530276/50‑1 (MQ=255)
cggcgACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGtggtgctggt                                                    <  1:1696653/49‑1 (MQ=255)
 ggcgACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGtggtgctggtgc                                                  <  1:7358212/50‑1 (MQ=255)
  gcgACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTgcg                                                 <  1:4071868/50‑1 (MQ=255)
  gcgACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGGGCTGGTgcg                                                 >  1:7238691/1‑50 (MQ=255)
   cgACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTgcgc                                               <  1:1454343/50‑1 (MQ=255)
    gACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCt                                              <  1:6228207/50‑1 (MQ=255)
    gACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCt                                              <  1:518181/50‑1 (MQ=255)
    gACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGGGCGCt                                              >  1:4021143/1‑50 (MQ=255)
    gACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGGGGTGGGGCGCt                                              >  1:5594583/1‑50 (MQ=255)
     aCGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCtt                                             <  1:6747558/50‑1 (MQ=255)
     aCGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGGGTTGGGGCGCtt                                             >  1:870187/1‑50 (MQ=255)
      cggcggcggCGATCCTCGGAACTATTGCGACGGTGGTGCTGGTGCGCttt                                            <  1:246145/50‑1 (MQ=255)
       ggcggcggcGATCCTCGGAACTATTGCGGCGGTGGTGGTGGGGCGCTTTg                                           >  1:1111426/1‑50 (MQ=39)
       ggcggcggcGATCCTCGGAACTATTGCGGCGGTGGGGCTGGGGGGCTTTg                                           >  1:1474187/1‑50 (MQ=38)
       ggcggcggcGATCCTCGGAACTATTGCGGCGGGGGTGCTGGGGGGCtttt                                           >  1:4498679/1‑49 (MQ=255)
       ggcggcggcGATCCTCGGAACTATTGCGGCGGCGGGGCTGGTGCGCtttc                                           >  1:1978169/1‑49 (MQ=255)
        gcggcggcgATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTgg                                          <  1:2929923/50‑1 (MQ=255)
          ggcggcgATCCTCGGAACTATTGCGGCGGTGGTTCGTGTGCGGTTTgg                                          >  1:5442371/1‑48 (MQ=37)
          ggcggcgATCCTCGGAACTATTGCGGCGGTGGTGCTGGGGCGCTTTgg                                          >  1:97112/1‑48 (MQ=255)
          ggcggcgATCCTCGGAACTATTGCGGCGGTGGGGCTGGGGCGCTTTgg                                          >  1:1801546/1‑48 (MQ=38)
           gcggcgATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAg                                       <  1:6201346/50‑1 (MQ=255)
           gcggcgATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAg                                       <  1:4138653/50‑1 (MQ=255)
            cggcgATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCGTTGGCAgg                                      <  1:640046/50‑1 (MQ=255)
            cggcgATCCTCGGAACTATTGCGGCGGTGGTGCTGGGGCGCCTTGGGAgg                                      >  1:2585208/1‑50 (MQ=38)
             ggcgATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGt                                     >  1:4848386/1‑50 (MQ=255)
             ggcgATCCTCGGAACTATTGCGGCGGTGGTGCTGGGGCGCTTTGGCAGGt                                     >  1:2108458/1‑50 (MQ=255)
               cgATCCTCGGAACTATTGCGGCGGTGGTGCTGGGGCGCTTTGGCAGGttt                                   >  1:5073630/1‑50 (MQ=255)
               cgATCCTCGGAACTATTGCGGCGGTGGGGCTGGGGCGCTTTTGCAGGttt                                   >  1:3985206/1‑50 (MQ=38)
                gATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTc                                  <  1:4835365/50‑1 (MQ=255)
                gATCCTCGGAACTATTGCGGCGGTGGGGCTGGGGCGCCTTCGGCGGGttt                                   >  1:2298474/1‑50 (MQ=19)
                 aTCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTcg                                 <  1:4718310/50‑1 (MQ=255)
                 aTCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTcg                                 >  1:3042400/1‑50 (MQ=255)
                 aTCCTCGGAACTATTGCGGCGGTGGGGCTGGTGCGCTTTGGCAGGTTTcg                                 >  1:1902867/1‑50 (MQ=255)
                  tCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTcgc                                <  1:4339965/50‑1 (MQ=255)
                   ccTCGGAACTATTGCGGCGGTGGTGCTGGGGCGCTTTTGCCGGTTTcgcg                               >  1:7277716/1‑50 (MQ=37)
                    cTCGGAACTATTGCGGCGGTGGGGCTGGTGCGCTTTTGCAGGTTTcgcg                               >  1:4366794/1‑49 (MQ=38)
                       ggAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGAt                            <  1:6552102/49‑1 (MQ=255)
                       ggAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGAt                            <  1:2928168/49‑1 (MQ=255)
                       ggAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGAt                            <  1:3061525/49‑1 (MQ=255)
                       ggAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATc                           <  1:6386393/50‑1 (MQ=255)
                       ggAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATc                           <  1:5601857/50‑1 (MQ=255)
                        gAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCa                          <  1:4206544/50‑1 (MQ=255)
                        gAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCa                          <  1:3088972/50‑1 (MQ=255)
                        gAACTATTGCGGCGGTGGAGCTTGGGCGCTTTGGCAGGTTTCGCGGGTc                           >  1:6457035/1‑49 (MQ=37)
                           cTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAAt                       <  1:7404820/50‑1 (MQ=255)
                              ttGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATgg                     <  1:1783541/49‑1 (MQ=255)
                                gcggcgGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCt                   <  1:3404648/49‑1 (MQ=255)
                                 cggcggTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCtt                  <  1:6143578/49‑1 (MQ=255)
                                  ggcggTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTg                >  1:5923451/1‑50 (MQ=255)
                                   gcggTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGc               <  1:1641634/50‑1 (MQ=255)
                                    cggTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGcc              >  1:5240907/1‑50 (MQ=255)
                                    cggTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGcc              <  1:6936632/50‑1 (MQ=255)
                                    cggTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGcc              <  1:3040736/50‑1 (MQ=255)
                                     ggtggtGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCt             >  1:243515/1‑50 (MQ=255)
                                      gtggtgCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCtt            <  1:4453120/50‑1 (MQ=255)
                                      gtggtgCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCtt            >  1:6067132/1‑50 (MQ=255)
                                      gtggtgCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCtt            <  1:3735710/50‑1 (MQ=255)
                                      gtggtgCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCtt            <  1:1719177/50‑1 (MQ=255)
                                       tggtgctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCttt           <  1:3791211/50‑1 (MQ=255)
                                        ggtgctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTa          >  1:4247537/1‑50 (MQ=255)
                                         gtgctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTa          <  1:5258521/49‑1 (MQ=255)
                                         gtgctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTa          <  1:2929247/49‑1 (MQ=255)
                                         gtgctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTAt         >  1:3139182/1‑50 (MQ=255)
                                         gtgctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTAt         <  1:2525292/50‑1 (MQ=255)
                                         gtgctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTAt         <  1:682369/50‑1 (MQ=255)
                                          tgctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTAt         >  1:2794341/1‑49 (MQ=255)
                                            ctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGa       <  1:5275179/49‑1 (MQ=255)
                                            ctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGa       <  1:3253222/49‑1 (MQ=255)
                                            ctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGAt      <  1:3263086/50‑1 (MQ=255)
                                            ctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGAt      <  1:1101439/50‑1 (MQ=255)
                                               gtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATc     <  1:6625961/48‑1 (MQ=255)
                                               gtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCa    <  1:1858516/49‑1 (MQ=255)
                                               gtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCAc   >  1:6341400/1‑50 (MQ=255)
                                               gtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCAc   <  1:6763817/50‑1 (MQ=255)
                                               gtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCAc   <  1:7080426/50‑1 (MQ=255)
                                               gtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCAc   <  1:7327484/50‑1 (MQ=255)
                                                tgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCAcc  <  1:2449136/50‑1 (MQ=255)
                                                tgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCAcc  >  1:4027265/1‑50 (MQ=255)
                                                tgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCAcc  <  1:7302301/50‑1 (MQ=255)
                                                tgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCAcc  <  1:2457512/50‑1 (MQ=255)
                                                |                                                  
CGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACC  >  NC_000913/897307‑897404

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