Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	498,917	T→C	F288S (TTC→TCC)	hemH →	ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	498,917	0	T	C	98.4%	216.0 / NA	62	F288S (TTC→TCC)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (17/44);  minor base G (0/1);  total (17/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGG‑T‑TTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGG  >  NC_000913/498671‑499158                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     attccccGCAACGTTGCCGCACAACGACTCGTGAACGGGCTGCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTc                                                                                                                                                                                                                                                  <  2:494636/247‑1 (MQ=255)    ccccGCAACGTGGCCGCTCAACGACGCGTGAATTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTGTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTTCCCGGGCTTTGCTGCGGACTGTCTGGAGACGCTGGAAGAGATTGCCGTGCAAAACCGTGAGGTCTGCCTCGGt                                                                                                                                                                                                                                               <  2:459540/250‑1 (MQ=255)       cGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGAGTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGcc                                                                                                                                                                                                                                            <  1:135964/251‑1 (MQ=255)        gCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCcg                                                                                                                                                                                                                                           >  1:318153/1‑251 (MQ=255)         cAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCcgg                                                                                                                                                                                                                                          <  1:337215/251‑1 (MQ=255)           aCGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCcggcg                                                                                                                                                                                                                                        <  2:42234/251‑1 (MQ=255)               tGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGa                                                                                                                                                                                                                                     <  2:171918/250‑1 (MQ=255)                gCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGaa                                                                                                                                                                                                                                    <  2:498393/250‑1 (MQ=255)                          gACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAat                                                                                                                                                                                                                         <  1:262206/251‑1 (MQ=255)                                 gAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTc                                                                                                                                                                                                                   <  2:350418/250‑1 (MQ=255)                                 gAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGGCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTcc                                                                                                                                                                                                                  <  2:280846/251‑1 (MQ=255)                                  aaCTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCg                                                                                                                                                                                                                 <  1:371824/251‑1 (MQ=255)                                      ggCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgcg                                                                                                                                                                                                              <  2:412922/250‑1 (MQ=255)                                      ggCTTCCGCATGGGGGATGGCACCGGAAAAAGGGATGATGACCTTTCAGTCGCGCTGTGGGCGTGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCTTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCgtc                                                                                                                                                                                                               <  2:37565/249‑3 (MQ=255)                                          tCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCTCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAAGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTa                                                                                                                                                                                                          <  1:275131/250‑1 (MQ=255)                                          tCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCtt                                                                                                                                                                                                           <  1:26738/249‑1 (MQ=255)                                           ccGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAt                                                                                                                                                                                                        >  1:415343/1‑251 (MQ=255)                                             gCATTGGGGATGGCACCGGAAAAAGTGATGATGATCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGc                                                                                                                                                                                                      <  2:389755/251‑1 (MQ=255)                                               aTTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGaaa                                                                                                                                                                                                                                   <  2:394660/220‑1 (MQ=255)                                                      aTGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGGCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGa                                                                                                                                                                                             <  1:341587/251‑1 (MQ=255)                                                          cACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTGTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGCCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACat                                                                                                                                                                                         <  1:280143/251‑1 (MQ=255)                                                             cGGAAAAAGTGATGATGACCTGTCAGTCGCGCGTTGGTCGGGAACCCTGGCTGAGGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGGGCCGGCGGGAAAAAATATGAATCTATTCCGGCGCTTAATGCCACGCCGGAACATATc                                                                                                                                                                                      <  2:72728/251‑1 (MQ=255)                                                              ggcaaaaGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATAAGAATATATTTCGGCGCTTTATGCCCCCCCGGAACAAATCg                                                                                                                                                                                     >  1:484140/4‑251 (MQ=255)                                                               gAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACataatg                                                                                                                                                                                     >  1:384340/1‑247 (MQ=255)                                                                     gTGATGATGCCCGTGCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCGTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAatgat                                                                                                                                                                              <  2:345224/251‑1 (MQ=255)                                                                          gatgaCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTa                                                                                                                                                                         <  1:352781/251‑1 (MQ=255)                                                                              aCCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATACCACGCCGGAACATATCGAAATGATGGCTAATC‑t                                                                                                                                                                    >  1:284434/1‑251 (MQ=255)                                                                               ccTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATAAGAATATATTCCGGCGCTTAATAACACCCCGGAACATATAGAAatggtaggtcatctt                                                                                                                                                                    >  1:284645/1‑239 (MQ=255)                                                                                  ttCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTTCGGCGCTTAATCCCACGCCGGAACAACTCGAAATGCTGGCTAAGC‑ttgtt                                                                                                                                                                >  1:285923/1‑251 (MQ=255)                                                                                    caggcgcgCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑ttgttg                                                                                                                                                               <  2:352697/246‑1 (MQ=255)                                                                                    cAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGGAAAAatat                                                                                                                                                                                                                             >  1:354618/1‑189 (MQ=255)                                                                                    cAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAatat                                                                                                                                                                                                                             <  2:354618/189‑1 (MQ=255)                                                                                     aGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGGACATATCGAAATGATGGCTAATC‑TTGTTGcc                                                                                                                                                             >  1:200796/1‑251 (MQ=255)                                                                                        cgcgCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCcgcg                                                                                                                                                          <  1:291365/251‑1 (MQ=255)                                                                                           gcTTTGGTCGGGAACCCTGGCTGATGCCTTANACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTAt                                                                                                                                                       <  2:250289/251‑1 (MQ=255)                                                                                                     ggAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCt                                                                                                                                             <  2:456813/251‑1 (MQ=255)                                                                                                       aaCCCTGGCTGATGCCTTATTCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGATGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCGAATC‑TTGTTGCCGCGTATCGCTAAAGCTGa                                                                                                                                           <  2:74939/251‑1 (MQ=255)                                                                                                              gCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTaa                                                                                                                                    <  1:197559/251‑1 (MQ=255)                                                                                                                      cTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTa                                                                                                                                                                         <  2:387500/207‑1 (MQ=255)                                                                                                                      cTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAAAATGAATATATTTCGGCGCTTAATGCCACGCCGGGACATATCGAAATGATgggta                                                                                                                                                                         >  1:387500/1‑204 (MQ=255)                                                                                                                        tataCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCGCTTTATGCCACGCCGGGACATATCGAAATGATGGCTAATC‑TTTTTTCCGCGTATTGCTAAAACTTAGGGGGAAAGAACTGAgg                                                                                                                          >  1:425120/1‑250 (MQ=255)                                                                                                                             cgacgaAACGCTGAAAATGCTCGGAGAAAAANGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGt                                                                                                                     <  2:226146/251‑1 (MQ=255)                                                                                                                                  aaaCGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAAAATTAATATATTCCGGCGCTTAATTCCACGCCGGAACATATCGGAAAGATGGCAAATC‑TTGGTGCCGCGTATCGCTAAAGCTCAGCGGGAAAGAAATGAGCGCCGTcgcgg                                                                                                                >  1:368579/1‑250 (MQ=255)                                                                                                                                     cGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCTCGTATCGCTAAAGCTGAGCGTTAAAGAACTGATCGCCGTCGCGAaga                                                                                                             <  1:389074/251‑1 (MQ=255)                                                                                                                                               gCTCGGAGAAAAAGGCGTAGGTCATATTCAGGGGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcgg                                                                                                   <  2:425120/251‑1 (MQ=255)                                                                                                                                                     agaAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGGAAAAATATGGATATATTTCGGCGGCTAATGCCCCCCCGGAACATATCGAAATGATGGGCAATCCTTTTTACC‑CGGCTCGCCAAAGCTGAGCGGGAAAGAACTGAGCGCCGTCGCGAAGAGcccccgcggcgggtt                                                                                             >  1:37324/1‑238 (MQ=255)                                                                                                                                                       aaaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGG‑T‑TTtc                                                                                         <  1:355428/251‑1 (MQ=255)                                                                                                                                                         aaaGGCGTAGGTCATAGTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCGCGTCGGCGG‑T‑TTtcat                                                                                       <  2:474921/251‑1 (MQ=255)                                                                                                                                                               gTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTTATGCCACGCCGGGACATATCGAAATGATGGCTAATC‑TTTTTGCCGCGTATCGCTCAAACTGAGCGGTAAAGAACTGAGCGCCGTCGCGGAGAGCCTCGTCGGCGG‑G‑TTTCCTCATccc                                                                                 >  1:449567/1‑250 (MQ=255)                                                                                                                                                                  ggTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGG‑t‑tt                                                                                           <  1:436592/238‑1 (MQ=255)                                                                                                                                                                   gTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGG‑T‑TTTCAGCATCCGTGaa                                                                             <  2:37324/251‑1 (MQ=255)                                                                                                                                                                        atTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGG‑T‑TTTCATCATCCGTGAATAATg                                                                        <  2:174511/251‑1 (MQ=255)                                                                                                                                                                             ggTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAAAATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGGTAATC‑TTGTTGCCCCGGATGGCCAAAAATGAGCGGGAAAGAAATGAGAGCCGCCCCGAAGAGCCCCGGCGGCGG‑TGTTTC‑TCATCACTGaataatagaaa                                                                    >  2:381494/1‑245 (MQ=255)                                                                                                                                                                                      cccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGGACATATCGCAATGATgg                                                                                                                                                                            >  1:12193/1‑140 (MQ=255)                                                                                                                                                                                      cccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATgg                                                                                                                                                                            <  2:12193/140‑1 (MQ=255)                                                                                                                                                                                                       gctcttccgatctgCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGGAAAAATATGAAAATATTCCGGCGCTTAATGCCACGCCCGAAAATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAAATGAGGGCCGTCGCGAAGAGCCTCGGCGGCGGGT‑TTTC‑TTATCCCTGGATAAATCAAAAAAACGTGTAGCGGGCCTGGGGagagg                                         >  2:401425/14‑250 (MQ=255)                                                                                                                                                                                                         tGTCTGGAGACGTGGGAATAGATTGCCGAGCAAAACCGTGAGGGCTCCCTCGGGGCCTGCGGGAAAAAATATGAATATATTCCGGCGCTTAAGGCCACGCCGGAACATATCGAAATGATGGCTGATC‑TTGTTGCCGCGTATCGCTAAAGCTGATCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGG‑T‑TTTCATCATCCGTGAATAATGCAAAAAGGCGTGCTGCGTTCCTGGGTAGAGtt                                        <  2:291541/250‑1 (MQ=255)                                                                                                                                                                                                                                   cGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGG‑T‑TTTCATCATcc                                                                                  <  1:169844/182‑1 (MQ=255)                                                                                                                                                                                                                                       cAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGG‑T‑TTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTaa         <  1:398190/251‑1 (MQ=255)                                                                                                                                                                                                                                         aaaCCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGG‑T‑TTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACg       <  1:401425/251‑1 (MQ=255)                                                                                                                                                                                                                                             cGTGAGGTCGCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGG‑T‑TTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTg   <  1:246867/251‑1 (MQ=255)                                                                                                                                                                                                                                              gtggggTCTCCCGCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGG‑T‑TTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTgg  <  2:488743/247‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     ATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATC‑TTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGG‑T‑TTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGG  >  NC_000913/498671‑499158

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