Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,562,0850CT37.5% 116.5 / 145.2 120A133T (GCA→ACA) ddpAputative D,D‑dipeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (31/44);  new base T (22/23);  total (53/67)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 124 total aligned reads displayed.

TTCACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGTTTT  >  NC_000913/1562036‑1562134
                                                 |                                                 
ttCACCGTATGTTCGTCGGGAGCATCAAACTTTAAATCTTTGGGAAATGc                                                   >  2:1194719/1‑50 (MQ=255)
  cACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                                 >  2:3583618/1‑50 (MQ=255)
  cACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCtt                                                 >  2:4303535/1‑50 (MQ=255)
  cACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCtt                                                 >  1:4380166/1‑50 (MQ=255)
  cACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCtt                                                 >  1:4860533/1‑50 (MQ=255)
   aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                                 <  2:4309275/49‑1 (MQ=255)
   aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCtt                                                 <  1:1925859/49‑1 (MQ=255)
    ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                               >  2:2703718/1‑50 (MQ=255)
    ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                               >  1:5263168/1‑50 (MQ=255)
    ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTct                                               <  1:4379064/50‑1 (MQ=255)
    ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTct                                               >  1:1266520/1‑50 (MQ=255)
    ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTct                                               <  2:1199937/50‑1 (MQ=255)
    ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTct                                               >  1:3528208/1‑50 (MQ=255)
     cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTctg                                              <  2:3768376/50‑1 (MQ=255)
      gTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgc                                             >  2:1610177/1‑50 (MQ=255)
      gTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTctgc                                             >  2:2825554/1‑50 (MQ=255)
       tATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgct                                            <  2:3023305/50‑1 (MQ=255)
       tATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTctgct                                            >  2:2701276/1‑50 (MQ=255)
        aTGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgctg                                           >  1:2581706/1‑50 (MQ=255)
         tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTctgctg                                           <  2:1781786/49‑1 (MQ=255)
         tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTctgctg                                           <  2:1529427/49‑1 (MQ=255)
         tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTgg                                          >  2:4770135/1‑50 (MQ=255)
           ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcc                                        <  1:2538716/50‑1 (MQ=255)
           ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGc                                         <  2:1031439/49‑1 (MQ=255)
           ttCGTCGGGAGCAACAATCTTTAAATCTTTGTGAAATGCTTCTGCTGGcc                                        <  2:3853180/50‑1 (MQ=255)
             cgtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGccc                                       <  2:4603345/49‑1 (MQ=255)
             cgtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGccc                                       <  1:2314287/49‑1 (MQ=255)
             cgtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGccc                                       <  2:4533310/49‑1 (MQ=255)
             cgtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGccc                                       <  2:2923132/49‑1 (MQ=255)
              gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCt                                     <  2:1282963/50‑1 (MQ=255)
              gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCt                                     <  2:4980748/50‑1 (MQ=255)
              gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCt                                     <  2:4783327/50‑1 (MQ=255)
              gtcgGGAGCATCAATCTTTAAATCGTTGGGAAATGCTTCTGCTGGCCCCt                                     <  2:3008012/50‑1 (MQ=255)
              gtcgGGAGCATCAATCTATAAATCTTTGGGAAATGCTTCTGCTGGCCCCt                                     <  2:401310/50‑1 (MQ=255)
                cgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                                   >  1:4517438/1‑50 (MQ=255)
                cgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                                   >  2:3436324/1‑50 (MQ=255)
                cgGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTgg                                   >  1:1198593/1‑50 (MQ=255)
                cgGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTgg                                   >  2:1008586/1‑50 (MQ=255)
                 gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                                   <  1:2276626/49‑1 (MQ=255)
                 gggAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGc                                  >  2:2077970/1‑50 (MQ=255)
                 gggAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCAgg                                   <  2:4071575/49‑1 (MQ=255)
                 gggAGCATCAAACTTTAAATCCTTGGGAAATGCTTCTGCTGGCCCCTgg                                   <  1:4623524/49‑1 (MQ=255)
                  ggAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGcc                                 <  2:5343422/50‑1 (MQ=255)
                  ggAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGcc                                 >  1:1637956/1‑50 (MQ=255)
                   gAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCg                                <  1:3172755/50‑1 (MQ=255)
                   gAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCg                                <  2:4340368/50‑1 (MQ=255)
                   gAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGcc                                 >  2:2860697/1‑49 (MQ=255)
                   gAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGcc                                 >  1:1491153/1‑49 (MQ=255)
                   gAGCATCAAGCTTTAAAGCATTGGGAAATGCTTCTGCTGGCCCCTGGCCg                                <  1:2436326/50‑1 (MQ=255)
                     gCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAt                              <  1:1133227/50‑1 (MQ=255)
                     gCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAt                              >  2:767336/1‑50 (MQ=255)
                     gCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAt                              <  1:2155730/50‑1 (MQ=255)
                     gCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAt                              >  1:60224/1‑50 (MQ=255)
                      cATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtt                             <  2:3171436/50‑1 (MQ=255)
                      cATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGAtt                             <  2:4694356/50‑1 (MQ=255)
                       aTCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttt                            >  2:4048383/1‑50 (MQ=255)
                        tCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGAtttt                           <  2:1001564/50‑1 (MQ=255)
                         cAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                           >  2:5153419/1‑49 (MQ=255)
                          aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttttt                          <  2:4526984/49‑1 (MQ=255)
                          aaTCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTta                         <  1:3606125/50‑1 (MQ=255)
                          aaTCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTta                         <  1:200773/50‑1 (MQ=255)
                           aTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                        >  1:723842/1‑50 (MQ=255)
                           aTCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTtag                        >  1:260401/1‑50 (MQ=255)
                             cTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTtag                        >  2:5245723/1‑48 (MQ=255)
                              tttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                     >  2:2834236/1‑50 (MQ=255)
                              tttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                     >  2:710702/1‑50 (MQ=255)
                              tttAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTtagtag                     >  2:4974460/1‑50 (MQ=255)
                                tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGcc                   >  2:4549216/1‑50 (MQ=255)
                                tAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTtagt                       <  1:860611/46‑1 (MQ=255)
                                tAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTtagt                       <  2:3657371/46‑1 (MQ=255)
                                tAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGcc                   >  1:4625223/1‑50 (MQ=255)
                                 aaaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGcc                   <  2:2541222/49‑1 (MQ=255)
                                 aaaTCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGcc                   <  2:1178489/49‑1 (MQ=255)
                                 aaaTCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg                  >  2:338933/1‑50 (MQ=255)
                                   aTCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt                <  1:4498196/50‑1 (MQ=255)
                                   aTCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt                <  1:931050/50‑1 (MQ=255)
                                   aTCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt                >  1:981609/1‑50 (MQ=255)
                                   aTCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt                <  1:3584818/50‑1 (MQ=255)
                                   aTCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt                <  1:3307755/50‑1 (MQ=255)
                                      tttGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa              >  2:3932144/1‑49 (MQ=255)
                                       ttGGGAAATGTTTCTGCTGGCCCCTGGCCGCTTTTTAGTAGCCGCTCaa             >  1:2581530/1‑49 (MQ=255)
                                        tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTc               <  1:4213888/46‑1 (MQ=255)
                                        tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa              <  2:1407446/47‑1 (MQ=255)
                                        tGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaa            <  2:3301403/49‑1 (MQ=255)
                                        tGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa             <  2:1714844/48‑1 (MQ=255)
                                         gggAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaaa           <  2:4837843/49‑1 (MQ=255)
                                         gggAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa             >  1:3381801/1‑47 (MQ=255)
                                         gggAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa             >  1:3305056/1‑47 (MQ=255)
                                          ggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaag          >  1:1878870/1‑49 (MQ=255)
                                          ggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaag          >  1:2073027/1‑49 (MQ=255)
                                          ggAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga         <  2:5358442/50‑1 (MQ=255)
                                            aaaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga         >  2:1942397/1‑48 (MQ=255)
                                            aaaTGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaa       <  1:5113933/50‑1 (MQ=255)
                                             aaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaag      >  1:4353090/1‑50 (MQ=255)
                                             aaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaag      <  1:3806279/50‑1 (MQ=255)
                                             aaTGCTTCTGCTAGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaag      <  2:2277736/50‑1 (MQ=255)
                                              aTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGt     <  2:4378946/50‑1 (MQ=255)
                                              aTGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGt     >  1:4950158/1‑50 (MQ=255)
                                                gCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGttt   >  1:4618661/1‑50 (MQ=255)
                                                 cTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGtttt  <  2:348381/50‑1 (MQ=255)
                                                 |                                                 
TTCACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGTTTT  >  NC_000913/1562036‑1562134

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