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%	187.5 / NA	57	N307T (AAC→ACC)	fdoG	formate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (21/36);  total (21/36)

GAGGCAGTTTTATCGTGAGCACTGGTTTCTGCGATGTATTCGCAGACTTTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGTGGTTGGT  >  NC_000913/4084609‑4085183                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                           tatgCAGTTTTCTCGTGAGCATTGGTTTCTGCGATGTATTCTCAGACTTTCAGGAACGCGTCTGTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACCTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCggtggt                                                                                                                                                                                                                                                                                        <  2:227832/296‑1 (MQ=255)  aGGCAGTTTTATCGTGAGCACTGTTTTCTGCGATGTATTCGCAGACTTTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGtat                                                                                                                                                                                                                                                                                      <  1:2572/300‑1 (MQ=255)  aGGCAGTTTTATCGTGAGCACTGGTTTCTGCGATGTATTCGCAGACTTTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGtat                                                                                                                                                                                                                                                                                      <  1:1840/300‑1 (MQ=255)      aGTTTTACCGTGATCACTGGTTTCTTCGATGTATTCGCAGACTTTCATGAACGCGTCTTTTGGCGGACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGc                                                                                                                                                                                                                                                                                  <  2:56788/300‑1 (MQ=255)       gttgtATCGTGAGCACTGTTTTCTGCGATGTATTCGCAGACTTTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGcc                                                                                                                                                                                                                                                                                 <  2:596752/296‑2 (MQ=255)       gTTTTATCGTGAGCACTGGTTTCTGCGATGTATTCGCAGACTTTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCt                                                                                                                                                                                                                                                                                 <  1:17423/300‑1 (MQ=255)       gTTTTATCGTGAGCACTGGTTTCTGCGATGTATTCGCAGACTTTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCt                                                                                                                                                                                                                                                                                 >  2:592813/1‑300 (MQ=255)                  aGCACTGTTTTCTGCGATGTATTCGCAGACTTTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGTAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTAtt                                                                                                                                                                                                                                                                       <  2:106768/299‑1 (MQ=255)                   gCACTGGTTTCTGCGATGTATTCGCAGACTTTCAGGAACGCGTCTTGTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTcg                                                                                                                                                                                                                                                                     <  1:358836/300‑1 (MQ=255)                                             gACTTTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTcagcag                                                                                                                                                                                                                                           <  1:309353/300‑1 (MQ=255)                                               cTTTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGt                                                                                                                                                                                                                                         >  1:524405/1‑300 (MQ=255)                                                tttCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCAGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTa                                                                                                                                                                                                                                        >  2:369187/1‑300 (MQ=255)                                                      gAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTGGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCaa                                                                                                                                                                                                                                  <  2:422841/300‑1 (MQ=255)                                                            gTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTNTCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCNTAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGcc                                                                                                                                                                                                                            >  2:492570/1‑300 (MQ=255)                                                                 ttGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACa                                                                                                                                                                                                                       <  1:2903/300‑1 (MQ=255)                                                                                        aCCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGAGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGAGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAATACCGGaa                                                                                                                                                                                                >  2:48560/1‑300 (MQ=255)                                                                                            caTCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGaa                                                                                                                                                                                            >  2:404277/1‑300 (MQ=255)                                                                                            caTCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGAAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAAAGaa                                                                                                                                                                                            >  2:131378/1‑300 (MQ=255)                                                                                                        aGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCAgagt                                                                                                                                                                               <  1:580446/300‑1 (MQ=255)                                                                                                           gggAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCATGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCATGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCAtagtag                                                                                                                                                                             <  1:569922/299‑1 (MQ=255)                                                                                                            ggAAACGTGTTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTc                                                                                                                                                                           <  1:408947/300‑1 (MQ=255)                                                                                                            ggAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTc                                                                                                                                                                           <  1:409020/300‑1 (MQ=255)                                                                                                               aaCGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCATGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑c                                                                                                                                                                       <  1:35018/300‑1 (MQ=255)                                                                                                                                 gTTCCACACGCAGCGCNGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAANTAGGGGCATAGTAGTCAGNCTAccgccgccgggcc                                                                                                                                                        >  1:537231/1‑295 (MQ=255)                                                                                                                                 gTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGAAAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTAACACCACCGTACGCGt                                                                                                                                                     >  2:172067/1‑300 (MQ=255)                                                                                                                                      acacGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAg                                                                                                                                                 <  2:514934/299‑1 (MQ=255)                                                                                                                                            aGCGCGGATGTTGCAGCGTGGTATCGCGTTTGTCGAAGCCGTTTTCGTCCAGTTCATAAGGCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAGAGCCGTAATCCTCACGCACGATCATGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGTTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGgat                                                                                                                                          <  2:319369/300‑1 (MQ=255)                                                                                                                                            aGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGgat                                                                                                                                          <  1:600182/300‑1 (MQ=255)                                                                                                                                            aGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCAGCCGTACGCGTAAAGCGAGgat                                                                                                                                          <  1:152183/300‑1 (MQ=255)                                                                                                                                             gcgcGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGc                                                                                                                                                                                                                                                                                  <  1:618769/165‑1 (MQ=255)                                                                                                                                             gcgcGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGc                                                                                                                                                                                                                                                                                  >  2:618769/1‑165 (MQ=255)                                                                                                                                                      ttGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAAt                                                                                                                                >  2:553585/1‑300 (MQ=255)                                                                                                                                                      ttGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGCGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAAt                                                                                                                                >  1:222562/1‑300 (MQ=255)                                                                                                                                                       tGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATc                                                                                                                               <  2:517640/300‑1 (MQ=255)                                                                                                                                                       tGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATc                                                                                                                               <  1:592813/300‑1 (MQ=255)                                                                                                                                                           gCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGc                                                                                                                                                                                                                                                                                  >  1:625010/1‑151 (MQ=255)                                                                                                                                                           gCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGc                                                                                                                                                                                                                                                                                  <  2:625010/151‑1 (MQ=255)                                                                                                                                                           gCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCt                                                                                                                           <  1:378565/300‑1 (MQ=255)                                                                                                                                                                gTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGCTCACAATCAGCTTcgc                                                                                                                       >  1:240031/1‑299 (MQ=255)                                                                                                                                                                gTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTc                                                                                                                         >  1:638085/1‑297 (MQ=255)                                                                                                                                                                                                   aTAAGTCCAGGATGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGTCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCgg                                                                                     <  2:169151/298‑1 (MQ=255)                                                                                                                                                                                                      aGTCCAGGAGGATTTATCGTACTGGCGTTTTTCCGCGTCGGAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGTACGATCAGGCTTTCGGTGGTATATGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAAccc                                                                                <  2:307145/300‑1 (MQ=255)                                                                                                                                                                                                       gTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAAccc                                                                                >  1:107602/1‑299 (MQ=255)                                                                                                                                                                                                               ggaTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGAACGGGt                                                                        >  2:352348/1‑299 (MQ=255)                                                                                                                                                                                                                   ttATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCt                                                                   <  2:360730/300‑1 (MQ=255)                                                                                                                                                                                                                        gTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAgc                                                              <  1:511932/300‑1 (MQ=255)                                                                                                                                                                                                                                          gtcgtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATcac                                            <  1:359110/300‑1 (MQ=255)                                                                                                                                                                                                                                             gtAGCCGGGGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCTAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATcaccac                                         <  2:165683/300‑1 (MQ=255)                                                                                                                                                                                                                                               aGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTAAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCacga                                       <  1:475755/300‑1 (MQ=255)                                                                                                                                                                                                                                                      gAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCAAGACGAGGtt                                >  1:237642/1‑300 (MQ=255)                                                                                                                                                                                                                                                          aGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCCTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGc                             >  2:143517/1‑299 (MQ=255)                                                                                                                                                                                                                                                                 cTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGa                     >  2:518631/1‑300 (MQ=255)                                                                                                                                                                                                                                                                     aaaGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTNAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTc                 <  1:469866/300‑1 (MQ=255)                                                                                                                                                                                                                                                                         ccGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAANGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGAACGGGTGAGCTTCAGCGGCGTTACCGCCAATCACCACGACGAGGTTGGCGTTCTTGATGTCGAcc             >  2:188514/1‑300 (MQ=255)                                                                                                                                                                                                                                                                               tcctcacgcgtgattaggtcggCGGTGGTATAGGCTTCGGTGTATTCGCGTTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCTACCGGTTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGtgg       <  2:295404/280‑1 (MQ=255)                                                                                                                                                                                                                                                                               tCCTCGCGCACGACCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCTGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCACCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCATtgg       <  2:296252/300‑1 (MQ=255)                                                                                                                                                                                                                                                                               tCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACNGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGtgg       <  1:131378/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                    acgcacgATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGtggttggt  <  1:620602/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                           GAGGCAGTTTTATCGTGAGCACTGGTTTCTGCGATGTATTCGCAGACTTTCAGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAA‑TAGGGGCATAGTAGTCAG‑CTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACCCAGTGGTTGGT  >  NC_000913/4084609‑4085183

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