Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,070,4620TC28.3% 178.0 / 95.1 113intergenic (‑61/+53)yihR/yihSputative aldose 1‑epimerase YihR/sulfoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base T (34/47);  new base C (21/11);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 3.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 125 total aligned reads displayed.

TGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCA  >  NC_000913/4070413‑4070508
                                                 |                                               
tGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                                 >  2:4224591/1‑50 (MQ=255)
tGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                                 <  1:279725/50‑1 (MQ=255)
 gAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                >  2:3137301/1‑50 (MQ=255)
  gggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCctt                                               >  1:3179133‑M2/30‑50 (MQ=255)
  gggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCctt                                               <  1:5734429‑M2/21‑1 (MQ=255)
  gggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCctt                                               >  1:3807565‑M2/30‑50 (MQ=255)
  gggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCctt                                               <  1:1109976‑M2/21‑1 (MQ=255)
  aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                <  1:3562819/49‑1 (MQ=255)
   ggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTg                                             >  2:1459628‑M2/29‑50 (MQ=255)
    gctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTgg                                            >  2:1416370‑M2/28‑50 (MQ=255)
    gCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                            <  1:4384093/50‑1 (MQ=255)
     ctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg                                           >  2:1072216‑M2/27‑50 (MQ=255)
     ctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg                                           >  1:4539803‑M2/27‑50 (MQ=255)
     cTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                           <  1:5149063/50‑1 (MQ=255)
      tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCACGCCCCTTTggg                                           >  1:1469182/1‑49 (MQ=255)
       gcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg                                           >  2:5496144‑M2/25‑48 (MQ=255)
       gcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg                                           >  1:4416053‑M2/25‑48 (MQ=255)
       gATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                           >  1:1821485/1‑48 (MQ=255)
        caagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTgggg                                          <  1:743489‑M2/25‑1 (MQ=255)
        aTGTCAAACCCAAAAACTGCTCCGGACGGTCCCCTCGCCCCTTTGGGgag                                        <  2:3486940/50‑1 (MQ=255)
        aTGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgag                                        <  2:5526767/50‑1 (MQ=255)
        aTGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgag                                        <  1:3551624/50‑1 (MQ=255)
        aTGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgag                                        <  1:1867989/50‑1 (MQ=255)
         tGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                       <  2:4407643/50‑1 (MQ=255)
         tGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                       <  1:2741259/50‑1 (MQ=255)
         tGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                         >  2:1438074/1‑48 (MQ=255)
         tGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                         >  1:1124340/1‑48 (MQ=255)
         aagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGga                                         >  1:4733763‑M2/23‑48 (MQ=255)
          agttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGagag                                      <  2:1072501‑M2/29‑1 (MQ=255)
           tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                                      <  2:634311/49‑1 (MQ=255)
           tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                                      <  2:3162040/49‑1 (MQ=255)
           tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                     >  2:5419907/1‑50 (MQ=255)
           tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                     <  1:4019592/50‑1 (MQ=255)
           tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                     >  2:5236909/1‑50 (MQ=255)
           gttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAgg                                     <  2:443662‑M2/30‑1 (MQ=255)
           gttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAgg                                     <  2:4756310‑M2/30‑1 (MQ=255)
           gttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAgg                                     <  1:2314067‑M2/30‑1 (MQ=255)
            ttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAggg                                    >  1:1411786‑M2/20‑50 (MQ=255)
            cAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                    >  1:5839500/1‑50 (MQ=255)
             aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                    <  2:869169/49‑1 (MQ=255)
             aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                                   >  2:4880676/1‑50 (MQ=255)
             aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                                   <  1:1142336/50‑1 (MQ=255)
              aaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                  <  1:91652/50‑1 (MQ=255)
               cagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGtt                                  <  2:5203458‑M2/33‑1 (MQ=255)
               cagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt                                   >  2:3509031‑M2/17‑48 (MQ=255)
               aCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTa                                 <  2:2402160/50‑1 (MQ=255)
               aCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTa                                 <  2:1765102/50‑1 (MQ=255)
               aCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTa                                 <  1:1799953/50‑1 (MQ=255)
                cccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGAAGAGGGTTAg                                <  1:74794/50‑1 (MQ=255)
                agaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAg                                <  1:1419013‑M2/35‑1 (MQ=255)
                 gaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAgg                               <  1:5765111‑M2/36‑1 (MQ=255)
                 ccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                               <  1:5032647/50‑1 (MQ=255)
                  cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                              <  2:407162/50‑1 (MQ=255)
                  cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                              <  2:4412235/50‑1 (MQ=255)
                  cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                              <  2:4861838/50‑1 (MQ=255)
                   ccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt                             <  2:893577‑M2/38‑2 (MQ=255)
                   aaaaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                               >  1:3413273/1‑48 (MQ=255)
                   aaaaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGga                             <  1:726617/50‑1 (MQ=255)
                    cgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTtagggtg                            <  2:343967‑M2/39‑3 (MQ=255)
                    cgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt                             >  2:487512‑M2/12‑48 (MQ=255)
                    cgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt                             >  1:1066655‑M2/12‑48 (MQ=255)
                    cgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt                             >  2:1678946‑M2/12‑48 (MQ=255)
                    aaaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgag                            <  2:3225180/50‑1 (MQ=255)
                    aaaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgag                            <  1:1206325/50‑1 (MQ=255)
                     gttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt                             >  1:1096566‑M2/11‑47 (MQ=255)
                     gttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt                             >  1:3669654‑M2/11‑47 (MQ=255)
                     gttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt                             >  2:4205453‑M2/11‑47 (MQ=255)
                     aaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga                           <  1:2007120/50‑1 (MQ=255)
                     aaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga                           <  2:5294350/50‑1 (MQ=255)
                     aaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGga                             >  2:3732481/1‑48 (MQ=255)
                     aaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGga                             >  2:901521/1‑48 (MQ=255)
                       aCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                         <  1:5738812/50‑1 (MQ=255)
                        cTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg                        <  1:3779589/50‑1 (MQ=255)
                            aCCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaa                    <  2:4236907/50‑1 (MQ=255)
                            aCCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaa                    <  2:4148573/50‑1 (MQ=255)
                            aCCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaa                     >  1:95865/1‑49 (MQ=255)
                            aCCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaa                     >  2:392298/1‑49 (MQ=255)
                            aCCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaa                     >  1:1597373/1‑49 (MQ=255)
                              cGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                  >  2:2626116/1‑50 (MQ=255)
                              cGGACGGTCCCCTCGCCCCTTTAGGGAGAGGGTTAGGGAGAGGGGAAAAc                  <  2:1840768/50‑1 (MQ=255)
                                   ggTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGCGAGGGGAAAACCCGCa             >  2:2115430/1‑50 (MQ=255)
                                    gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc            >  1:2999008/1‑50 (MQ=255)
                                     tCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAcc           >  1:696314/1‑50 (MQ=255)
                                      ccccTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAccc          >  1:1058057/1‑50 (MQ=255)
                                       cccTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAccc          >  2:1747645/1‑49 (MQ=255)
                                       cccTCCCCCTTGTGGGGAGGGGGTGGGGGGGGGGGGAAAACCCGCACCCt         <  1:291705/50‑1 (MQ=12)
                                         cccGCCCCTTTGGGGGGTGGGTAGGGGAGAGGGCAAAACCCGCACCCTAt       <  2:747527/48‑1 (MQ=25)
                                         cTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt       >  1:3997313/1‑50 (MQ=255)
                                            gCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt    >  1:2641604/1‑50 (MQ=255)
                                            cccccTTTGGGGGAGGGGGGGGGGGAGGGGGGAAAACCCGCACCCTAttt     <  2:3468384/49‑1 (MQ=255)
                                             ccccTTTGGGGGAGGGGGGGGGGGAGGGGGTAAAACCCGCACCCTAtttt    <  2:2901636/50‑1 (MQ=255)
                                             ccccTTTGAGGAGAGGGTAAGGGAGGGGGGAAAACCCGCGCCCTATTTTc   <  2:1411024/50‑1 (MQ=25)
                                               ccTTTGGGGGGGGGGGTGAGGGGGGGGGGAAAACCCGCACCCTATTTTCa  <  1:1230818/50‑1 (MQ=255)
                                               ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt    >  2:2586157/1‑47 (MQ=255)
                                               ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt    >  1:2196314/1‑47 (MQ=255)
                                               ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt    >  2:4093021/1‑47 (MQ=255)
                                                cTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt    >  1:4153632/1‑46 (MQ=255)
                                                cTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt    >  1:5238214/1‑46 (MQ=255)
                                                 tttGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt    >  2:5824624/1‑45 (MQ=255)
                                                 tttGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt    >  2:4667149/1‑45 (MQ=255)
                                                 |                                               
TGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCA  >  NC_000913/4070413‑4070508

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