Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,084,903	T→G	N307T (AAC→ACC)	fdoG ←	formate dehydrogenase O subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,084,903	0	T	G	100.0%	324.6 / NA	92	N307T (AAC→ACC)	fdoG	formate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (46/46);  total (46/46)

CAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGAC  >  NC_000913/4084659‑4085147                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     cAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGta                                                                                                                                                                                                                                                 <  1:256383/250‑1 (MQ=255)      aCGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCt                                                                                                                                                                                                                                           >  1:102877/1‑251 (MQ=255)         cgTCTTTTGGCGTACCACAGAGGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCATCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTTTCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTCCTTGCGTTTTTCCTCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGGGGTATAGTCTTCg                                                                                                                                                                                                                                        <  2:444987/251‑1 (MQ=255)           tCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCggg                                                                                                                                                                                                                                      >  1:186804/1‑250 (MQ=255)              tttGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGgtgt                                                                                                                                                                                                                                    <  1:337385/250‑1 (MQ=255)                 ggCGTACCCCATATGTTTTCAACCCCATCTGGCTTGTAGCGGGAAACGGGCTGTTTCAGCAAGTTCCCCACGCAGCGCTGATGTTGCAGCGTGGTAGCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGATGGCTGTGGTATAGGCTTCGGTGTATTc                                                                                                                                                                                                                                <  2:426454/251‑1 (MQ=255)                      aCCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGt                                                                                                                                                                                                                           <  1:464061/251‑1 (MQ=255)                      aCCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGCGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGt                                                                                                                                                                                                                           >  1:191632/1‑251 (MQ=255)                       ccACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGt                                                                                                                                                                                                                           <  1:217006/250‑1 (MQ=255)                         acaGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGa                                                                                                                                                                                                                        <  1:22455/251‑1 (MQ=255)                           aGATGTTTTCAACCACAACTGGCGTGTAGCGGGAAACGTACTGTTTCAGCAAGTTCCACACGAAGCGCGGCTGTTGCAGCGTGGTATCGCGTTTGGCGCAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGCTCAGGCTGGCGGTGGTATATGATTCGGTGTATTCGCGGTTggtt                                                                                                                                                                                                                      >  2:139510/1‑248 (MQ=255)                             aTGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAAttt                                                                                                                                                                                                                    <  2:243054/251‑1 (MQ=255)                                      aCCACATCGGGCGTGGGGCGGGGAACGTGCTGTTTCAGCAAGTCCCCCACGCAGCGCGGATGTTGCAGCGTGGTATCGCTTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCAtt                                                                                                                                                                                                              <  2:389951/248‑1 (MQ=255)                                              tGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAggg                                                                                                                                                                                                   >  1:382942/1‑250 (MQ=255)                                               ggCGTGTAGCGGGAAACGTGCGTTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTGCGGTGTATTCGCGGTGGAATTTTTCATTGTTCAGCAGGt                                                                                                                                                                                                   <  2:261384/250‑1 (MQ=255)                                                gcgggtAGCGGTAAACGGGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGGTTTCGTCCAGTTCATA‑GTCCAGGAGGATTTATCGTACTGGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACa                                                                                                                                                                                                <  2:479852/247‑1 (MQ=255)                                                            aaaCGTGCTGGTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCt                                                                                                                                                                                     <  2:399576/251‑1 (MQ=255)                                                                  gCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTGCCGCGTCGTAGCCGGTTAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAg                                                                                                                                                                                <  1:380071/250‑1 (MQ=255)                                                                      tttCAGCAAGTTCCACCCGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGCGGATTTATCGTACTTGCGTTTTTCCGAGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGAGGGATAGGCGTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCCGCAGGTACAGCCATACGCCTGACAGGaa                                                                                                                                                                             >  2:426864/1‑249 (MQ=255)                                                                      tttCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGAATCCCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGCATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAggga                                                                                                                                                                             >  1:99930/1‑247 (MQ=255)                                                                          aGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAAt                                                                                                                                                                       >  1:24599/1‑251 (MQ=255)                                                                          aGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAAAAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAAt                                                                                                                                                                       >  1:216948/1‑251 (MQ=255)                                                                            cAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAACGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGACTTTTTCAtt                                                                                                                                                                                                              >  1:269625/1‑210 (MQ=255)                                                                            cAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGGTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGGGGTATAGGCTTCGGTGTATTCGTGGTTGAATTTTTCAtt                                                                                                                                                                                                              <  2:269625/210‑1 (MQ=255)                                                                               gTTCCACACGCAGCGCGGATGATGCAGCGTGGTATCGCGTTTGGCGAAGCCGCTTTCGTCCAGTTCATAAGTCCAGGAGGCTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACCGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGGGGGATTGGCTTCGGTGTATTCGCGGTTGACTTTTTCATTGTTCATCAGGTACAGCAATACGCCCGACAGGAAAGCAATGTc                                                                                                                                                                    >  2:257436/1‑249 (MQ=255)                                                                                ttCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCggtggt                                                                                                                                                                                                                                                  <  1:382788/170‑1 (MQ=255)                                                                                    acacGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCg                                                                                                                                                             >  1:267164/1‑251 (MQ=255)                                                                                        gCAGCGCGGATGTTGCAGCGTGGTAGCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCTCGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGaaa                                                                                                                                                         <  2:462033/251‑2 (MQ=255)                                                                                               ggatttttcatcGTGGTAGCGCTTTTGGCGAAGCCGTTCTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCATTACCGGAACGAATAgg                                                                                                                                                  <  2:277554/240‑1 (MQ=255)                                                                                               ggATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGGACGAATAgg                                                                                                                                                  >  1:508923/1‑251 (MQ=255)                                                                                                 agtttGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGGTCATAAGTCCAGGAGGATTTATCGTTCGTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGTATAgggg                                                                                                                                                <  2:180330/248‑1 (MQ=255)                                                                                                       cAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCAtagt                                                                                                                                          >  1:475561/1‑251 (MQ=255)                                                                                                        aGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAAAAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCAtagt                                                                                                                                          >  1:308105/1‑250 (MQ=255)                                                                                                         gcttGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTTCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGGGTATTCGCGGTTGAAGTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGCACGAATAGGGGTAtagt                                                                                                                                          <  2:252423/246‑1 (MQ=255)                                                                                                           gTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTc                                                                                                                                      >  1:402615/1‑251 (MQ=255)                                                                                                               tATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCt                                                                                                                                  >  1:430855/1‑251 (MQ=255)                                                                                                               tATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCt                                                                                                                                  >  1:368990/1‑251 (MQ=255)                                                                                                                     gTTTGGCGACGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGGGGTATAGGCTTAGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGGACAGCa                                                                                                                                                                                             >  2:169756/1‑186 (MQ=255)                                                                                                                     gTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCa                                                                                                                                                                                             <  1:169756/186‑1 (MQ=255)                                                                                                                        tGGCGAAGCCGTTTTCGTACCGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTCGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGCGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGGCAGCCATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgccgcc                                                                                                                         >  2:243678/1‑251 (MQ=255)                                                                                                                         ggCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCTGTTGAATTTTTCATTTTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAAGGTCAGTa                                                                                                                                                                <  2:303625/211‑1 (MQ=255)                                                                                                                         ggCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTAATTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTa                                                                                                                                                                >  1:303625/1‑211 (MQ=255)                                                                                                                         ggCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTAATTGCGTTTTTCCGCGTCGTAGACGGTGAACAGGCCATC‑TTAAAAGCCGTAATCCTC‑AAGCACGATCAGGCTGGCGGTGGGATAGACTTCGGGGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGAACCGGAACGAATAGCGGCAATGTAGTCAGCTAccgccgccg                                                                                                                        >  1:156570/1‑251 (MQ=255)                                                                                                                              aGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTCCTTCCGTCTTTACGCGTCGTAGCCGGTAAACAGGCAATC‑TTCAACGCCGTACTCCTC‑ACGCACCATCAGGCTGGCGGAGGTATAGGCATCGGTGTATTCGCGGATGACTTTTTCATTGTTCCGCAGGTACAGCCATACGCCTGACAGGAAAGCAATGTCAGTACCCGCACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAcgc                                                                                                                   >  1:328834/1‑251 (MQ=255)                                                                                                                               gCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAcgcg                                                                                                                  >  1:71965/1‑251 (MQ=255)                                                                                                                                 cGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCTTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACTATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTa                                                                                                                <  2:102877/251‑1 (MQ=255)                                                                                                                                      tCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGc                                                                                                            <  1:406206/250‑1 (MQ=255)                                                                                                                                        gTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAGAGCGAg                                                                                                         <  1:147003/251‑1 (MQ=255)                                                                                                                                           cAGTTCATAAGTCCAGGAGGATTTATCGTACTTGAGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCCTCGGTGAATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGGAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgat                                                                                                      >  1:506981/1‑251 (MQ=255)                                                                                                                                           cAGTTCATAAGCCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCGGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgat                                                                                                      <  2:267164/251‑1 (MQ=255)                                                                                                                                                cATAAGACAAGGCGGATTTCTCGTACTTGCGTTTTTCCGCGTCGTAGCCGGAGAACAGGCGATC‑TTCAAAGCCGCAATCCTC‑ACGCACGATCAGGCTGGCGGAGGACTAGGCTTCGGGGGCTTCGCGGTTGGATTTTTCATTGGTCCGCAGGCACCGCAATACGCCCGACAGGAAAGCAATGTAAGTAACGGAACGAATAGGGGCATAGGAGTCAGCTACCGCCGCCGTACGCGGAAAGCGAGgatcgatc                                                                                                 >  2:47224/1‑251 (MQ=255)                                                                                                                                                  tAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATcac                                                                                               >  1:34960/1‑251 (MQ=255)                                                                                                                                                   aaGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATcac                                                                                               <  1:67612/250‑1 (MQ=255)                                                                                                                                                       ccAGGCGGATTTATCGTACTTGCGTTTTTCCGCGTCGTCGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCATC‑ACGCTCGATCAGTCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTTTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCa                                                                                          <  2:186804/251‑1 (MQ=255)                                                                                                                                                       ccAGGAGGATTTATCGTTCTTGCTTTTTTCCTCGTCGTATCCGGTGATCAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTcagc                                                                                                                                                                                                       <  1:228403/142‑1 (MQ=255)                                                                                                                                                       ccAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCTTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGTAAAGCAATGTCAGTACCGGAACGAATAGGGGCAGAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCa                                                                                          <  1:221870/251‑1 (MQ=255)                                                                                                                                                         aggaggATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGGTCACAATCAGc                                                                                        >  1:319699/1‑251 (MQ=255)                                                                                                                                                          ggaggaTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGGTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGACAAGGGGCATAGTAGTCAGCTAACGCCGCCGTACGCGTAAAGCGAGGATCGATCACAAACAGCt                                                                                       >  2:288077/1‑251 (MQ=255)                                                                                                                                                             ggaTTTTTCGTACTTGCGTTTTCCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTTACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcg                                                                                    <  1:426864/251‑1 (MQ=255)                                                                                                                                                               aTTTATCGTACTTGCGTTTGTCCGCGTCGTAGCCGGTGAACAGGACATC‑TGCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGTTTCGGTGTATTCGCGTTTGAATTTTGCATTGTTCAGCAGGTACAGCAATACGCGTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAAGCGTAAAGCGAGGATCGATCACAATCAGCTTcgcg                                                                                  <  2:408595/251‑1 (MQ=255)                                                                                                                                                                        aCTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTtgtg                                                                          <  1:278836/250‑1 (MQ=255)                                                                                                                                                                           tGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAAt                                                                       <  2:319699/250‑1 (MQ=255)                                                                                                                                                                           tGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGC‑ATCCTTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGAAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAAt                                                                       <  2:402615/250‑1 (MQ=255)                                                                                                                                                                              gTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTTTATTCGCGGTTGAATTTTTCATTTTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGCAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTg                                                                   <  2:153626/251‑1 (MQ=255)                                                                                                                                                                                 tttCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTa                                                                                                                      >  1:441257/1‑197 (MQ=255)                                                                                                                                                                                 ttcccgcgccttAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGGGGTGGTATAGGCTCCGGTGTATCCGGGGTTAAATTTTTCATTGTTCAGAAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGAAACAAAAAGGGGAATAGTAGCCAGCTCCCGCCGCCGTa                                                                                                                      <  2:441257/186‑1 (MQ=255)                                                                                                                                                                                   cccGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTc                                                              <  1:259863/250‑1 (MQ=255)                                                                                                                                                                                    ccGCGTCGTAGCCGGCGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAACCGAGGATCGATAAAAACCAGCTTCGCGCCGTTGTGAATTTTGGCTTcc                                                             >  1:399900/1‑251 (MQ=255)                                                                                                                                                                                     cgcgTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCa                                                            <  1:341063/251‑1 (MQ=255)                                                                                                                                                                                     cgcgTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGGACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTcc                                                             >  2:358069/1‑250 (MQ=255)                                                                                                                                                                                        gtcgtAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGAGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCg                                                         >  1:99963/1‑251 (MQ=255)                                                                                                                                                                                                    gAACAGGCCATC‑TTAAACCCCGTAATCGTCCAC‑CACGAACAGGATGGCGGCGCTCTCGGTTTCGGTGACTTCGCGGGTAACTTCTTCATTGATCCGCAGCTTCAGCAATACGCCTGACAGGATAGCAAGGTCAGTACCGGAACGCATAGGGGAATCTGAGTCAGCTACCGCCGCCGTCCGCGTCAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTAGTGAATTTTGGCTTCAATCGCCCAGCgaatac                                             >  2:321639/1‑246 (MQ=255)                                                                                                                                                                                                     aaCAGGCAATC‑TTCAAAGCCGTAATCCTC‑ACACACGATCAGGCTGGCGGCGGTATAGACTTCGGTGTATTCGCTGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGGCAGGCAAGCAATGTCAGTACCGGAACGAATACGGGCCTAGAAGACAGCCACCGCCGCCGTACGCGTAAAGCGAGGAACGATCACAATCAGCTTCGCGCCGATGTGAATTTTGGCTTCCATCGCCCAGCGGAAcac                                            >  2:356136/1‑249 (MQ=255)                                                                                                                                                                                                      aCAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcgc                                                                                   <  1:202358/211‑1 (MQ=255)                                                                                                                                                                                                       cAGGCCCTC‑TTCAAAGCCGTAATCCTC‑ACGCACGCTCAGGCTGGCGGTGGCATAGGATTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCCATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCACCGCCGTACGCGTAAAGCGAGGATCGATCACAATCCGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCCGCGGAAACCgc                                          >  2:336461/1‑250 (MQ=255)                                                                                                                                                                                                        aGGCCTTC‑TTCAAAGCCGTAATCCCC‑TCGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTAGTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGGACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGAc                                         <  2:382942/251‑1 (MQ=255)                                                                                                                                                                                                         ggCCATT‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGAcc                                        >  1:139856/1‑251 (MQ=255)                                                                                                                                                                                                         ggCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGAcc                                        >  2:4040/1‑251 (MQ=255)                                                                                                                                                                                                               cctgaaaaGCCGTAATCCTC‑ACGAACGATCAGGCTGGCGGGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTGCATTGTTCAGCAGGTACAGCAATACGTCTGACAGGAAAGCAAGGCCAGTAACGGACCGAATGGGGGCAGAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGTCCAGCGGAACCCGACCGGGTg                                   <  2:222235/246‑1 (MQ=255)                                                                                                                                                                                                                  tCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTTTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCATCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGt                                    <  2:24599/247‑1 (MQ=255)                                                                                                                                                                                                                   aaaaGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATGGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCt                               <  2:71965/250‑1 (MQ=255)                                                                                                                                                                                                                      aGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTAAAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTc                             >  1:434075/1‑250 (MQ=255)                                                                                                                                                                                                                         cGTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAACGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAgcg                         <  1:162692/251‑1 (MQ=255)                                                                                                                                                                                                                          gTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCa                                                                                                                                                                         >  1:494585/1‑106 (MQ=255)                                                                                                                                                                                                                          gTAATCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGGATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCa                                                                                                                                                                         <  2:494585/106‑1 (MQ=255)                                                                                                                                                                                                                            aaTCCTC‑ACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACCGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAgcggcg                      >  1:516302/1‑251 (MQ=255)                                                                                                                                                                                                                            aaTCCTC‑ACGCACGATCAGGCTGGCGGTGGAATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAgcggcg                      >  1:414069/1‑251 (MQ=255)                                                                                                                                                                                                                                     cgcacgATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGcc              >  1:264734/1‑251 (MQ=255)                                                                                                                                                                                                                                        ccgATCAGGCTGGCGGTGGTATAGGCTGCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCAt           <  1:139510/250‑1 (MQ=255)                                                                                                                                                                                                                                           aTCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATccc        >  1:369384/1‑249 (MQ=255)                                                                                                                                                                                                                                             cAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACACTCAGCTTCGCGCCGTTGTGAATTTTGGATTCCCTCGCCCAGCGGAACCCGACCGGGGGAGCTTCCGCGGCCTTACCGCCCATccccc      >  2:276370/1‑247 (MQ=255)                                                                                                                                                                                                                                              aGGCTGGCGGAGGTATAGGCTTCGGTCTATACGCAGTTGAATTTTTCATTGGTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTTAGTACCGGAACGAATAGGGGCCTAGTAGTCAGCTACCGCCGCCGTACGCGTAAAACGAGGATCGATAACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCCGCGGCGTTACCGCCCATcaccac     >  2:169239/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ggCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAAGGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCacgac  <  2:34960/249‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     CAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATC‑TTCAAAGCCGTAATCCTC‑ACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGAC  >  NC_000913/4084659‑4085147

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