Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,070,4620TC23.8% 127.7 / 42.0 63intergenic (‑61/+53)yihR/yihSputative aldose 1‑epimerase YihR/sulfoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base T (30/18);  new base C (9/6);  total (39/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTT  >  NC_000913/4070413‑4070506
                                                 |                                            
tGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                              <  1:3485448/50‑1 (MQ=255)
tGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                              >  1:3486421/1‑50 (MQ=255)
tGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                              >  1:4335674/1‑50 (MQ=255)
  aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt                                            >  1:2450375/1‑50 (MQ=255)
  aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt                                            <  1:497065/50‑1 (MQ=255)
  aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                             <  1:1119647/49‑1 (MQ=255)
  aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                             <  1:4985142/49‑1 (MQ=255)
   ggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTg                                           >  1:1319998‑M2/29‑50 (MQ=255)
   ggCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                           >  1:5581563/1‑50 (MQ=255)
   ggCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                           <  1:5278314/50‑1 (MQ=255)
   ggCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                           <  1:994994/50‑1 (MQ=255)
    gctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTgg                                          >  1:3514537‑M2/28‑50 (MQ=255)
    gCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtttag                                          >  1:83249/1‑48 (MQ=255)
    gCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                          <  1:2035663/50‑1 (MQ=255)
    gCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                          <  1:3687370/50‑1 (MQ=255)
      tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                         >  1:1557791/1‑49 (MQ=255)
      tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                         >  1:3919760/1‑49 (MQ=255)
        caagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg                                         >  1:3445342‑M2/24‑47 (MQ=255)
        caagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg                                         >  1:258996‑M2/24‑47 (MQ=255)
         tGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                     <  1:2416291/50‑1 (MQ=255)
         aagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGga                                       >  1:3711977‑M2/23‑48 (MQ=255)
          agttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGagag                                    >  1:1957106‑M2/22‑50 (MQ=255)
           gttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAgg                                   >  1:6473029‑M2/21‑50 (MQ=255)
            ttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAggg                                  <  1:1193157‑M2/31‑1 (MQ=255)
             tgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt                                 >  1:2736283‑M2/19‑50 (MQ=255)
             aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                  <  1:2201200/49‑1 (MQ=255)
             aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                  <  1:2767262/49‑1 (MQ=255)
             aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                                 >  1:3775917/1‑50 (MQ=255)
               cagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGtt                                <  1:6581203‑M2/33‑1 (MQ=255)
               aCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                                 >  1:2970277/1‑48 (MQ=255)
               aCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTa                               <  1:5222571/50‑1 (MQ=255)
                cccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAg                              <  1:4288368/50‑1 (MQ=255)
                agaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAg                              <  1:265911‑M2/35‑1 (MQ=255)
                 gaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAgg                             <  1:6096446‑M2/36‑1 (MQ=255)
                 ccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                             >  1:4325891/1‑50 (MQ=255)
                  cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                            <  1:5098670/50‑1 (MQ=255)
                  cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                            <  1:2116434/50‑1 (MQ=255)
                  cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                            <  1:2507810/50‑1 (MQ=255)
                  accgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAggg                            <  1:3650776‑M2/37‑1 (MQ=255)
                  accgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAggg                            <  1:4497640‑M2/37‑1 (MQ=255)
                  accgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAggg                            >  1:4328052‑M2/14‑50 (MQ=255)
                    aaaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgag                          <  1:3100375/50‑1 (MQ=255)
                    aaaaCTGCACAGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgag                          <  1:4978853/50‑1 (MQ=255)
                     aaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga                         <  1:1733066/50‑1 (MQ=255)
                     aaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga                         <  1:4647474/50‑1 (MQ=255)
                       aCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                       <  1:2470685/50‑1 (MQ=255)
                       aCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                       >  1:3218284/1‑50 (MQ=255)
                        cTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg                      >  1:3975081/1‑50 (MQ=255)
                        cTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg                      >  1:3625391/1‑50 (MQ=255)
                           cACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGa                    >  1:3248125/1‑49 (MQ=255)
                             ccGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                 >  1:4148579/1‑50 (MQ=255)
                               ggACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc               >  1:3503128/1‑50 (MQ=255)
                               ggACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc               >  1:1867221/1‑50 (MQ=255)
                                gACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAccc              >  1:4096055/1‑50 (MQ=255)
                                   ggTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCa           >  1:1933812/1‑50 (MQ=255)
                                   ggTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCa           >  1:4180839/1‑50 (MQ=255)
                                   ggTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCa           >  1:5987693/1‑50 (MQ=255)
                                     tCCCCTCGCCCGTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAcc         <  1:6122447/50‑1 (MQ=255)
                                     tCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAcc         >  1:6686380/1‑50 (MQ=255)
                                      ccccTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAccc        >  1:2167312/1‑50 (MQ=255)
                                         cTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt     >  1:4595174/1‑50 (MQ=255)
                                         cTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt     >  1:2510982/1‑50 (MQ=255)
                                         cTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt     >  1:5640823/1‑50 (MQ=255)
                                            gCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt  >  1:281061/1‑50 (MQ=255)
                                            gCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt  >  1:1454265/1‑50 (MQ=255)
                                            gCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt  >  1:2608200/1‑50 (MQ=255)
                                             ccccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt  >  1:4384351/1‑49 (MQ=255)
                                               ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt  >  1:6528431/1‑47 (MQ=255)
                                               ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt  >  1:2604943/1‑47 (MQ=255)
                                               ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt  >  1:3655193/1‑47 (MQ=255)
                                               ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt  >  1:2830837/1‑47 (MQ=255)
                                                 |                                            
TGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTT  >  NC_000913/4070413‑4070506

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