Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,361,5330AT20.9% 75.2 / 11.3 47I226N (ATC→AAC) rhmRputative DNA‑binding transcriptional regulator for the rhm operon
Reads supporting (aligned to +/- strand):  ref base A (19/15);  new base T (0/9);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 2.48e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCCAGCCCTGACGGCGAATACGCGCCAGTT  >  NC_000913/2361391‑2361674
                                                                                                                                              |                                                                                                                                               
tgtgcttttttGGTCGGGTCATGAGGTGTTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACGTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGACTTGCCCGCCAGATAATCACGGTATTCTTCATTGa                                                                                                                                                 <  2:448811/137‑1 (MQ=255)
 tttAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGAt                                                                                                                                                >  1:19171/1‑143 (MQ=255)
 tttAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGAt                                                                                                                                                >  2:569037/1‑143 (MQ=255)
  ttAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGAtt                                                                                                                                               <  2:431796/143‑1 (MQ=255)
   tAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGAttt                                                                                                                                              <  2:733032/143‑1 (MQ=255)
         ttGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGc                                                                                                                                        <  2:529270/143‑1 (MQ=255)
          tGGTAGGGTCATGATGTTCTCCATGGTTAGGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCg                                                                                                                                       <  1:357997/143‑1 (MQ=255)
           ggTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCACCAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATCACCACGGGCTGCTGCATTGACTGGTAGGCgg                                                                                                                                      >  2:314056/1‑143 (MQ=255)
           cgTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCgg                                                                                                                                      <  1:323590/142‑1 (MQ=255)
            gTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCGGATAATCCCGGCTTTCTTCCTTCGTTTGTAAACgga                                                                                                                                     >  1:10665/1‑142 (MQ=255)
                    atgatgTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAG‑‑CGacc                                                                                                                           >  1:465996/1‑143 (MQ=255)
                      gatgTTCTCAATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTCTTCTTCATTGATTGGTAAGCGGGTAC‑‑CGaccac                                                                                                                         >  2:570510/1‑143 (MQ=255)
                        tgTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATGTCCCTGGCGCAGGCACTGGCCTTGCCCGCCAGATAATCACGGTACTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCg                                                                                                                       >  2:285809/1‑143 (MQ=255)
                                 ttGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCTATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAg                                                                                                              <  1:483200/143‑1 (MQ=255)
                                   gTTATGGGGGTGGTAAGTGAAAGGGACTTGTCCATCGCAACAGACGTGAAATATTCCTGGCGCAGGCTATTGCCTTGCCCGCCAGATAAGCACGGTATTCTTCTTTTTTTTTTAATGGTGCGG‑‑CGACCACCGAGATCGCAGCg                                                                                                            <  2:306511/143‑1 (MQ=255)
                                           gCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTc                                                                                                    <  1:701672/143‑1 (MQ=255)
                                                  aGTCAATGGGACTTTTCCATCCCATCAGACGTGAATTATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCTTTTTTTTTTAATGATACGG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTg                                                                                             <  2:42360/143‑1 (MQ=255)
                                                   gTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGACCGCAGCGGGAAGCTCGTCATTgg                                                                                            >  2:828815/1‑143 (MQ=255)
                                                     tacaGGGACTTTCCCAGCCGAACAGACCTGAAATATCCCTCGCGCAGGCGATTGTCTGGCCCGTCAGATAAGGACGGTAGTGTCCTTTGTTTTTTAA‑‑TGGTACGGCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCg                                                                                          <  2:777367/140‑1 (MQ=255)
                                                        aGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTa                                                                                       <  1:816187/143‑1 (MQ=255)
                                                         gggACTGTTCCATCCCAACAGACGTTAAATATCCCTGGCGCAGGCAATGGCCTTGCCGGCCAGATTATCATGAGATTCTTCTTTTTTTTGTAATGATACGG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTaa                                                                                      <  2:695150/143‑1 (MQ=255)
                                                                  ccATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCCCGGTCTTCTTCACTGCTCGGCAAGAGGGTTC‑‑CGTCCTCCCAGGTCGCCGCGGTGAGCCCGTTCCTGGCGGCACACACAGgcg                                                                             >  2:834482/1‑143 (MQ=255)
                                                                     tCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAg                                                                          <  1:422423/143‑1 (MQ=255)
                                                                        ccgctgccGTGTAATATCGCTGGCGCAGTCAATGGCTTTGCCGGCCAGACCATCACCGTATTCTTCTTTGTTTTTTAA‑‑GGGTACGGCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGa                                                                       <  2:591904/136‑1 (MQ=255)
                                                                        cAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGa                                                                       >  1:346834/1‑143 (MQ=255)
                                                                          ataggcGTGAAATAGCCGTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGaca                                                                     <  1:777694/138‑1 (MQ=255)
                                                                          aCAGACGTGAAATATCCCGGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTTGTTATTGGCGTTAAACACCGGCGCAGCGaca                                                                     <  1:442619/143‑1 (MQ=255)
                                                                                 tGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTCGCGTTCAACACCGGCGCAGCGGCTCGACggg                                                              >  2:600524/1‑142 (MQ=255)
                                                                                 tGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGCAAGCTCGTTATCGGCGTTAAACCCCGGCGCAGCGCCGCAACGGa                                                              >  2:600531/1‑143 (MQ=255)
                                                                                  gagcgatcgggggCTCAGGCAATGGCCTTGCTCGGCAGATAATCACGGTATTCTTCATTTTTTTTTAAGGGTGCGG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGAc                                                             <  2:730737/132‑1 (MQ=255)
                                                                                  gaattttCCCTGGCGTTGTCTTTGGCCTTGCCCGCCAGATAATCACGTTATTCTTCTTTTTTTTTTAATGAGGCGG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGAc                                                             <  2:212104/137‑1 (MQ=255)
                                                                                  gAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTGATTGATTTGTAAGCGGATAC‑‑CGACCACCGAGATCGCAGCGGGAAGCTCGTTATGGGCGTTCGACTCTGGCGCAGCGCCACAGCGGAc                                                             >  2:357710/1‑143 (MQ=255)
                                                                                  gAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAACCACGGTCTTCTTCATTGCTTCGTACGCGGGTCC‑‑CCACCACCCCGCCCGCAGCAGCACGCCCGCCACCGGCGGTCAACACCCGCGCacccgcgccgctcga                                                             >  2:535360/1‑129 (MQ=255)
                                                                                             ggCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTCAGCCCCGGCGCAGCGACCCCACGGTCATCGGCGtcct                                                  >  1:37435/1‑140 (MQ=255)
                                                                                              gcgcAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAGCCCGGTAGTGGGGGTTAACAACGGCGCAGGGACACAACGGACGATGGGGTAATCt                                                >  1:772680/1‑143 (MQ=255)
                                                                                                   ggCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGGAAGCTCGTTATTGGCGTTAAAAACAGCTACTTCGTCACAACGGACATCGGCGTAATCTTCGc                                            >  1:432875/1‑143 (MQ=255)
                                                                                                       ctgcccTCGCCCGCCAGATGCTCACGCTTTTCTTCTTTTTTTTTTAAGGATACGG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGtt                                        <  2:264586/139‑1 (MQ=255)
                                                                                                        tGGCCTTGTCCGCCAGATAATCACGGTGTTCTTCTTTTTTTTTTAATGATACGG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTg                                       <  2:644922/143‑1 (MQ=255)
                                                                                                          gCCTTGCCCGCCAGATAATCACGGTTTTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGcccctgcc                                     >  1:616821/1‑137 (MQ=255)
                                                                                                                   gCCAGGTAAGCACGGTCTTCTTCATTGCTTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCCa                            <  1:221336/143‑1 (MQ=255)
                                                                                                                       gATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCCAGccc                        <  1:569037/143‑1 (MQ=255)
                                                                                                                               cGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCCAGCCCTGACGGCg                >  1:566556/1‑143 (MQ=255)
                                                                                                                                    ttcttcAATGATTTGTAAGCGGGTAG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCCAGCCCTGACGGCGAATAc           <  1:522932/143‑1 (MQ=255)
                                                                                                                                      cttcATTGATTTGTAAGCGGGTGC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCCAGCCCTGACGGCGAATAcgc         <  2:400020/143‑1 (MQ=255)
                                                                                                                                          aTTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCCAGCCCTGACGGCGAATACGCGCCa     <  1:517446/143‑1 (MQ=255)
                                                                                                                                          aTTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCACGTTATTAGCGTTAACCACCCGCCCCGCGCCACCACCGACATGGGGGTAAGCCTCGCCGTTGTCATAGCTCCAGGCGTTACGGCGCAAAcgccaca     >  2:339341/1‑139 (MQ=255)
                                                                                                                                           ttGATTTGCAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGCAAGACCGTTATTCGCGCTAAACACCCGCCCAGCGACACAACGGACCTCGGCGTAATCCTCGCAGCTCTCATAGCTCCAGGCCTGCAGGCGGATACACGccgg    >  2:705743/1‑141 (MQ=255)
                                                                                                                                             gATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTGGTCATCGCTCCAGCCCTGACGGCGAATACGCGCCAGtt  >  1:451390/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                               
GTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCCAGCCCTGACGGCGAATACGCGCCAGTT  >  NC_000913/2361391‑2361674

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