Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,761,230	G→A	S7F (TCC→TTC)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,761,230	0	G	A	100.0%	191.6 / NA	57	S7F (TCC→TTC)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/29);  total (28/29)

GTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGGCACAGGCAATGGCGATAGTCGCCAGATTGACCGTATCGGTTTTGAGT  >  NC_000913/3760952‑3761521                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                    gTAAAGCTCGCCCAAGTTTGGTATGCAGCACTGTTCCCGTCAGGTTGATCACCGGACGCAGCGCACTCTGCGCTTCTTTCGTCAACCGGGCATTGACTTCTTGCGCCCAGTTTTCACACCACGCATGCGGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGTAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAg                                                                                                                                                                                                                                                                                <  2:151404/300‑1 (MQ=255)               aGGTTGGTATGCAGCACGGTTCCCTTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGc                                                                                                                                                                                                                                                                      <  2:42791/296‑1 (MQ=255)               aGGTTGGTATGCAGCACGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAt                                                                                                                                                                                                                                                                  >  1:65903/1‑300 (MQ=255)                     gTATGCAGCACGGTTCCCTTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGc                                                                                                                                                                                                                                                            <  2:11551/300‑1 (MQ=255)                       aTGCAGCACGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCAAGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGccc                                                                                                                                                                                                                                                          >  1:417022/1‑300 (MQ=255)                         gcagcaCGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCctc                                                                                                                                                                                                                                                        >  1:455407/1‑300 (MQ=255)                          cagcaCGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAACGAAAGGAAGGAGCTATCGCGCACTAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGATTGTCATAGGAAACCTGGCTGATCAAgaccatct                                                                                                                                                                                                                                                       >  2:134395/1‑293 (MQ=255)                          cagcaCGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAGGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGCAACCTGTCTGATCACGGCCctct                                                                                                                                                                                                                                                       >  2:133627/1‑300 (MQ=255)                          cagcaCGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTAGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATTAATAGACGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGTCTGAt                                                                                                                                                                                                                                                                  >  2:133736/1‑289 (MQ=255)                              actgTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAATAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGTAAATGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTcac                                                                                                                                                                                                                                                    <  2:95959/296‑1 (MQ=255)                               cGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGCCTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAACGATCAATAGTCGGAAGTTGACTATAGAGGAAACACGTTTCGGTTGTAATAGGAAACCTGGCTGATCAAGGCCCTCTcaca                                                                                                                                                                                                                                                   >  2:35215/1‑299 (MQ=255)                               cGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGACGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATTAAGGCCCTCTaaca                                                                                                                                                                                                                                                   >  2:35264/1‑299 (MQ=255)                                 gTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCNGGCTGCCACGAATCACTTCTCGCGCTTCGTCGNGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTATCacac                                                                                                                                                                                                                                                  >  2:455831/1‑298 (MQ=255)                                   tCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGgag                                                                                                                                                                                                                                              <  1:414366/300‑1 (MQ=255)                                        tCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGgaga                                                                                                                                                                                                                                             <  2:601481/296‑1 (MQ=255)                                            gTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGtt                                                                                                                                                                                                                                     >  2:548713/1‑300 (MQ=255)                                                   aCCGTACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGTAGAAGGGCGTTTAACAt                                                                                                                                                                                                                               <  2:111283/299‑1 (MQ=255)                                                      ggACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAAc                                                                                                                                                                                                                            <  1:202024/299‑1 (MQ=255)                                                       gACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCa                                                                                                                                                                                                                          <  2:304601/300‑1 (MQ=255)                                                                 gcTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCAGGAAGTTGACTATNGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCCCACGGAGAAGGGCGTTTAACATAACCACGGATTGTaa                                                                                                                                                                                                                >  2:482590/1‑300 (MQ=255)                                                                 gcTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCAAGCAGGCAGCGTCTGGCTGCCGCGAATCACTTCTCGCGCTTCGTCGAGCACCTGAAGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGAAAAGAAAGGAAGGAGCTATAGCGCAATAAGCGATCAATAGCCGGAAGTCGACTACAGAGGAAACGAGTTTCGGTTGTAATAGGAAACCTAGCTGAta                                                                                                                                                                                                                                                                 >  2:173494/1‑250 (MQ=255)                                                                      gcgcTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAGCAATTCCACCACGCGGGTGTGACCATAAGTATCAAGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCAGGAAGTTGGGTATAGAGGAAACGTGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACATGGAGAAGGGCGTTTAAAATAACTACGGATTGTAACGTga                                                                                                                                                                                                           >  2:50681/1‑300 (MQ=255)                                                                                    aaCCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCaggag                                                                                                                                                                                              <  2:363136/299‑1 (MQ=255)                                                                                        gggCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACat                                                                                                                                                                                         <  1:537805/300‑1 (MQ=255)                                                                                                       gcgcCCAGTTTTCACACCACGCAGGCAGCGTCTTGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGTTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGcc                                                                                                                                                                          <  2:74606/300‑1 (MQ=255)                                                                                                             aGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGTTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGc                                                                                                                                                                    <  2:96818/300‑1 (MQ=255)                                                                                                              gTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAAACACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGc                                                                                                                                                                    >  2:427612/1‑299 (MQ=255)                                                                                                              gTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTAGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGACGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGc                                                                                                                                                                    >  2:24632/1‑299 (MQ=255)                                                                                                              gTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTATCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCAAGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCAGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTATCACACGGAGAAGGGCGTTTAACATAACCACGGAGTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGc                                                                                                                                                                    >  2:24622/1‑299 (MQ=255)                                                                                                                          cGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTAGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGANAGGAAGGAGCTATCGCGCAATAAGCGATCAATNGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGAAATATCGCGCATCAAGCATTTGGCGGTTCGGTGCGa                                                                                                                                                        >  2:470240/1‑299 (MQ=255)                                                                                                                                  gCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGATATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCaaaaatt                                                                                                                                               >  1:419840/1‑297 (MQ=255)                                                                                                                                        ggCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTcgc                                                                                                                                           >  2:590171/1‑298 (MQ=255)                                                                                                                                             ccACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTaaa                                                                                                                                    >  1:303715/1‑300 (MQ=255)                                                                                                                                                                  tcgtcgAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGCGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGAc                                                                                                               <  1:309911/300‑1 (MQ=255)                                                                                                                                                                  tcgtcgAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCAAGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGAc                                                                                                               >  2:553830/1‑300 (MQ=255)                                                                                                                                                                    gtcgAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTa                                                                                                             >  1:170762/1‑300 (MQ=255)                                                                                                                                                                    gtcgAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTa                                                                                                             <  2:416449/300‑1 (MQ=255)                                                                                                                                                                        aGCATCTGACGCAACAATTCCACCACTCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATNAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAATGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTCGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAgtg                                                                                                         <  1:125394/300‑1 (MQ=255)                                                                                                                                                                        aGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATNAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGTAGAAGGGCGTTTAACATAACCACGGATTGTAACGTTAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAgtg                                                                                                         <  1:125404/300‑1 (MQ=255)                                                                                                                                                                        aGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATNAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAgtg                                                                                                         <  1:126247/300‑1 (MQ=255)                                                                                                                                                                        aGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAgtg                                                                                                         <  1:179873/300‑1 (MQ=255)                                                                                                                                                                             cTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGc                                                                                                    >  1:389150/1‑300 (MQ=255)                                                                                                                                                                              tGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGTAAGGAGCTATCGCGCAATAAGCGAGCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCg                                                                                                   <  1:172127/300‑1 (MQ=255)                                                                                                                                                                                         ttCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACANAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACAcc                                                                                        <  2:452259/300‑1 (MQ=255)                                                                                                                                                                                                acgcgtgggggACCATAAGTATCACGCAAAGGAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCt                                                                                   <  2:321228/288‑1 (MQ=255)                                                                                                                                                                                                acGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGTAAACCTGGCTGATCAAGGCCCTCTCACACGTAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCt                                                                                   <  2:564420/298‑1 (MQ=255)                                                                                                                                                                                                    gggTGTGACCATAAGTATCACGCAAAGAAATGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTTTAACGTTAGATGGGTCATGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGgcgc                                                                              <  2:35019/299‑1 (MQ=255)                                                                                                                                                                                                          gACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCAGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAAcgc                                                                         >  1:160587/1‑298 (MQ=255)                                                                                                                                                                                                           aCCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACgcggc                                                                      >  1:412161/1‑300 (MQ=255)                                                                                                                                                                                                                  gTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCAGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTNATCAAGGCCCTCTAACCCGGAGAAGGGCGTTTNACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGTTAAACAGGTTTTCGAACCGTTTGACTAAGTGTGGGCGGTCGACACACAAACCTGGCGCAACGCGGCGGAAAtt                                                               >  2:459042/1‑300 (MQ=255)                                                                                                                                                                                                                       aCGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGat                                                          <  2:374426/300‑1 (MQ=255)                                                                                                                                                                                                                                        aGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGGCACAGGCa                                         <  1:593064/300‑1 (MQ=255)                                                                                                                                                                                                                                        aGCTATCGCGCAATAAGCGAGCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCGGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCTGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGGCACACGCa                                         <  1:236827/300‑1 (MQ=255)                                                                                                                                                                                                                                                   aaTAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGGCACAGGCAATGGCGATAGt                              <  1:568133/300‑1 (MQ=255)                                                                                                                                                                                                                                                                      gAAGTTGACTATAGAGGAACCGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTAACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATAGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGCTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGAGGAAATTGAGATATCCGAGGACACAGGCAATGGCGATAGTCGAAAGACTGGCCGTATc            >  2:159284/1‑299 (MQ=255)                                                                                                                                                                                                                                                                      gAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGGCACAGGCAATGGCGATAGTCGCCAGATTGACCGTATc            >  1:604010/1‑299 (MQ=255)                                                                                                                                                                                                                                                                               tataGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGGCACAGGCAATGGCGATAGTCGCCAGATTGACCGTATCGGTTTTGAgt  <  1:455831/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                    GTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGGCACAGGCAATGGCGATAGTCGCCAGATTGACCGTATCGGTTTTGAGT  >  NC_000913/3760952‑3761521

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