Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,297,469	C→T	G633D (GGC→GAC)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,297,469	0	C	T	98.4%	186.2 / NA	63	G633D (GGC→GAC)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (32/29);  minor base G (0/1);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGCCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGT  >  NC_000913/4297223‑4297712                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    cccatggCCGCTGCGCGCAGGCGAGTTATCAACTGGTTGTACTCGGCAATACCGTACGTCTCGGCGGCGCTCTGTTCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAAGCTGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑TTTGATAATTTTGTCttt                                                                                                                                                                                                                                                 <  2:383229/245‑1 (MQ=255)  ccAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTa                                                                                                                                                                                                                                                <  1:3189/251‑1 (MQ=255)         cGCTCCGCGCAGGCGAGGTTTCAACTTGTTGACCTCGTCGATCACGTACTGCTCGGCGGCGCGCCTATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGt                                                                                                                                                                                                                                         <  2:117313/251‑1 (MQ=255)                cgcAGGCGAGTTTTCAACTTGTTGTACTCGTCAGTCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTTCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCa                                                                                                                                                                                                                                  <  1:399898/251‑1 (MQ=255)                  cAGGCGAGTTTTCAACTTGTTGTACTCGTCAGTCACGTACTGCTCGGCGGAGAGCTGATAGGCGGTCGGCTCGCAGCAAAAGGCGCAGTATTTTTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGAAAACAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTTATGTAACTCGCCCCTTTGTTCGG‑ACGATCGCTGCCCTGCGCACG‑GGGGCGACTTTTGTCTTTACGCGAgggacaacaa                                                                                                                                                                                                                                 >  2:568813/1‑241 (MQ=255)                      cGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCa                                                                                                                                                                                                                             >  1:244195/1‑250 (MQ=255)                        aGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATg                                                                                                                                                                                                                          >  1:95254/1‑251 (MQ=255)                                            cGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCa                                                                                                                                                                                                        <  1:547576/249‑1 (MQ=255)                                                  tCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTg                                                                                                                                                                                                >  1:130969/1‑251 (MQ=255)                                                   cACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTgg                                                                                                                                                                                               <  1:467776/251‑1 (MQ=255)                                                    aCGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGCGATACTTTTGTCTTTACGCGCGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGCCGTTTGGc                                                                                                                                                                                              >  1:561047/1‑251 (MQ=255)                                                      gTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGATTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGCCCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGCGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCgg                                                                                                                                                                                            >  2:22826/1‑250 (MQ=255)                                                       tACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTc                                                                                                                                                                                           >  1:174303/1‑251 (MQ=255)                                                            cTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTGCAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTATCTTTACGCGAGTGCACCCAAACCAATGCCTCCTCTTCAATACCCAGACGTTTGGCGTCTTCgg                                                                                                                                                                                      >  1:297611/1‑251 (MQ=255)                                                               ggcggcGCGCTGATCGGCGATCGGCTCGAGGCGAACGGCGCAGTATTTGTACTCCGGCGGTTTCGTAATCGGGAGTAAGATTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑TCGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTAt                                                                                                                                                                                   <  2:219345/251‑1 (MQ=255)                                                                     gctcTGATCGGCGGTCGGCTCTACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGGCTTTACGCGAGTGCACCCAAACCAATGTCTCATCTTCAAGACCCAGACGTTTGGCGTCTTCGGTATTGAttt                                                                                                                                                                             <  2:559037/248‑1 (MQ=255)                                                                            tcggcgatcggcTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCCCCACTGTAAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGAACTGCGCACT‑GGTGATAATTTTGTCTTTACGCGAGTGCCCCCAAACAAATGCCTCATCTTCAATACCCCGACGTTTGGCGTCTTCGGTCTTGTTTTTTGCGt                                                                                                                                                                       >  2:257877/1‑250 (MQ=255)                                                                              ggcgatcggcTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGc                                                                                                                                                                    <  2:112592/251‑1 (MQ=255)                                                                                  aacggatCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAgg                                                                                                                                                                <  2:4439/245‑1 (MQ=255)                                                                                     ggcTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTATCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGGTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATACTTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCAAGCCGTTTTGCGTCTTCGGTATTGATTTGTGCGTAGCCAggctt                                                                                                                                                             >  2:262491/1‑248 (MQ=255)                                                                                      gcTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTtca                                                                                                                                                            <  1:450523/251‑1 (MQ=255)                                                                                             gcgAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATcagcca                                                                                                                                                     <  1:22826/251‑1 (MQ=255)                                                                                             gcgAACGGCGCAGTATTTGTACTCCGGCGCTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATcagcc                                                                                                                                                      >  1:487524/1‑250 (MQ=255)                                                                                                   ggCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGcgccg                                                                                                                                               >  1:46471/1‑251 (MQ=255)                                                                                                             ttgaaCTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGc                                                                                                                                     >  1:541394/5‑251 (MQ=255)                                                                                                                     cGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTAc                                                                                                                             >  2:406292/1‑251 (MQ=255)                                                                                                                      ggCGTTTTCGTAAACGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCCCCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCCCG‑GGTGATAATTTTGTCTTTACTCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGCCGTTTGGCGTCTTCGGTATTGATTTTTGCGTATCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTAcc                                                                                                                            <  1:568813/251‑1 (MQ=255)                                                                                                                        cGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCgg                                                                                                                          >  1:432976/1‑251 (MQ=255)                                                                                                                               gTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCg                                                                                                                    >  1:214062/1‑250 (MQ=255)                                                                                                                                     gggCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAAGCCACCACTGGTAGGTAATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCa                                                                                                             <  1:7064/251‑1 (MQ=255)                                                                                                                                          tAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGTCGATTGGCGTCTTCGGTATTGATTTGTGCGTAGCCTGGCTCATCAGCAAGAGCCGCCAGTGCCGCACAGTTACCGCTCATCgcatggcgagagg                                                                                                        >  2:213009/1‑239 (MQ=255)                                                                                                                                                        ccAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTCCCGGTCATCGAACGGCAAGAGTAGTGACCAACTTc                                                                                           >  2:382170/1‑250 (MQ=255)                                                                                                                                                        ccAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGCACGGCAAGGGTCGTGGCCAACttac                                                                                          >  2:271460/1‑249 (MQ=255)                                                                                                                                                               gtgaaagccACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTGCGG‑ACGAGCGCGGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGGGCGTAGCCAGGTTCAACATCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCg                                                                                   <  1:557334/243‑1 (MQ=255)                                                                                                                                                               gTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCg                                                                                   <  1:53775/251‑1 (MQ=255)                                                                                                                                                                  acatgccccaactccccagtGGTAGTTCATGATAATCGCCCCTTTGTTCGGTAC‑TTCGCTGGTCTGCGCACG‑GGTTATAATTTTGTCTTTACGCGTGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTTGTCGGTATTGATTTGTGCGTAGCAAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGtt                                                                                 <  2:390778/231‑1 (MQ=255)                                                                                                                                                                  aaagcccccAATCCCCCTCGGGTAGGTCATGTAAATTGCCTCTTTGTTCGG‑ACGATCGCTGCCCTGCGAACG‑GGTGATATTTTTGTCTTTCCGCGAGTGCACCCAAACCAATGCCTCTTCTGCAATCCCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCTCGCACCGtt                                                                                 <  2:389784/243‑1 (MQ=255)                                                                                                                                                                          caaaccaccaCGGGTGGGACATGTAAACCGCCCCTTTGGTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCTTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTAcc                                                                        <  1:452677/247‑1 (MQ=255)                                                                                                                                                                                accacTGGTAGGTCATGTAAATCGCCCCTGTGTTCGG‑ACGATCGCTGCCCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCggg                                                                  <  1:262491/251‑1 (MQ=255)                                                                                                                                                                                            tCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTg                                                      >  1:83449/1‑251 (MQ=255)                                                                                                                                                                                             cATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGCGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTg                                                      >  2:300498/1‑250 (MQ=255)                                                                                                                                                                                                  aaaTCGCCCTTTTTTGCGG‑ACGATCGCTGCCTTGCGC‑CGTGGTGATAATTTTGTCTTTACGCTAGTGCCCCCAAACCTATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTTTGCGTAGCCAGGTTCATCAGCCAGCTCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGAACAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTg                                                <  2:508488/251‑1 (MQ=255)                                                                                                                                                                                                    aTCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCAc                                                                                         <  1:444261/208‑1 (MQ=255)                                                                                                                                                                                                    aTCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGGGACCAACTTCAc                                                                                         >  2:444261/1‑208 (MQ=255)                                                                                                                                                                                                       gacccTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCACTACCCAGACGTTTGGCGTCTTCGGTCTTGGTTTGTGCGTAGCCAGGTTCATCAGCCAgcg                                                                                                                                                  >  2:430214/3‑148 (MQ=255)                                                                                                                                                                                                       gCCCCTTTGTGCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAATCCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGCTTCATCAGCCAgcg                                                                                                                                                  <  1:430214/148‑1 (MQ=255)                                                                                                                                                                                                           cTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCCCCGTTGCCAGTACCATCGGGTACTCGTCGGTGGGTTTGTCGTTTGgc                                       >  2:252340/1‑251 (MQ=255)                                                                                                                                                                                                             ttGTTCGT‑ACGAGCGCTGACCTGCGCACG‑GGTGATAATTTTGTCGTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGgcg                                      <  1:564440/250‑1 (MQ=255)                                                                                                                                                                                                             ttGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGCCAGTACCATCGGGTACTCGTCGGTGAGATTGTCGCTTGgcgc                                     >  1:395921/1‑251 (MQ=255)                                                                                                                                                                                                              tgttgcgggtgATCGCTGACCTGCGCCCG‑GGTGTGAATTTTGTCTTTACGCGAGTGCTCCCAATCCAATGCCTCATGTTCAAGACCCAGTCGTTTGGCGTCTTCGGTATGTATTTGTGGGTAGCCAGGGTCATCCGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCAGCGGGTACTCGTCGGTGAGTTTGTCGATTGgcg                                      <  2:561047/240‑1 (MQ=255)                                                                                                                                                                                                                        gATCGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCg                           >  1:141294/1‑251 (MQ=255)                                                                                                                                                                                                                           cGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGCACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCATCGGTGAGTTTGTCGATTGGCGCTACCCAGACGc                          >  2:227751/1‑249 (MQ=255)                                                                                                                                                                                                                           cGCTGACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGc                          >  1:131580/1‑249 (MQ=255)                                                                                                                                                                                                                               gACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAGCCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTg                     >  1:13270/1‑250 (MQ=255)                                                                                                                                                                                                                               gACCTGCGCACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTg                     <  1:252340/250‑1 (MQ=255)                                                                                                                                                                                                                                    gcgcTCG‑GGTGATAGTTTTGTCTTTCCGCGAGTGCACCAAAACCAATGCCTCATCTTAAATACCCAGACGTTTGGCGTCTTCGGTATTAATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCAATTGGCGCTACCCAGTCGCAGGTgaagaa                <  2:244195/250‑1 (MQ=255)                                                                                                                                                                                                                                      gcACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTGCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtcg                                                            <  1:524583/204‑1 (MQ=255)                                                                                                                                                                                                                                      gcACG‑GGTGATAATTTTGTATTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCCTCAGCCAGCGCCGACAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCCCGAACCGTTGACAGTACCATCGGGTACtcg                                                            >  2:524583/1‑204 (MQ=255)                                                                                                                                                                                                                                       cACG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAAAGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCAGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTAGGAGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGGGAAGAACTgg            >  1:209165/1‑250 (MQ=255)                                                                                                                                                                                                                                        aCG‑GGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGCGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTgcg           >  1:549073/1‑251 (MQ=255)                                                                                                                                                                                                                                         aggggggATAATTTTGTCTTTACGCGAGGGCCCCCAACCCAATGCCTCATCTTCAATTCCCAGACGTTTGGCTTCTGCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCAGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTg             <  2:46471/243‑1 (MQ=255)                                                                                                                                                                                                                                                ttAATTTTGCCTTTAAGCGAGTGCACCCCACCCAATGCCTCATCTTCAATACCCAGACTTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCCGCGCCGACAGTGCCGCACAGTTACAGATCATCGAACGGCAAGAGTAGTGACCAAATTCACGCACCGTTGACAATACCATCTGTTACTCGTCGGTGCGTTTGTAGATTGACGCTCCCAAGTCGccgtcgaagacctacgcccgacc    >  2:541252/2‑229 (MQ=255)                                                                                                                                                                                                                                                  aaTTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGt  <  1:227751/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    GCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGG‑ACGATCGCTGACCTGCGCACG‑GGTGATAATTTTGCCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGT  >  NC_000913/4297223‑4297712

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: ‑