Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,361,5350AT27.9% 62.8 / 10.1 44I226N (ATC→AAC) rhmRputative DNA‑binding transcriptional regulator for the rhm operon
Reads supporting (aligned to +/- strand):  ref base A (20/11);  new base T (0/12);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 9.64e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCCAGCCCTGA  >  NC_000913/2361394‑2361657
                                                                                                                                             |                                                                                                                            
tttAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAAGAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAGTCCCGGTATTCTTCCTTGAt                                                                                                                             >  1:460901/1‑143 (MQ=255)
tttAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACTGTATTCTTCCTTGag                                                                                                                             >  1:591535/1‑142 (MQ=255)
  tgatgtttGGTAGGGTCATGATGCTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGAttt                                                                                                                           <  2:423746/138‑1 (MQ=255)
  tAATTTTTGGTAGGGTCAGGATGTTCTCCATTGTTACGAGGCGTCTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGAttt                                                                                                                           <  1:161315/143‑1 (MQ=255)
      ttttGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCTTAACAGATGTGAAATATCCCTGGAGCAGGCAATGGACTTGACCGCCAGATAATCACGGTATTATTCATTGATTTGTaa                                                                                                                       >  2:699009/1‑143 (MQ=255)
             aGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTa                                                                                                                <  1:149996/143‑1 (MQ=255)
                 tCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCGAGACCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAAGCACGGTAATCTGCGTTGATGTGTAGGCGGGGAC‑‑CGa                                                                                                          >  2:297588/1‑143 (MQ=255)
                  cATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGac                                                                                                         <  1:292577/143‑1 (MQ=255)
                     gatgTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTGC‑‑CGaccac                                                                                                      >  1:350489/1‑143 (MQ=255)
                            tcCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAACCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATc                                                                                               >  1:844703/1‑143 (MQ=255)
                             ccATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCACTGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCg                                                                                              >  1:201762/1‑143 (MQ=255)
                                  gTTATGAGGCTTGTAATTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAAGCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCg                                                                                         <  1:263357/143‑1 (MQ=255)
                                      tGACGCTTGTAAGTCCGCGCGACTTTTCCATCGCACCCGTCGTGAACTGTCCATGGCGCAGGTTATGGCCTTGCTCGCCAGATAATCACCGGATTCTTCTTTTTTTTTTAATGATACGG‑‑CGACCACCGAGATCGCAGCGGTaa                                                                                     <  2:271624/143‑1 (MQ=255)
                                       gcgggttGGAAGGCAAAGGGACTGGTCCGTCCCATCGGACGTGAAATATCGCTGGCGGTTGTGTGGTCGTTGCTCGCCAGATAGCCACCGTAGTCTTCTTTGTTTTTTAA‑‑TGATACGGCGACCACCGAGATCGCAGCGGTAAg                                                                                    <  2:643052/138‑1 (MQ=255)
                                          gCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGACAATGGCCTTGCCCGCCAGCAAATCACGGCATTCTTCATTGCGTTGTAGACGGGTAC‑‑AGAGCACCGAGAGAGCAACgggaggcac                                                                                 >  2:595049/1‑136 (MQ=255)
                                           cTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGCAAAGCTCg                                                                                >  1:18538/1‑143 (MQ=255)
                                            ttGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGAGCAGGCAATGGCCATGGCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCGCAGAGCTAGGAGCGGTAAGCTCgg                                                                               >  1:303142/1‑142 (MQ=255)
                                                     acggtgtcTTTGCTATCCCGTCAGAAGTGAATGATCCATGGCGCAGGTAATGGCCTTGGCCGCCAGATAATCACTGTGTGCTTCTTTGTTTTTTAATGAGGCGG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGt                                                                      <  2:671871/136‑1 (MQ=255)
                                                         gtgctcttCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCATGCGG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTaaa                                                                  <  2:845696/137‑1 (MQ=255)
                                                                ttcgttccAGCAGACGTGAAATTTCCCTGGCGCAGGCCTTGTCTTTGCCCGCCAGATAATCACGCTATTCTTCTTTTTTTTTTAATGATACGG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGgc                                                           <  2:710602/137‑1 (MQ=255)
                                                                tacATCCCGACATACGGTTAATTTCCCTGGCGCAGGCTGTTGCGTTGCGCGCCAGATAATCTCTGTATTATTTTTTTTTTTTTAA‑‑GGGTACGGCGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGgc                                                           <  2:56291/141‑1 (MQ=255)
                                                                  cATCCCAACAGACGTGAAATATCCCTAGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTCAGCGGGTAC‑‑CGACCACCGAGATCGCACCGGTAACCTCGATATAAGCGTTAAAAACCAgcgc                                                         >  2:653005/1‑143 (MQ=255)
                                                                   accccAACAGACGTGAAATCTCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGCTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCTAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCa                                                        <  1:430153/141‑1 (MQ=255)
                                                                   aTCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGTGATCGCAGCGGTACGCTCGTTATTGGCGTTATacactcgatat                                                        >  1:115939/1‑136 (MQ=255)
                                                                       cAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGGTATTGGCGTTAAACACCGGCGTAGCGa                                                    >  1:276428/1‑143 (MQ=255)
                                                                            gACCTGAACTCTCCCGCGCGCAGGCGCTGCCCTCGTCCGCCAGATAGTCACGGTCTTCTTCTTTTTTTTGTAAGCAGACAG‑‑CGACCACCGAGAGCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACaa                                               <  2:644308/143‑1 (MQ=255)
                                                                              cGTTAAATATCCCTGGCGCAGGCACTGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTTTTTTGTAATGGTGCGG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACg                                             <  2:682434/143‑1 (MQ=255)
                                                                               gttatctatCCCTCGTTCAGGCATTGGCCTTGCCCGCCAGATAATCACGGTCTTCTTCTTTGTTTTTTAATCATGCGG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACgg                                            <  2:400414/137‑1 (MQ=255)
                                                                               gTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGCACGCAGCGGAAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACgg                                            >  1:752425/1‑143 (MQ=255)
                                                                                   aaTATCCCTGGCGCAGGCAATGACCTCGCGCACCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCGCGTTATCGGCGTTAAACACCGGGGCAGCGACACAACGGACAt                                        >  2:313719/1‑143 (MQ=255)
                                                                                     tatCTCTGGCGCAGTGATTGGATTTGCCCGCTAGATAGTCGCGGTATGCTTCATTGTTTTGTAAGGGGACGG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCg                                      <  2:702480/143‑1 (MQ=255)
                                                                                          cgggCGCTGGCAATGGCCTTGCTCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTa                                 <  1:382227/141‑1 (MQ=255)
                                                                                          cTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTa                                 >  1:662526/1‑143 (MQ=255)
                                                                                             gcgcAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATc                              >  1:597208/1‑143 (MQ=255)
                                                                                                cagcgaaTGGCCGCGCCCGCCAGAGAAGCACGGGAGTCCGCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTc                           <  1:22903/138‑1 (MQ=255)
                                                                                                          cttggCCGGCCAGTTCATCACGGTATTCTTTGTTGTTTTGTAAGGAGGCAG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCa                 <  2:606477/139‑1 (MQ=255)
                                                                                                           cTTGCCCGCCAGATGTTCACGGTATTCTTCTTTGTTTTGTAATGGGGCGC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCAt                <  2:750605/143‑1 (MQ=255)
                                                                                                           cTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCAt                <  1:179798/143‑1 (MQ=255)
                                                                                                            ttGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCTGCGTAATCTTCGCCGTTGTcgtc               >  1:210420/1‑140 (MQ=255)
                                                                                                              gCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACGCAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGc             >  1:546655/1‑143 (MQ=255)
                                                                                                              gCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACCACGGACATCGGCGTAATCTTCGCCGTTGTCATAGc             >  1:861193/1‑143 (MQ=255)
                                                                                                                 cGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGCCACCGCCGTACTCTTCGCCGTTGTCCTAGCTcc          >  1:754257/1‑143 (MQ=255)
                                                                                                                 cGCCACACAAGCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTcc          >  1:205408/1‑143 (MQ=255)
                                                                                                                   cgaGATAATCACCGTATTCTTCTTTGATTTTTAAGCGTGCGG‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCCAg        <  2:380417/141‑1 (MQ=255)
                                                                                                                         aaTCACGGTATTCTTCTTTGTTTTTTAATGGTGCGC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCCAGCCCTGa  <  2:72015/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                            
TTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTAC‑‑CGACCACCGAGATCGCAGCGGTAAGCTCGTTATTGGCGTTAAACACCGGCGCAGCGACACAACGGACATCGGCGTAATCTTCGCCGTTGTCATAGCTCCAGCCCTGA  >  NC_000913/2361394‑2361657

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: ‑