Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,037,5861.T68.8% 39.0 / 13.3 16coding (567/1359 nt)hprSsensory histidine kinase HprS
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base T (2/9);  total (3/13)
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.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TAGTTCTCTCGGTAATGCGCTAACCTCAACAGGCTCTCGGCTATCGTTATAATTCAGCGCTTCCGTTACACCACTCAACTTTTTGATCTCTCGTAATCCCGTTCTGATTAACAGCGGACTTAATACTGAGCAAAGTACAATGGCGACAATGCAAAT  >  NC_000913/2037508‑2037663
                                                                               |                                                                             
tAGTTCCCCCGGTAATGCGCTAACCTCAACAGGCTCTCGGCTATCGTTATAATTCAGCGCTTCCGTTACACTACTCAACtt                                                                             >  2:335529/1‑81 (MQ=255)
  gTTCTCTCGGTAATGCGCTAACCTCAACAGGCTCTCGGCTATCGTTATAATTCAGCGCTTCCGTTACACCACTCAACtttt                                                                           <  1:335495/81‑1 (MQ=255)
  gTTCTCTCGGTAATGCGCTAACCTCAACAGGCTCTCGGCTATCGTTATAATTCAGCGCTTCCGTTACACCACTCAACtttt                                                                           <  2:335496/81‑1 (MQ=255)
          ggTAATGCGCTAACCTCAACAGGCTCTCGGCTATCGTTATAATTCAGCGCTTCCGTTACACCACTCAACTTTTTTGAtctc                                                                    <  1:335497/81‑1 (MQ=255)
              aTGCGCTAACCTCAACAGGCTCTCGGCTATCGTTATAATTCAGCGCTTCCGTTACACCACTCAACTTTTTTGATCTCTCg                                                                 <  2:335498/80‑1 (MQ=255)
              aTGCGCTAACCTCAACAGGCTCTCGGCTATCGTTATAATTCAGCGCTTCCGTTACACCACTCAACTTTTTGATCTCTCGt                                                                >  2:335527/1‑80 (MQ=255)
                  gcTAACCTCAACAGGCTCTCGGCTATCGTTATAATTCAGCGCTTCCGTTACACCACTCAACTTTTTTGATCTCTCGTAATc                                                            <  2:335499/81‑1 (MQ=255)
                      aCCTCAACAGGCTCTCGGCTATCGTTATAATTCAGCGCTTCCGTTACACCACTCAACTTTTTTGATCTCTCGTAATCCCGt                                                        <  1:335500/81‑1 (MQ=255)
                                gCTCTCGGCTATCGTTATAATTCAGCGCTTCCGTTACACCACTCAACTTTTTTGATCTCTCGTAATCCCGTTCTGATTAAc                                              <  1:335501/81‑1 (MQ=255)
                                 ctctCGGCTATCGTTATAATTCAGCGCTTCCGTTACACCACTCAACTTTTTTGATCTCTCGTAATCCCGTTCTGATTAACa                                             <  2:335502/81‑1 (MQ=255)
                                              ttATAATTCAGCGCTTCCGTTACACCACTCAACTTTTTTGATCTCTCGTAATCCCGTTCTGATTAACAGCGGACTTAATAc                                <  2:335503/81‑1 (MQ=255)
                                                   aTTCAGCGCTTCCGTTACACCACTCAACTTTTTGATCTCTCGTAATCCCGTTCTGATTAACAGCGGACTTAATACTGAGCa                          <  1:335504/81‑1 (MQ=255)
                                                         cgcTTCCGTTACACCACTCAACTTTTTTGATCTCTCGTAATCCCGTTCTGATTAACAGCGGACTTAATACTGAGCAAAGTa                     <  1:335505/81‑1 (MQ=255)
                                                              ccGTTACACCCCTCCACTTTTTGGTCTCTCCTAATCCCGTTCTGATTAACAGCGGACTTAATACTGAGCAAAGTACAATg                <  1:335506/80‑1 (MQ=255)
                                                                  tACACCACTCAACTTTTTTGATCTCTCGTAATCCCGTTCTGATTAACAGCGGACTTAATACTGAGCAAAGTACAATGGCGa            >  1:335531/1‑81 (MQ=255)
                                                                   acacCACTCAACTTTTTTGATCTCTCGTAATCCCGTTCTGATTAACAGCGGACTTAATACTGAGCAAAGTACAATGGCGa            <  1:335507/80‑1 (MQ=255)
                                                                       cacTCAACTTTTTTGATCTCTCGTAATCCCGTTCTGATTAACAGCGGACTTAATACTGAGCAAAGTACAATGGCGACAATg       >  2:335532/1‑81 (MQ=255)
                                                                            aaCTTTTTGATCTCTCGTAATCCCGTTCTGATTAACAGCGGACTTAATACTGAGCAAAGTACAATGGCGACAATGCAAAt  <  2:335509/80‑1 (MQ=255)
                                                                            aaCTTTTTGATCTCTCGTAATCCCGTTCTGATTAACAGCGGACTTAATACTGAGCAAAGTACAATGGCGACAATGCAAAt  <  1:335508/80‑1 (MQ=255)
                                                                               |                                                                             
TAGTTCTCTCGGTAATGCGCTAACCTCAACAGGCTCTCGGCTATCGTTATAATTCAGCGCTTCCGTTACACCACTCAACTTTTTGATCTCTCGTAATCCCGTTCTGATTAACAGCGGACTTAATACTGAGCAAAGTACAATGGCGACAATGCAAAT  >  NC_000913/2037508‑2037663

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