Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,298,199	C→T	A390T (GCA→ACA)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,298,199	0	C	T	98.2%	180.6 / NA	57	A390T (GCA→ACA)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/39);  minor base G (0/1);  total (17/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGAT  >  NC_000913/4297905‑4298490                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                    gggTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCgg                                                                                                                                                                                                                                                                                                 <  1:447251/299‑1 (MQ=255)   gTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGNCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCAATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGc                                                                                                                                                                                                                                                                                              <  1:189149/300‑1 (MQ=255)   gTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGc                                                                                                                                                                                                                                                                                              <  1:2761/300‑1 (MQ=255)    tGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGc                                                                                                                                                                                                                                                                                              <  1:94808/299‑1 (MQ=255)                tgatgaTTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCNGAACGATAACCAGTTCCAGATCTTCAAAGGCTNTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCTGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGa                                                                                                                                                                                                                                                                                 <  2:641449/300‑1 (MQ=255)                tgatgaTTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCNGAACGATAACCAGTTCCAGATCTTCAAAGGCTNTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGa                                                                                                                                                                                                                                                                                 <  2:641583/300‑1 (MQ=255)                     atTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGa                                                                                                                                                                                                                                                                            <  1:403844/300‑1 (MQ=255)                     atTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTGGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGa                                                                                                                                                                                                                                                                            <  1:403973/300‑1 (MQ=255)                                   tttCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATgcgc                                                                                                                                                                                                                                                              <  2:581157/300‑1 (MQ=255)                                   tgtCAGATCCCATTTCGGTTAAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATgcgc                                                                                                                                                                                                                                                              <  2:14190/298‑1 (MQ=255)                                        gATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGNCGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCa                                                                                                                                                                                                                                                         <  2:666441/300‑1 (MQ=255)                                             cATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCtt                                                                                                                                                                                                                                                    >  2:425057/1‑300 (MQ=255)                                             cATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTCCATGCTCGCCCCACGACGTTGACGGTAAAATAGCACATGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTTCAGATTTTCAAAGGCTTTACGTACTGACGACAGCTCCGCGTTAGTTTGTCGCGGCTCTTCGCCCATAATGTACGCGGCACGCACTTCGACATGCGCTGTGCGGTGCGGCAGCTCGCTGAt                                                                                                                                                                                                                                                                                >  2:425050/1‑272 (MQ=255)                                                   ggTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTTTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCAc                                                                                                                                                                                                                                               <  2:164176/299‑1 (MQ=255)                                                            gCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCt                                                                                                                                                                                                                                     <  1:356858/300‑1 (MQ=255)                                                            gCCTGGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCt                                                                                                                                                                                                                                     <  2:488998/300‑1 (MQ=255)                                                                                                  aaaCACGCCTTCATGCTCGCCCCACGACGTTGATGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTa                                                                                                                                                                                               >  1:590037/1‑300 (MQ=255)                                                                                                  aaaCACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGCCAGCTCCGAGTCAGTTTGTCGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGTGGTTAGCTGGATCTTTCaagt                                                                                                                                                                                                >  2:439718/1‑296 (MQ=255)                                                                                                  aaaCACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGAGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGTAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCAGGATCTTTCACGTa                                                                                                                                                                                               >  2:439586/1‑300 (MQ=255)                                                                                                  aaaCACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATATCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTAACGTa                                                                                                                                                                                               >  2:439738/1‑300 (MQ=255)                                                                                                              aTGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCg                                                                                                                                                                                    <  1:286481/299‑1 (MQ=255)                                                                                                                     ccccACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATg                                                                                                                                                                                                                                                                 <  2:140428/215‑1 (MQ=255)                                                                                                                     ccccACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTAGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTa                                                                                                                                                                            >  1:229359/1‑300 (MQ=255)                                                                                                                            acgTTGACGGTAAAATAACATCCGCCGCCTACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCNTGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGc                                                                                                                                                                      <  2:596276/299‑1 (MQ=255)                                                                                                                                          aTAACATCCTCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAACTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCg                                                                                                                                                       <  2:229359/300‑1 (MQ=255)                                                                                                                                          aTAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCg                                                                                                                                                       <  1:542879/300‑1 (MQ=255)                                                                                                                                               aTCCGCCGCCGACGCGGTTTGGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAgg                                                                                                                                                   <  2:566776/299‑1 (MQ=255)                                                                                                                                                          aCGCGGTTTTGGTCATAAGAACGTCCTGAACGATAACCAGTTCCAGATCGTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCg                                                                                                                                                       <  2:63544/284‑1 (MQ=255)                                                                                                                                                          aCGCGGTTTTGGTCATAAATATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTATGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTTCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCTGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCAccc                                                                                                                                              <  2:41542/293‑1 (MQ=255)                                                                                                                                                                           aaGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAAcc                                                                                                                      >  1:356871/1‑300 (MQ=255)                                                                                                                                                                             gATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGACATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATCTCGCAGGCACCCTGAACGTTGTTCTGACCACGAAcc                                                                                                                      >  2:583169/1‑298 (MQ=255)                                                                                                                                                                                gTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGtt                                                                                                                 >  2:834060/1‑300 (MQ=255)                                                                                                                                                                                    tGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAaca                                                                                                             <  2:482344/300‑1 (MQ=255)                                                                                                                                                                                                         tCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTTTATCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATTGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAgg                                                                                        <  1:549832/300‑1 (MQ=255)                                                                                                                                                                                                                       aCGTACTGCCGACAGCTCCTCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCAt                                                                          <  2:354341/300‑1 (MQ=255)                                                                                                                                                                                                                       aCGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCATGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCAt                                                                          <  1:491445/300‑1 (MQ=255)                                                                                                                                                                                                                             tGCCGACAGCTCCGCTTCAGTTTGTAGCGGATCTTCTCCCATAATGTACGCGACACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTCCTCGCGGTTAGCCGGATCTTTCATGTACGGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGa                                                                     <  2:356871/299‑1 (MQ=255)                                                                                                                                                                                                                                      cTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCATGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGc                                                                  <  1:440732/293‑1 (MQ=255)                                                                                                                                                                                                                                                          gATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGc                                       >  1:255445/1‑300 (MQ=255)                                                                                                                                                                                                                                                          gATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACAACCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGc                                       >  1:257206/1‑300 (MQ=255)                                                                                                                                                                                                                                                                     aTAATGTACGCTGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGa                              <  1:425050/298‑1 (MQ=255)                                                                                                                                                                                                                                                                     aTAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGa                              <  1:425057/298‑1 (MQ=255)                                                                                                                                                                                                                                                                      tAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCANGCTTTCCACGCCCCATGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACt                           <  1:611829/300‑1 (MQ=255)                                                                                                                                                                                                                                                                          gtacgataaacgcacTTCTCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGt                       <  2:455093/291‑1 (MQ=255)                                                                                                                                                                                                                                                                          gtacactgcacgaacTGCGCCATGCGCTGTGCGTTGCGGCAGCTTACGGATACGATAGCCGGTATGCGCGGGCAGGATGTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCTCGGATACGTATCCGGCAGCGCGCCCATATCACAGGCACCCTGAACGTTGTTCTGACCACGAATCGGGTTAACACCCGCATGCGGCTTACCGAGGTGACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGTTAGAACTGGGt                       <  2:455753/293‑1 (MQ=255)                                                                                                                                                                                                                                                                          gTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGt                       <  1:20146/300‑1 (MQ=255)                                                                                                                                                                                                                                                                          gTACGCGGCACCCACTTCGCCATGCGCTGTTCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTAGGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGTGAATTTCTCGCGGTTAGCCGAACCTTGCACGTACTGGTATCCCGGATACGTATCCGGTAGCGCGCTTATATCg                                                                                                                                                       <  2:63640/172‑1 (MQ=255)                                                                                                                                                                                                                                                                                   acgcacTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACNGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCAt              >  1:157091/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                   acgcacTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACNGAGGTTACCGGTCAGCATCGCGAAGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCAt              >  1:157105/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                       acTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCANGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGcc           <  1:173182/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                       acTTCGCCATGCGCTGGGCGGTGCGGCAACCAGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTGTCCACTCCCCAGGCTTTGGCGAATTGCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTAGCCGGCAGCGCGCCCATATCACAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGTTAGAACTGGGTTACACCCATGccc          <  2:456772/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                          tCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCcac       >  1:683947/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                            gCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGTGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGAAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGAATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTACACGCCCTGGTAGAACTGGGTTACAACCATGCCCcaca      >  2:194225/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                            gCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCcaca      >  1:680648/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                               aTGCGCTGTGCGGTGCGGCAGCTCTCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAgg    <  1:439586/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                               aTGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGAt  <  1:834060/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                               aTGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCCGGCAGGCTTTCCACGCCCCATGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCAGGAGGAGGAAAGAGGTTACCGGTCAGCATCGCTAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGTTAGAACTGGGTTACACCCATGCCCCACAgg    <  1:439738/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                                ggcgctgtgcGGTGCGGCAGCGCGCGGATGCGAGAGCCGGGATACGCTGTCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCACAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGAt  <  2:472572/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                    GGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGAT  >  NC_000913/4297905‑4298490

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