Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,888,8200TC20.7% 54.1 / 10.4 29Q427R (CAG→CGG) cysIsulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding
Reads supporting (aligned to +/- strand):  ref base T (11/12);  new base C (3/3);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCA  >  NC_000913/2888679‑2888961
                                                                                                                                             |                                                                                                                                             
aTCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGt                                                                                                                                        >  4:23653/1‑149 (MQ=255)
      atatGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCCGCAGCGTGaaaac                                                                                                                                  <  1:30786/149‑5 (MQ=255)
           cTCATCGCTGACACCATGTTTCGCCATTAAATTATCGAAGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcattc                                                                                                                              <  8:41919/148‑1 (MQ=255)
            tCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACACTGCGGCGTGACGGcattcat                                                                                                                            >  5:157252/1‑149 (MQ=255)
                         cATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGctc                                                                                                                  >  1:6051/1‑146 (MQ=255)
                         cATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGctc                                                                                                                  <  2:6051/146‑1 (MQ=255)
                               tCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGc                                                                                                          <  3:23653/148‑1 (MQ=255)
                               tCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGc                                                                                                          <  6:139062/148‑1 (MQ=255)
                               tCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACACTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGc                                                                                                          <  6:157252/148‑1 (MQ=255)
                                  ccATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGAtc                                                                                                      >  8:20139/1‑149 (MQ=255)
                                        aaTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTtctc                                                                                                 <  7:20139/148‑1 (MQ=255)
                                               gATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtc                                                                                          <  2:190665/148‑1 (MQ=255)
                                                                                   ttCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGAtcttct                                                                                                  >  8:58579/1‑104 (MQ=255)
                                                                                   ttCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGAtcttct                                                                                                  <  7:58579/104‑1 (MQ=255)
                                                                                                             ggAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCtt                                                                                                              <  5:25441/66‑1 (MQ=255)
                                                                                                             ggAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCtt                                                                                                              >  6:25441/1‑66 (MQ=255)
                                                                                                                  gACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCtt                       <  6:145663/148‑1 (MQ=255)
                                                                                                                  gACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCtt                       >  5:145663/1‑148 (MQ=255)
                                                                                                                     acGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGC‑TTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTtgtg                   >  2:80670/1‑148 (MQ=255)
                                                                                                                         aaGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATc               >  2:39251/1‑149 (MQ=255)
                                                                                                                               tCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGc                                                                >  4:108776/1‑94 (MQ=255)
                                                                                                                               tCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGc                                                                <  3:108776/94‑1 (MQ=255)
                                                                                                                                   gTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGa                 <  5:56751/137‑1 (MQ=255)
                                                                                                                                   gTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGa                 >  6:56751/1‑137 (MQ=255)
                                                                                                                                     cgtCACGCCGCCGCGGGACGGCATTCATTAACCCGCGCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTc                                             <  8:124970/105‑1 (MQ=255)
                                                                                                                                     cgtCACGCCGCAGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTc                                             >  7:124970/3‑107 (MQ=255)
                                                                                                                                     cgtCACGCCGCAGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGAtctc    >  1:24041/3‑148 (MQ=255)
                                                                                                                                         aCGCCGCAGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCca  <  5:67693/146‑1 (MQ=255)
                                                                                                                                         aCGCCGCAGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCca  >  6:67693/1‑146 (MQ=255)
                                                                                                                                             |                                                                                                                                             
ATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCA  >  NC_000913/2888679‑2888961

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑