Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,9511.G74.5% 297.0 / 106.4 102coding (65/1386 nt)yifKputative transporter YifK
Reads supporting (aligned to +/- strand):  ref base . (15/11);  new base G (34/42);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 126 total aligned reads displayed.

CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCT  >  NC_000913/3980902‑3980999
                                                   |                                                 
ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                    <  1:311617/50‑1 (MQ=255)
ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                    <  2:2054659/50‑1 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑g                                                 >  1:4321034/1‑50 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑g                                                 <  1:3717044/50‑1 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑g                                                 <  2:4867874/50‑1 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑g                                                 <  1:1337746/50‑1 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑g                                                 <  1:1136976/50‑1 (MQ=255)
  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gg                                                >  1:3829941/1‑50 (MQ=255)
  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gg                                                <  2:50756/50‑1 (MQ=255)
   gagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggg                                               <  1:2913417/50‑1 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggg                                              >  2:3647412/1‑50 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggg                                              >  2:2249762/1‑50 (MQ=255)
      cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggggg                                            <  1:2328840/50‑1 (MQ=255)
      cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggggg                                            >  1:663257/1‑50 (MQ=255)
      cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑gggggg                                            >  1:2761664/1‑50 (MQ=255)
       tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCTCCCT‑‑ggggggg                                           >  1:4086086/1‑49 (MQ=255)
       tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggggggg                                           >  1:1254296/1‑49 (MQ=255)
        aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggggggg                                           >  2:3293980/1‑48 (MQ=255)
        aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑ggggggg                                           >  2:2761018/1‑48 (MQ=255)
        aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCa                                          >  2:2634010/1‑50 (MQ=255)
         cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGc                                           <  2:3296667/48‑1 (MQ=255)
         cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCa                                          <  1:4456399/49‑1 (MQ=255)
          aGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCAc                                         >  1:4094364/1‑49 (MQ=255)
            cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                       >  2:2894814/1‑50 (MQ=255)
            cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                       <  2:4537508/50‑1 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                      >  2:1460649/1‑50 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                      >  2:1177230/1‑50 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCAtt                                     <  1:1361571/50‑1 (MQ=255)
              tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTg                                    <  2:952927/50‑1 (MQ=255)
              tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTg                                    >  1:4094963/1‑50 (MQ=255)
              tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTg                                    <  2:5090561/50‑1 (MQ=255)
                ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                                   <  1:4652474/50‑1 (MQ=255)
                ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                                   >  1:5269781/1‑50 (MQ=255)
                ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                    <  1:2606120/49‑1 (MQ=39)
                ggCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTgg                                   <  2:4654046/49‑1 (MQ=255)
                 gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                                  >  2:3835420/1‑50 (MQ=255)
                 gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                                  <  1:526763/50‑1 (MQ=255)
                  cTGGAAGCTCGACATATCGAACTCATTGCCCTGGGGGGGCACCATTGGCg                                 >  1:4460627/1‑50 (MQ=255)
                  cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                                 >  2:1266857/1‑50 (MQ=255)
                  cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                                 <  2:106456/50‑1 (MQ=255)
                   tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                                >  2:556172/1‑50 (MQ=255)
                    ggAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                               >  2:616946/1‑50 (MQ=255)
                    ggAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCg                              >  1:398228/1‑50 (MQ=255)
                    ggAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCg                              >  2:5279974/1‑50 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                               <  2:2194685/49‑1 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                              >  1:2142731/1‑50 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                              <  2:4292699/50‑1 (MQ=255)
                       aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGc                            >  2:5232799/1‑50 (MQ=255)
                       aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGc                            <  1:3198310/50‑1 (MQ=255)
                        gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGcc                           >  2:1674479/1‑50 (MQ=255)
                         cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                          <  2:3469713/50‑1 (MQ=255)
                         cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                          <  1:4829/50‑1 (MQ=255)
                           cGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                        >  1:1162022/1‑50 (MQ=255)
                            gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                       <  2:1719424/50‑1 (MQ=255)
                            gACATATCGAACTCAGCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                        <  2:3106090/49‑1 (MQ=38)
                            gACAGATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                       <  1:3464357/50‑1 (MQ=255)
                             aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                      <  2:4091208/50‑1 (MQ=255)
                             aCATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTa                     >  1:2081253/1‑50 (MQ=255)
                              cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                     <  1:2529128/50‑1 (MQ=255)
                              cATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTAt                    <  1:4993940/50‑1 (MQ=255)
                               atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                    >  2:2188048/1‑50 (MQ=255)
                               atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                    >  1:4054246/1‑50 (MQ=255)
                               atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                    <  1:5196079/50‑1 (MQ=255)
                                tatCGAACTCATCGCCCTGGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                    <  1:3142322/50‑1 (MQ=255)
                                tatCGAACTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATgg                  >  2:4544842/1‑50 (MQ=255)
                                 atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg                  >  2:3968971/1‑50 (MQ=255)
                                 atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg                  >  1:320703/1‑50 (MQ=255)
                                  tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                 >  2:531976/1‑50 (MQ=255)
                                  tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg                 <  1:2537937/50‑1 (MQ=255)
                                  tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg                  <  1:2875542/49‑1 (MQ=39)
                                  tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg                  <  2:2483979/49‑1 (MQ=39)
                                  tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg                  <  1:1034526/49‑1 (MQ=39)
                                   cGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                >  2:3414332/1‑50 (MQ=255)
                                   cGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg                >  2:1292261/1‑50 (MQ=255)
                                    gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg               <  1:493072/50‑1 (MQ=255)
                                    gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg               >  2:426718/1‑50 (MQ=255)
                                    gAACTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATggggg               <  2:4846369/49‑1 (MQ=255)
                                     aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc              <  2:3998801/50‑1 (MQ=255)
                                     aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc              <  2:5149049/50‑1 (MQ=255)
                                     aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc              >  1:3147026/1‑50 (MQ=255)
                                     aaCTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc             <  2:2866802/50‑1 (MQ=255)
                                      aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc             <  1:834293/50‑1 (MQ=255)
                                      aCTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccg            >  1:3285662/1‑50 (MQ=255)
                                        tCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc           >  2:503127/1‑50 (MQ=255)
                                        tCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc           >  2:1657003/1‑50 (MQ=255)
                                         cATCGCCCTT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc           >  1:1031973/1‑49 (MQ=39)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc          <  2:2593793/49‑1 (MQ=39)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc          <  2:4085751/49‑1 (MQ=39)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc          <  2:911544/49‑1 (MQ=39)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc          <  1:226453/49‑1 (MQ=39)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc           <  1:2547512/48‑1 (MQ=39)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa         <  2:4695083/50‑1 (MQ=255)
                                           tCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg        <  1:4297356/50‑1 (MQ=255)
                                           tCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg        >  1:376778/1‑49 (MQ=255)
                                            cGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAg        >  2:2557885/1‑49 (MQ=39)
                                             gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt       <  1:1792824/49‑1 (MQ=39)
                                             gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt       <  1:55935/49‑1 (MQ=39)
                                             gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa      <  2:378187/50‑1 (MQ=255)
                                             gCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAc     >  1:3533892/1‑50 (MQ=255)
                                                  t‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCt  >  2:2111687/1‑49 (MQ=255)
                                                   |                                                 
CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCT  >  NC_000913/3980902‑3980999

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