Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,502,2580AC57.1% ‑1.3 / 22.6 35G56G (GGT→GGG) cysG_1Siroheme synthase
Reads supporting (aligned to +/- strand):  ref base A (8/7);  new base C (15/5);  total (23/12)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGGC  >  KT2440/4502127‑4502388
                                                                                                                                   |                                                                                                                                  
cccGGCAGGCTAAGGGCCGGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGGCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCt                                                                                                                                 <  5:182546/134‑1 (MQ=255)
aCCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCcccc                                                                                                                                  >  4:139184/1‑134 (MQ=255)
 ccGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCCCGCACCAGCCCCGCCCCGCccccgt                                                                                                                                >  3:220910/1‑133 (MQ=255)
     cAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCCGCAACCGCCCAGGTCGGCTGCCGGCTAGCCACGCCCCCGCCCCTCACCGCCCCCTTCCCg                                                                                                                            >  3:285248/1‑135 (MQ=255)
        gCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCCCCCGACCACGTCGGCGGCCAGATAGCCACGCCCCCCCCCCTCCCCCCCCCCATCCCCGGc                                                                                                                         >  7:258086/1‑135 (MQ=255)
        gCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCCCCGTCCCGGGc                                                                                                                         >  3:182402/1‑135 (MQ=255)
                     gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCCACCGCCCAGGTCGGCTGCCTGATAGCCCCGCCCCCGCCCCTCCCCGCCCCCTTCCCGGGCCCCCGCGGccccc                                                                                                            >  1:122586/1‑132 (MQ=255)
                     gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCCCGCCCCCGCCCCTCACCGCCCCCATCCCGGGCCACCGCGGCCacc                                                                                                            >  1:217039/1‑133 (MQ=255)
                     gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCCCGCACCAGCCCCTCACCGCCACCGTCCCGGGCCAACGCGGCCCAc                                                                                                            >  5:171471/1‑135 (MQ=255)
                      cATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCCACCGCCCACGTCGGCTGCCGGCTAGCCCCGCCCCCGCCCCGCCCCGCCCCCATCCCCGGCCCACCCCGCCCAca                                                                                                           >  3:195346/1‑134 (MQ=255)
                      cATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCCCCAGCCCCTCCCCGCCCCCATCCCGGGCCCACGCGGCCAACt                                                                                                           >  2:273672/1‑135 (MQ=255)
                            ccGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGGCTGCTAGCCACGCACCAGCACCTCCCCGCCCCCGTCCCGGGCCACCGCGGCCAACAGGCCgg                                                                                                     >  4:264729/1‑134 (MQ=255)
                                 acCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGCTCACCAGCCGGCCACCGCCCAGGTCGGCGGCCGGCTAGCCCCGCCCCCGCCCCGCCCCGCCCCCGTCCCGGGCCACCGCGGCCCACTGGCCGGCGCCg                                                                                                >  8:208470/1‑135 (MQ=255)
                                        cATTCAGCCCAGGGTCGTCAGTCGCCGCGCTCACCCGCCGGCCCCCGCCCAGGTCGGCGGGCTGATAGCCCCGCCCCCCCACCGCACCGCCCCCGTCCCGGGCCCCCGCCGCCCACTGGCCGGCGCCGGCCGgcg                                                                                         >  5:82562/1‑135 (MQ=255)
                                        cATTCAGCCCAGGGTCGTCAGTCGCCGCGCTCACCAGCCGGCCACCGCCCAGGGCGGCTGCCTGATAGCCACGCCCCCGCACCTCCCCGCCCCCGTCCCGGGCCCACGCGGCCCCCAGGCCGTCCGCGGCCGgcg                                                                                         >  3:181582/1‑135 (MQ=255)
                                         aTTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCCCCTCACCGCCCCCTTCCCGGGCCCACGCGGCCAACTGGCCGTCGACGTCCGgcgc                                                                                        >  7:239549/1‑135 (MQ=255)
                                         aTTCAGCCCAGGGTCGTCAGTAGCCGCGATAACCATCAGGCAACCCACCCGGTCGGCTGCCAGCGAGACACCCACCATCAACTAACCGCACCCATACCCGCACAACGCGGCCCACAGGCCGTCGAAGACCGgcgc                                                                                        >  8:21951/1‑135 (MQ=255)
                                                ccAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCCACCGCCCAGGTCGGCTGCCGGATAGCCCCGCCCCAGCACCTCACCGCCCCCGTCCCCGGCCCACCCGGCCCACAGGCCGGCGCCCGCCCGCGcccccccc                                                                                 >  5:97558/1‑132 (MQ=255)
                                                ccAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCCCCAGCCCCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCAcg                                                                                 >  7:40184/1‑135 (MQ=255)
                                                    ggTCGGCAGGCGCCGCGATCACCCCCCGGCAACCCACCCGGTGGGCTGCCTGATAGCCACGCACCAGCACCCCACCCCCACCTTCCCGGCCCACCGCGCCCAACTGGCCGTCGACGTCCGGCGCCACCACGCgca                                                                             <  4:220910/135‑1 (MQ=255)
                                                            gTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCCCCCCCTTCCCGGGCCAACCCGGCCAACAGGCCGTCGACGTCCGGCGCCACCACCCGCAGCAcgccg                                                                     >  5:187295/1‑135 (MQ=255)
                                                                ccgcGATCACCAGCCGGCAACCGACCAGGGCGGCTGCCTGATAGCCCCGCACCAGCACCTCACCGCCCCCTACCCGGGCCACCGCGGCCAACTGGCCGGCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCgg                                                                 <  8:186516/135‑1 (MQ=255)
                                                                       cACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGc                                                          <  7:66342/135‑1 (MQ=255)
                                                                                 cacccGACCTGGGCGGCTGCCTGATAGCCCCCCACCCCCACCCCACCCCCCCCTCCCCGGCCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGgcgc                                                <  4:285248/132‑1 (MQ=255)
                                                                                 cAACCGACCCGGTCGGCGGCCTGATAGCCCCGCCCCAGCACCTCCCCCCCCCCTTCCCGGGCCAACGCGGCCAACGGGCCGGCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGgcgc                                                <  7:208470/135‑1 (MQ=255)
                                                                                       aCCCGGTGGGCGGCCTGATACCCCCGCACCCGCACCTCACCCCCCCCTCCCCGGCCCACCGCGGCCAACTGGCCGCCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTg                                          <  7:301993/135‑1 (MQ=255)
                                                                                          aGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCCCCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACCTCCGGCGCCCCCCCGCGCCGCCCGCCGCCGGCCTCGGCCAGCAGGCGCGCCGTgcgc                                       >  7:316012/1‑135 (MQ=255)
                                                                                               ggCTGCCAGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCCCCACTCGCAGCCCGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCCAGGc                                  >  1:21643/1‑135 (MQ=255)
                                                                                                         tAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGACGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCCGGCCCGCCTTGCGCCAGGCGATCTCACCg                        >  8:193568/1‑135 (MQ=255)
                                                                                                           gcaccccACCCGCTCCTCCCCGCCACATTCCCGGCCCACCGCGGCCACGTGGCGGTCTACGTCGGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCTGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGcc                      <  7:329487/129‑1 (MQ=255)
                                                                                                            ccACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGcc                      <  8:215402/134‑1 (MQ=255)
                                                                                                                 cACCAGCACCTCACCGCCACCTCCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGACAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACcgac                <  4:73688/135‑1 (MQ=255)
                                                                                                                         ccTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCATGCGCAAGGCGATCTCCCCGCCACCGACGCccacca        >  1:309054/1‑135 (MQ=255)
                                                                                                                         ccTCACCGCCACCTTCCCGCGCCAACGCGGCCACCTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTAGCGCAAGGCGATCTCACCGCCACCCACGACCAGCa        >  6:191321/1‑135 (MQ=255)
                                                                                                                           tCACCGCCCCGGCCCGGGGCCACCGCGGCCAAGTGGCCGTGGAGGTCGGGCGCCACCACGCGCAGCACGCCGCCGGCATGGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACTACCAGCAcg      <  7:260613/135‑1 (MQ=255)
                                                                                                                            ccccGCCACCTCCCCGGGCCACCGCGGCCAACGGGCCGGCGACGTCCGGCGCCCCCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCAcgc     <  7:126795/133‑1 (MQ=255)
                                                                                                                               cGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAACGCGATCTCACCGCCCCCGACGCCCAGCACGCGgc  >  1:284105/1‑135 (MQ=255)
                                                                                                                                   |                                                                                                                                  
ACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGGC  >  KT2440/4502127‑4502388

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