Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,304,647	T→C	I391V (ATT→GTT)	mtr ←	tryptophan:H(+) symporter Mtr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,304,647	0	T	C	100.0%	258.8 / NA	74	I391V (ATT→GTT)	mtr	tryptophan:H(+) symporter Mtr
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (44/30);  total (44/30)

CTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAATCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCG  >  NC_000913/3304405‑3304873                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                   cTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCa                                                                                                                                                                                                                             >  1:307181/1‑250 (MQ=255)         gTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGc                                                                                                                                                                                                                    >  1:31211/1‑251 (MQ=255)              aGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGGACCAGCGCAATCATCGGCTTGcc                                                                                                                                                                                                               <  1:243366/251‑1 (MQ=255)                  ggCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTAGCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCAccc                                                                                                                                                                                                           <  2:101244/251‑1 (MQ=255)                               aaGGAAGAGTTGTTGGGATAGTGTGTCGTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAAt                                                                                                                                                                                               <  1:230804/250‑1 (MQ=255)                                      agTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCt                                                                                                                                                                                       <  1:278502/251‑1 (MQ=255)                                             gggATAGTGTTTTCTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATATACCTTATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCTCTGATTACTGATACTCCTGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGGTGCCGTCGCCAAACACCAGATCCAGCGCAATCATCTGCTTGCCGCCCGAGACTCGGAATGTCGGGGTGCCaaa                                                                                                                                                                                 <  2:351239/250‑1 (MQ=255)                                             gggATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCaa                                                                                                                                                                                  >  2:234401/1‑249 (MQ=255)                                             gggATAGTGTGGCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTTCATAGACCTGATAAGCGAAGCGCAGCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACGCCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGAGGCGGAATTTCGGGCTGCCaaa                                                                                                                                                                                 <  2:535381/250‑1 (MQ=255)                                              ggATAGTGTGTTATAACTGGGGATGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGAAAAGCGAAGCGCATCAGGCATTGTGGAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGAGAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGTACCAGCGCAATCATCGTCTTGCAACCCCAGACGCGGAATTTCGGGCTGCCAAAGc                                                                                                                                                                               <  2:257339/251‑1 (MQ=255)                                                aTAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCCGACGCGGAATTTCGGGCAGCCAAAGCGt                                                                                                                                                                             >  2:354044/1‑251 (MQ=255)                                                  aGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCTCCCCAGACGCGGAATTTCGGGCTGCCAAAGCGtt                                                                                                                                                                            <  2:148378/250‑1 (MQ=255)                                                   gtgtgtCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGAGAAGCGAAGCGCATCAGGCATTGTGTAGGAAGCAGAAATGTCGGATAAGGCCCCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGCACCAGCGCAATCATCGGCTGGACACCCCAGACGCTGAATGTCGGGCTGCCCAAGCGttt                                                                                                                                                                           >  2:481491/1‑250 (MQ=255)                                                   gtgtgtCTTAACGCGGGAAGCCTTACCCGAGATGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCAGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCTGCGCAAGCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTa                                                                                                                                                                          <  1:544441/251‑1 (MQ=255)                                                     gtgtCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCACCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACCCCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTAc                                                                                                                                                                         <  2:250393/250‑1 (MQ=255)                                                       gtCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGaa                                                                                                                                                                                                <  1:585966/225‑1 (MQ=255)                                                       gtCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGaa                                                                                                                                                                                                >  2:585966/1‑225 (MQ=255)                                                           tAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATg                                                                                                                                                                   <  1:322763/250‑1 (MQ=255)                                                             aCGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACATGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACCCCGGCAGTAAATTAAAGCTCGATAAAAGATGCACCAGTGCGTTGCCGTCGCCAAACACCAGAACCAGCGCAATCAACGGCTTGCCACCCCAGCCTCGGAATTTCGGGCTGCCAAAGCGGTTACGCGATGccc                                                                                                                                                                <  2:133734/251‑3 (MQ=255)                                                                 gagAAGGCTTATCCGCGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGCACTGTCGGATAAGGCAACGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGACTTTAGGGCTGCCAAAGCGTTTACGCGATGCACGGGc                                                                                                                                                            >  2:175656/3‑251 (MQ=255)                                                                  ggAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCt                                                                                                                                                           >  1:549177/1‑251 (MQ=255)                                                                   gAAGCCTTATCCGAGCTGGCAACGCTGTCCTACAGAGACCTGATAAGCGAATCGCATCAGGCAGTGTGTAGGCAGCAGAAATGGCGGATAAGGCACCGCTGATTACTGAGACACCGGCAGTAAATTAAAGCTCGATAAAAGATGCACCAGTGCGTGGCCGACGCCAAGCACGAGAACCAGCGCAATCATCGGCTTGCGACCCCAGACGCGGAATTTGGGGGTGGCAAAACGTTTAccggatgtacaggcta                                                                                                                                                          <  2:530075/251‑16 (MQ=255)                                                                    aaGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATACGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTACAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCAGGCGGCCAAAGCGTTTACGCGATGCACGGGCTa                                                                                                                                                          >  2:52262/1‑250 (MQ=255)                                                                        cTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGACTTTCGGGCTGCCAAAGCGTTTACGCGGTGCACGGGCTAACAgc                                                                                                                                                     >  2:230229/1‑251 (MQ=255)                                                                          tATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAgcgc                                                                                                                                                   >  1:230149/1‑251 (MQ=255)                                                                               gAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTAGGATAAGGCACCGCTGATTACTGATACACCGGCAGTCAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGa                                                                                                                                               >  1:551045/1‑250 (MQ=255)                                                                                 gCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGGCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCGCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGTTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGTGCCTGAACa                                                                                                                                            <  2:295833/251‑1 (MQ=255)                                                                                    ggCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGGCGCCAAACACCAGAACCAGAGCAATCATCGGCTTGGCACCCCAGACGCGGAATTTCGGGATGCCAAAGCGTTTACGAGATGCACGGGATAACAGCGCCGGAACAAtt                                                                                                                                         >  2:239661/1‑251 (MQ=255)                                                                                     gCAACGCGGTCCTACATAGACATTATAAGTGAAGCGAATCAGGCCTTGTGTAGGCAGCAAACATGTCGGGTAAGGCACCGATGATTACTGATACACCGGCAGTAAATTACAGCTAGATATACTATGCACCAGTGCGTTGCCGACACCAAACAACAGAACCAACGCAATAATCGGCTTGCCACCCCAGACGCGGAATTGCGGGCAGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGGACaactg                                                                                                                                        >  2:570861/1‑248 (MQ=255)                                                                                         cGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCAACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTgccgc                                                                                                                                    >  1:367441/1‑251 (MQ=255)                                                                                             gTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATACAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGccc                                                                                                                                  >  1:243310/1‑249 (MQ=255)                                                                                                      aGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCATCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGcc                                                                                                                       <  2:215243/251‑1 (MQ=255)                                                                                                         ccTGATACGCGAAGCGCATCAGGCGTTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGAGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCt                                                                                                                    >  2:345951/1‑251 (MQ=255)                                                                                                               aagcgaagcgCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATGAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCa                                                                                                              <  1:354044/251‑1 (MQ=255)                                                                                                                agcgaagcgCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATGTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAg                                                                                                             <  2:178470/251‑1 (MQ=255)                                                                                                                 gcgaagcgCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGc                                                                                                            <  1:264924/251‑1 (MQ=255)                                                                                                                   gaagcgCATCAGGCATTGAGTAGGCAGCAGATATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAATTAAGCACCAGAGCGTTGCCGACGCCACACACCAGCACCAGCGCAATCCTCGGCTTGCCACCCCCGACGCGAACTTTCGGGCTGCCACAGCGTTTACGCGATGCACGGGCAAACAGCGCAGGAAAAGTTACCGCCACGAATGTAGCCGCAAAAACAGCAt                                                                                                          >  2:221653/1‑251 (MQ=255)                                                                                                                        gcATCAGGCATTGGGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCAAGATGGTAGCCGCTAAACCAGCATAACCa                                                                                                     <  2:496130/251‑1 (MQ=255)                                                                                                                                cATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACTCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACACCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATAGCGt                                                                                              >  1:251728/1‑250 (MQ=255)                                                                                                                                   tgtgTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTAc                                                                                            >  1:227704/1‑249 (MQ=255)                                                                                                                                    gtgtAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGCACAATTGCCGCCCAGATGGTAGCCGCGAAAACAGCATAACCAAAGACGTACAgg                                                                                         >  2:158110/1‑251 (MQ=255)                                                                                                                                     tgtAGGCAGCAGAAATGTCGGCTAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATACAATATGCACCAGTGCGTTGCCGACGCCAAACACAAGAACCAGCGCAATCATCGGCTTGCCACCCCCGCAGCGGAATTTCGGGCTGACAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCACAGATGGTAGCCGCTACACCAACATAAACAACAGCGTACAgg                                                                                         >  2:89831/1‑250 (MQ=255)                                                                                                                                     tgtAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGATAAACCAGCATAACCAATGGCGTACAgg                                                                                         >  1:411347/1‑250 (MQ=255)                                                                                                                                        agccagcagAAAAGAAGGATAATGCACCGCTGATTACTGATACCCCGGCAGAAAAGTAAAGCTCGATAAAATATGCCCCAGTGCGTTCCCCACGCCAAACACCAGAACGAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAAGTCCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCACAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAAt                                                                                      <  2:119552/247‑1 (MQ=255)                                                                                                                                         ggCAGCAGAAATGACGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGACCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATcc                                                                                    >  1:328196/1‑251 (MQ=255)                                                                                                                                                   aTGTCGGATAAGGCAACGCTGCCTACTGATACACCGGCAGAAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGACGACACCAAACACCAGAACCAGCGCCATCATCTGCTTGCCACCACAGACGCGGAATTTCGGGCTGCCAAATCGTTTACGCGATGCACAGGCTACCAGCGCCGGAACAATTGCCGCCCAGATGGGAGCAGCTAAACAAGCATAAACAATGGCGTACAGGAATCCGTTCGGTaac                                                                          >  2:309388/1‑251 (MQ=255)                                                                                                                                                        ggATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTAGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGAGCCGGAACAATTGCCGCCCAGATGGCAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGaacaac                                                                       >  1:358392/1‑249 (MQ=255)                                                                                                                                                          aTAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGcc                                                                   >  1:484848/1‑251 (MQ=255)                                                                                                                                                               gCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGccccccc                                                              >  1:57005/1‑251 (MQ=255)                                                                                                                                                                cACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGccccccc                                                              >  1:159276/1‑250 (MQ=255)                                                                                                                                                                             tGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGc                                                >  1:107188/1‑251 (MQ=255)                                                                                                                                                                                 acacCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGGGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGGAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACGGGCGGGGCAAAg                                            <  2:307181/251‑1 (MQ=255)                                                                                                                                                                                    ccGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTc                                         <  1:122633/251‑1 (MQ=255)                                                                                                                                                                                    ccGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTc                                         <  2:227704/251‑1 (MQ=255)                                                                                                                                                                                       gCAGTAAATTAAAGCTCGATAAAATATGCCCCCGTGCGTTTGCGACGCCAAACACCAGAACCAGCGCTATCATCGGCTTGCCCCCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACATCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACACCATCCCCCCCACAACTGGCGGTGCAAAGGTCAGc                                      <  1:155775/251‑1 (MQ=255)                                                                                                                                                                                       gCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGc                                      >  1:52587/1‑251 (MQ=255)                                                                                                                                                                                               ttAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGGACGCAATGCCGt                              >  2:468040/1‑251 (MQ=255)                                                                                                                                                                                                      tCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGAAATGCCGTTTTCa                         >  1:540727/1‑249 (MQ=255)                                                                                                                                                                                                        gATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCg                     >  1:474098/1‑251 (MQ=255)                                                                                                                                                                                                          aaaaaaaTGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACATCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGc                   <  2:334563/250‑1 (MQ=255)                                                                                                                                                                                                                atGCACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGACCAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCcac               >  2:323624/1‑249 (MQ=255)                                                                                                                                                                                                                atGCACCAGTGCGTTGCAGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCCCCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATCCACAGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACACCTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCACACGGCCcac               >  2:270826/1‑249 (MQ=255)                                                                                                                                                                                                                 tGCACCAGTGCGTTGCCGACGCCACACACCAGAACCAGCGCACTCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCAGCCAAAGCGTTTACGCGCTGCACGAGCTAACCGCGCCGGAACAATTGACGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGATTCCAAGCGGCCcacacc            >  2:32324/1‑249 (MQ=255)                                                                                                                                                                                                                  gctccAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGcc           <  1:564396/248‑1 (MQ=255)                                                                                                                                                                                                                   cACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATgg                                                                                                                            >  1:358805/1‑137 (MQ=255)                                                                                                                                                                                                                   cACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCg          >  1:410834/1‑251 (MQ=255)                                                                                                                                                                                                                   cACCAGTGCGTTGCCGACGCCAAACACCAGAACCAGCGCAAGCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCGAACAGCGCCGGAACAATTGCCGCCCAGATgg                                                                                                                            <  2:358805/137‑1 (MQ=255)                                                                                                                                                                                                                   aaCCAGTGAGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGATAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCACAccgccg          >  2:414648/2‑251 (MQ=255)                                                                                                                                                                                                                       aGTGCGTTGCCGACGCCAAACCCCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAgtc      <  1:328529/251‑1 (MQ=255)                                                                                                                                                                                                                         tGCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAgtcgt    >  1:192739/1‑251 (MQ=255)                                                                                                                                                                                                                          gCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAgtcgtc   >  1:504484/1‑251 (MQ=255)                                                                                                                                                                                                                          gCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAgtcgtc   >  1:281108/1‑251 (MQ=255)                                                                                                                                                                                                                          gCGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCAGCTTTCACGCGGCCCCCAGCCGAgccgtc   >  1:342526/1‑251 (MQ=255)                                                                                                                                                                                                                           cGTTGCCGACGCCAAACACCAGAACCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGtcgtcg  <  2:220979/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                   CTGCATTCGTTTGAGCTGGCGATGGATGTTAAGGAAGAGTTGTTGGGATAGTGTGTCTTAACGCGGGAAGCCTTATCCGAGCTGGCAACGCTGTCCTACATAGACCTGATAAGCGAAGCGCATCAGGCATTGTGTAGGCAGCAGAAATGTCGGATAAGGCACCGCTGATTACTGATACACCGGCAGTAAATTAAAGCTCGATAAAATATGCACCAGTGCGTTGCCGACGCCAAACACCAGAATCAGCGCAATCATCGGCTTGCCACCCCAGACGCGGAATTTCGGGCTGCCAAAGCGTTTACGCGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCG  >  NC_000913/3304405‑3304873

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