Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,553,864	A→T	13.5%	L214F (TTA→TTT)	hemF →	coproporphyrinogen III oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,553,864	0	A	T	13.5%	203.4 / 4.1	71	L214F (TTA→TTT)	hemF	coproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base A (12/47);  new base T (8/2);  total (20/51)
Fisher's exact test for biased strand distribution p-value = 4.60e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTATTA  >  NC_000913/2553627‑2554102                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     tGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTa                                                                                                                                                                                                                                                      >  2:1028287/1‑238 (MQ=255)                tGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCg                                                                                                                                                                                                                                       <  2:900414/237‑1 (MQ=255)                        acgggCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCg                                                                                                                                                                                                                             <  1:458087/236‑1 (MQ=255)                        aCTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCg                                                                                                                                                                                                                             <  1:248200/239‑1 (MQ=255)                          tGGCATCTCGCCGCCCGAGACCTGTGCCTGCCATTGGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAt                                                                                                                                                                                                                           <  1:836253/239‑1 (MQ=255)                          tGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCCGATTTCGACCGCCGTTTTGCCTTTATGCAGGCGGGA‑GGCAAAGGGCACCACCACGCCTATTTTACAAATGT‑ccgggcacgggaaagcat                                                                                                                                                                                                                           >  1:243140/1‑222 (MQ=255)                          tGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAt                                                                                                                                                                                                                           <  1:636993/239‑1 (MQ=255)                               tCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCt                                                                                                                                                                                                                     <  2:1047778/239‑1 (MQ=255)                                  cACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTTATGACCTGAACACCCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTT‑GGCAAAGGCCACACCGCCGCCTCTTTACCAATTGT‑CCGGCGACGGCAAGCGAT‑gggcctcg                                                                                                                                                                                                                  >  1:22584/1‑233 (MQ=255)                                              cttgcttccgtcATTTGGGAGAGACGTTTATTCCCGTTATGAAAAGTGAGGAGTCTCATACTACTGCCTTCACCATGGCAGCAAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGGCCTGAACACGACAGGTTTCGACCGCTGTCTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGGCGGAAAGCGAT‑GGCCTACGgctagctc‑gagc                                                                                                                                                                                                     <  1:874730/228‑1 (MQ=255)                                               cTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGcgagcgc‑gagcg                                                                                                                                                                                                    <  1:960261/239‑1 (MQ=255)                                                tgtgCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCgagcgc‑gagcgc                                                                                                                                                                                                   <  2:798773/239‑1 (MQ=255)                                                  tgcctgccATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGGTGGTCTTTGATGAACTGAACACGCCAGATTTCGGCCGCTGGTTTGCCTTTATGCAGGAGGGG‑GGCCAAAGCAAAACCGACGCTTATTTTCCCAATTT‑TTTGCGACGGGAAGCCGT‑GGGCTACCgccggggggcaccgca                                                                                                                                                                                                  >  1:704263/1‑225 (MQ=255)                                                  tgcctgccATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTTTGCAGGCGGGT‑GGGAAAGGCGAAACCCACCCCTATTTACCCAATTT‑CTAGGGGCGGGAAGCGAT‑GGGCCACgggcggcccggcgccca                                                                                                                                                                                                  >  1:222395/1‑224 (MQ=255)                                                  tgcctgccATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTTTGCAGGCGGGA‑GGCAAAGGCTAAACCCACGCTTATTTACCCAATTT‑TGAGCGACGGAAAGCGGA‑GGCCTCCGCGGAGGGC‑CGGCGC‑a                                                                                                                                                                                                 >  1:798773/1‑238 (MQ=255)                                                     ctgccATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AAtt                                                                                                                                                                                              <  2:163892/238‑1 (MQ=255)                                                         cATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTTTGCAGGCGGGA‑GGCAAAAGGCACAACCACCGCTCTTTACCAATTTTCCGGGGGGCGGAAA‑CGAT‑GGGCCAAGGCGAGGGC‑GAGCGC‑CAATTCCAg                                                                                                                                                                                         >  1:539780/1‑239 (MQ=255)                                                            ttGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGGTGACCCGAACACCCCAGATTTCCACCGCTCGTTTTCCTTTATGCAGGCGGTAGGGCAAAGGC‑ACCACGCAGCTTATTTTACAATTTT‑CGCAGGGCGGAAAGCGGT‑GGGACACGGCGCGGGC‑GCGAGC‑AAATTCCacgta                                                                                                                                                                                      >  1:961477/1‑235 (MQ=255)                                                                   aGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATAtctt                                                                                                                                                                                <  1:816517/238‑3 (MQ=255)                                                                   aGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATAtcgt                                                                                                                                                                                <  2:320763/238‑1 (MQ=255)                                                                      cGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATAtcgt‑cg                                                                                                                                                                             <  2:1052399/237‑1 (MQ=255)                                                                      cGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGACTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑cgc                                                                                                                                                                            <  2:238325/238‑1 (MQ=255)                                                                            tCCCCGGTACAAAAAGTGGTGCGACGAATACTTCTACCGCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGtt                                                                                                                                                                     <  1:409116/239‑1 (MQ=255)                                                                                 gTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTTATGACCTTAACACGCCAGATTTCGGCCGCTGTTTTGCCTTTATTCCGGGGGGA‑GGGAAAAGGTTCCCCCGCGCCTTTTTTCCAAATTT‑CGAGGGGCCGGAAGAGGTGGGCCTACGG‑GGGCGCGGGGCGA‑AATTTCCCGTTA‑ATTGT‑CGCGGGCGGTTTGTc                                                                                                                                                                >  1:472435/1‑239 (MQ=255)                                                                                      aaaaaGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGGA‑GGGCAAGGCTACACCGACGCCTTTTTTCCAATTTT‑CGAGGGACGGGAAACGGT‑GGCCTACGGCGAGCGC‑GAGCGC‑CATTTTCCATTTTATCTTCCGCGGGCGTTATGtggggt                                                                                                                                                            >  1:64915/1‑234 (MQ=255)                                                                                         aaGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCaa                                                                                                                                                        <  1:931340/239‑1 (MQ=255)                                                                                          aGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAGCAGCGCGGGATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGGCCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAAt                                                                                                                                                       <  2:166582/239‑1 (MQ=255)                                                                                          aGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAAt                                                                                                                                                       <  1:29312/239‑1 (MQ=255)                                                                                          aGTGGTGCGACGAATACTTCTACCGCAACCAGCGCGACGAACAGCTCGGGATTGGCGGGCTGTTCTTTGAGGACCTGAACACGCCAGATTTCGACCGCTGTTTGGCCTTTATGCTGGCGTTA‑GGCAAAGGCTACACCGACGCTTATTTCCCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTGCGGGGAGCGC‑GGGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAAt                                                                                                                                                       <  1:1013749/239‑1 (MQ=255)                                                                                             ggtgCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAAtc                                                                                                                                                      <  2:528139/237‑1 (MQ=255)                                                                                             ggtgCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAAtc                                                                                                                                                      <  2:1185099/237‑1 (MQ=255)                                                                                               tgCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTTATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAAGCTACACCGACGCCTATTTTCCCATTTT‑CCAGCGACGGGAAACGAT‑GGGCTCAGGCGGGCGC‑GAGCGCAAATTTCCAG‑TATATCGT‑CGCGGGCGTTATTTCGGGTTCCAtttgg                                                                                                                                                   >  1:724861/1‑238 (MQ=255)                                                                                                gCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAAtctggtc                                                                                                                                                 <  2:472435/239‑1 (MQ=255)                                                                                                  gacgaATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAAtctggtct                                                                                                                                                <  2:340262/238‑1 (MQ=255)                                                                                                       aTACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATc                                                                                                                                          <  2:724861/239‑1 (MQ=255)                                                                                                                        ttctgaaCGACCAGCGCGGTATTGTGGGTCTGTTCTTTTATGTGCTGAGCACGCCAGATTTCGGCCGCTGTTTCGCCTTTATGCAGGCGGTG‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGG‑CGAGCGACGGAAAGCGAT‑GGCCTACGGTGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGG‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGc                                                                                                                        <  1:471396/234‑1 (MQ=255)                                                                                                                          cGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCt                                                                                                                      <  2:222395/239‑1 (MQ=255)                                                                                                                               cGAACAGCGCGGTATTGGCTGGCTGTTCTTTGATGACCTGAACACGCCGGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACTCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGGT‑GTTTGGCCTGCaaa                                                                                                                 <  2:206171/239‑1 (MQ=255)                                                                                                                               cGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCaaa                                                                                                                 <  2:22584/239‑1 (MQ=255)                                                                                                                                       gcgGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTgg                                                                                                             <  1:1028287/235‑1 (MQ=255)                                                                                                                                             ttGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCAcc                                                                                                   <  2:1084559/239‑1 (MQ=255)                                                                                                                                             ttGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGG‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCAcc                                                                                                   <  2:209101/239‑1 (MQ=255)                                                                                                                                              tGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGGA‑GGCAAAGGCTACACCCACGCTTTTTTTCCCATTTT‑CGAGCGACGGGAAGCGAT‑GGCCCACGGCGAGCGC‑GAGCGC‑CATTTCCAGTTATATTGT‑CGCGGTCGGTATTTCGAGGTTAATCTGTGTTGGGGTTGCGGCCCGCG‑GGTTGGCCTGCACAATggcggggggccc                                                                                                   >  1:977397/1‑232 (MQ=255)                                                                                                                                                  gggCTGTTCTTTGATGACCTGAACACGCCAGATTTCGTCCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGt                                                                                               <  2:624570/238‑1 (MQ=255)                                                                                                                                                  gggCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGt                                                                                               <  2:366903/238‑1 (MQ=255)                                                                                                                                                     cTGTTCTTTGATGACCTGAACACGCCCGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGGCAAAGCTACACCGACGCTTATTTTCCAATTGT‑CGGGCGACGGAAAGCGAT‑GGGCTACGGCCAGCGC‑GAGCGC‑GAATTCCAGTTATATCCC‑TGCGGTCGGTGTGTCTAGTTTAATCTGGGCTGGGATTGGCGACCGCTGGTTGGGCCGCCAACCTGACGGGGGCACCggtctt                                                                                             >  1:230974/1‑233 (MQ=255)                                                                                                                                                         tCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCt                                                                                        <  2:539780/238‑1 (MQ=255)                                                                                                                                                             tgatgaCCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGt                                                                                               <  2:290145/227‑1 (MQ=255)                                                                                                                                                                  aCCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGcgagcgc‑gagcg                                                                                                                                                                                                    >  2:688550/1‑124 (MQ=255)                                                                                                                                                                  aCCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGcgagcgc‑gagcg                                                                                                                                                                                                    <  1:688550/124‑1 (MQ=255)                                                                                                                                                                  aCCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCAc                                                                                                    >  1:535596/1‑217 (MQ=255)                                                                                                                                                                  aCCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCAc                                                                                                    <  2:535596/217‑1 (MQ=255)                                                                                                                                                                      gAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATgccgc                                                                          <  2:64915/239‑1 (MQ=255)                                                                                                                                                                       aaCACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGGCTATCCTGATGTCAATgccgc                                                                          >  2:667957/1‑238 (MQ=255)                                                                                                                                                                        acacGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCa                                                                        <  2:961477/239‑1 (MQ=255)                                                                                                                                                                              cAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGCGTCTATCCTGCTGGCTAAGCCGCCCCCGGt                                                                   >  1:954486/1‑238 (MQ=255)                                                                                                                                                                              cAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTa                                                                  <  2:132680/239‑1 (MQ=255)                                                                                                                                                                              cAGATTGCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATACCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCCCGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTa                                                                  <  2:427290/239‑1 (MQ=255)                                                                                                                                                                                  tttCGACCTCTGTTTTGCCTTTATGAAGGGGGTA‑GGCAAGGGCTACCCCGAGGCTTATTTACCATTTGT‑CTAGCGAGGGAAAGCGAG‑GGCCTACGGCGACGGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCATTCTGGTCTGGGATCGCGGCGCGCT‑GTTTGGCCGGCAAACTGGCGGGGGCAACGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCt                                                              <  1:813642/239‑1 (MQ=255)                                                                                                                                                                                                 tGCTTTTATGCAGTCGGTA‑GGCAAAGGTTACAACAACGCTTTTTTTCCGATTGT‑CGAGGGACGGAAAGCGAT‑GTGCTACGGTGAGCGC‑GAGCTC‑AATTACCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGAGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGCCAATGCCGCCACTGGTACGCTGGGAATATGATTATc                                               <  1:841541/239‑1 (MQ=255)                                                                                                                                                                                                             ggCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGTGTCTATCCTGCTGTCACTGTCCCCCCTGGGACGCTGGGCAAATGATTATTAGCCAAAAGATg                                   >  1:885743/1‑239 (MQ=255)                                                                                                                                                                                                               cGGTA‑GGCAAAGGCTACACCGACGTTTATTTGCGAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGAGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGATTATCAGCCAAAAGATgg                                  <  1:1088480/238‑1 (MQ=255)                                                                                                                                                                                                               cGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGATTATCAGCCAAAAGATgg                                  <  2:977397/238‑1 (MQ=255)                                                                                                                                                                                                                         aggcgtccacCGCCGGTTATTTTCCATTTGT‑CGCGCGACGGAAAGCTAT‑GGCGTACGGCGAGTTC‑GAGCGC‑AATTTGCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCTCGGCACGCT‑GTCTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGTAATATGATTATCAGCCAAAAGATGGCAGCCCAgca                        <  1:892471/232‑3 (MQ=255)                                                                                                                                                                                                                                cacCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGATTAt                                                >  2:682221/1‑208 (MQ=255)                                                                                                                                                                                                                                cacCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGATTAt                                                <  1:682221/208‑1 (MQ=255)                                                                                                                                                                                                                                 acCGACGCTTATTTACCAATTGG‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGtt                <  2:704263/239‑1 (MQ=255)                                                                                                                                                                                                                                  ccGACGCCTATTTACCAATTGT‑CGAGCGACGAAAAGCGAT‑GGCCTACGGCGAGCGC‑CAGCGC‑AATTTCCAGTTATATCGT‑CGCGTTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGAAACCTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTa               <  1:1079886/239‑1 (MQ=255)                                                                                                                                                                                                                                   tgACGCTTATTTACCAATTGT‑CGAGCGAGGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCATTTATATCTG‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTaa              <  2:165998/238‑1 (MQ=255)                                                                                                                                                                                                                                     aCGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTAAGCCGGGAATATGATTATCAGCCAAAAGGTGGCAGCCCAGAAACGGCGTTAAgg            >  2:892288/1‑238 (MQ=255)                                                                                                                                                                                                                                      cGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGGCGCCACTGGTACGCTGGGAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTg           >  1:664797/1‑239 (MQ=255)                                                                                                                                                                                                                                         ttATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAg         <  2:954486/238‑1 (MQ=255)                                                                                                                                                                                                                                               aCCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTttatta  <  2:243140/239‑1 (MQ=255)                                                                                                                                                                                                                                               aCCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTttatta  <  2:885032/239‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     TGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTGCCATTTGGCGAAGACGTTTATCCCCGTTACAAAAAGTGGTGCGACGAATACTTCTACCTCAAACATCGCAACGAACAGCGCGGTATTGGCGGGCTGTTCTTTGATGACCTGAACACGCCAGATTTCGACCGCTGTTTTGCCTTTATGCAGGCGGTA‑GGCAAAGGCTACACCGACGCTTATTTACCAATTGT‑CGAGCGACGGAAAGCGAT‑GGCCTACGGCGAGCGC‑GAGCGC‑AATTTCCAGTTATATCGT‑CGCGGTCGTTATGTCGAGTTCAATCTGGTCTGGGATCGCGGCACGCT‑GTTTGGCCTGCAAACTGGCGGGCGCACCGAGTCTATCCTGATGTCAATGCCGCCACTGGTACGCTGGGAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTATTA  >  NC_000913/2553627‑2554102

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