Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,361,5500CG28.0% 110.7 / 15.6 50G221A (GGT→GCT) rhmRputative DNA‑binding transcriptional regulator for the rhm operon
Reads supporting (aligned to +/- strand):  ref base C (18/18);  new base G (0/14);  total (18/32)
Fisher's exact test for biased strand distribution p-value = 6.81e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCC  >  NC_000913/2361451‑2361649
                                                                                                   |                                                                                                   
ggACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAc                                                                                                     >  1:308076/1‑100 (MQ=255)
ggACTTTTCCATCCCAACAGACGTGAAATATCACTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAc                                                                                                     >  1:1492580/1‑100 (MQ=255)
 gACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAcc                                                                                                    >  1:704680/1‑100 (MQ=255)
  aCTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCg                                                                                                   <  2:1349795/100‑1 (MQ=255)
        ccATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGaccacc                                                                                             <  1:792734/100‑1 (MQ=255)
          aTCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCga                                                                                           <  1:925023/100‑1 (MQ=255)
           tCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCgag                                                                                          <  1:732591/100‑1 (MQ=255)
            cccAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCgaga                                                                                         >  1:920755/1‑100 (MQ=255)
             ccAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCgaga                                                                                         >  2:840815/1‑99 (MQ=255)
             ccAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCgag                                                                                          >  1:623187/1‑98 (MQ=255)
              cAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATc                                                                                       <  1:831056/100‑1 (MQ=255)
                aCAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGc                                                                                     <  1:545046/100‑1 (MQ=255)
                     cGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAGCGACCACCGAGATCGCAGCgg                                                                                <  2:1353924/100‑1 (MQ=255)
                     cGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCgg                                                                                <  1:139654/100‑1 (MQ=255)
                      gTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGt                                                                               >  1:1404526/1‑100 (MQ=255)
                       tGAAATATCCCTGGCGCAGGTAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCTTTGTTTTTTAATGATACGGCGACCACCGAGATCGCAGCGGTa                                                                              <  2:1092870/100‑1 (MQ=255)
                        gAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGCACCGACCACCGAGATCGCAGCGGTaa                                                                             <  2:1556523/100‑1 (MQ=255)
                          aaTATCCCTGTCGCAAGCATTGTCCTCTCTCGCCAGATAACCACGGTATTCTTCTTTTTTTTTTAAGGAGACGGCGACCACCGAGATCGCAGCGGTAAGc                                                                           <  2:1497145/100‑1 (MQ=255)
                              tCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGt                                                                       >  1:1241167/1‑100 (MQ=255)
                               cccTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGtt                                                                      >  2:1196478/1‑100 (MQ=255)
                               cccTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCATCATTGATTTGTAAGCGGGTACCGACCACCGAGAACGCAGCGGTAAGCTCGtt                                                                      >  2:279281/1‑100 (MQ=255)
                                ccTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTa                                                                     >  1:1002290/1‑100 (MQ=255)
                                  tGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCTTTTTTTTTTAATGATACGGCGACCACCGAGATCGCAGCGGTAAGCTCGTTAtt                                                                   <  2:1429415/100‑1 (MQ=255)
                                   ggCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTg                                                                  >  1:292392/1‑100 (MQ=255)
                                    gcgcAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTgg                                                                 <  1:1278270/100‑1 (MQ=255)
                                     cgcAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGc                                                                <  1:417790/100‑1 (MQ=255)
                                        aGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGtt                                                             <  2:920747/100‑1 (MQ=255)
                                        aGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGtt                                                             <  1:1208594/100‑1 (MQ=255)
                                         ggCAATGTCCTTGCCCGCCAGATAATCACGGTATTCTTCTTTGTTTTTTAATGGTACGGCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTa                                                            <  2:1515291/100‑1 (MQ=255)
                                           cAATGGCCCTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTaaa                                                          >  1:977286/1‑100 (MQ=255)
                                            aaTGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTCaa                                                           >  2:1069287/1‑98 (MQ=255)
                                             aTGGTCTTGCCCGCCAGATAATCACGGTATTCTTCATTTATTTGTAAGCGGGTAGCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAaca                                                        <  2:345838/100‑1 (MQ=255)
                                             aTGGCCTTGCCCGCCAGATAATCACGGTATGCTTCTTTGTTTTTTAATGATACGGCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAaca                                                        <  2:946669/100‑1 (MQ=255)
                                                gCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCg                                                     >  2:641364/1‑100 (MQ=255)
                                                  cTTGCCCGCCAGATAATCACGGTATTCTTCTTTTTTTTTTAAGGATGCGGCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGgc                                                   <  2:1261017/100‑1 (MQ=255)
                                                  cTTGCCCGCCAGATAATCACGGTATTCTTCTTTGTTTTTTAATGATACGGCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCgg                                                    <  2:219227/99‑1 (MQ=255)
                                                  cTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGgc                                                   <  1:585880/100‑1 (MQ=255)
                                                       ccGCCAGATAATCACGGTATTCTTCTTTGTTTTGTAATGATGCGGCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGc                                              <  2:1416689/100‑1 (MQ=255)
                                                        cGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCg                                             >  2:884055/1‑100 (MQ=255)
                                                         gCCAGATAATCACGGTATTCTTCATTTTTTTGTAATCGGACAGCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCa                                                <  2:938648/96‑1 (MQ=255)
                                                            aGATAATCACGGTATTCTTCTTTTTTTTGTAATGATGCGGCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGacac                                         <  2:172500/100‑1 (MQ=255)
                                                               tAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAAc                                      >  2:1418483/1‑100 (MQ=255)
                                                                aaTCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCCGCGCAGCGACACAACg                                     >  2:373529/1‑100 (MQ=255)
                                                                     cGGTATTTTTCTTTGTTTTTTAATCGTGTGCCGACCACCGAGACCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACa                                 <  2:1492570/99‑1 (MQ=255)
                                                                        tATTCTTCTTTTATTTGTAATGATACAGCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCgg                             <  2:483471/100‑1 (MQ=255)
                                                                            cttcATTGATTTGTAAGCAGACACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCg                           <  2:704677/98‑1 (MQ=255)
                                                                              tcATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAAt                       <  1:1368202/100‑1 (MQ=255)
                                                                                 ttGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGGATAATGtt                    >  2:359206/1‑100 (MQ=255)
                                                                                    aTTTGTAAGCGGGCAGCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCg                  <  2:575692/99‑1 (MQ=255)
                                                                                     tttGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGcc                <  1:279282/100‑1 (MQ=255)
                                                                                      ttGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCg               >  1:1390141/1‑100 (MQ=255)
                                                                                        gTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTGATACGCGGCAAAAACCGGCGCCGCTGCCCAACCGCCATCCCCGTCAGCGTCGCCgc              >  2:1418667/1‑98 (MQ=255)
                                                                                            gCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCa         <  1:641367/100‑1 (MQ=255)
                                                                                                   ccGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAAAGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTcc  >  1:1102370/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCC  >  NC_000913/2361451‑2361649

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑