Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,449,779	G→A	S35S (AGC→AGT)	rpsC ←	30S ribosomal subunit protein S3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,449,779	0	G	A	100.0%	298.9 / NA	76	S35S (AGC→AGT)	rpsC	30S ribosomal subunit protein S3
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/37);  total (39/37)

AGCAACCAGTTTTGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAAC  >  NC_000913/3449533‑3450023                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     aGCAACCAGTTTTGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGt                                                                                                                                                                                                                                                  >  2:356626/1‑251 (MQ=255)  gCAACCAGTTTTGCGGCCCGTTACGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCCGCGATGTCCGCTACGACCTTACGCAGTTTGTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTACCTTAGTCAGGTACTGACGTACTTTAAAATCACTGtn                                                                                                                                                                                                                                                 <  2:289069/251‑2 (MQ=255)    aaCCAGTTTTGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTcc                                                                                                                                                                                                                                                >  1:374173/1‑250 (MQ=255)       cAGTTTTGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGt                                                                                                                                                                                                                                            <  1:287668/251‑1 (MQ=255)         gTTTTGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTg                                                                                                                                                                                                                                          >  2:106680/1‑251 (MQ=255)            ttGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTc                                                                                                                                                                                                                                        <  1:129998/250‑1 (MQ=255)             tGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCa                                                                                                                                                                                                                                       <  2:628647/250‑1 (MQ=255)              gCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAg                                                                                                                                                                                                                                      <  1:113181/250‑1 (MQ=255)                  ccAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAAt                                                                                                                                                                                                                                 >  1:319396/1‑251 (MQ=255)                    aGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTc                                                                                                                                                                                                                               >  2:380792/1‑251 (MQ=255)                            cTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGgtgt                                                                                                                                                                                                                       >  2:583034/1‑251 (MQ=255)                               aCGAGCTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTc                                                                                                                                                                                                                     <  2:104577/250‑1 (MQ=255)                               aCGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCg                                                                                                                                                                                                                    <  2:585970/251‑1 (MQ=255)                               aCGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGCACTTACGGACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTc                                                                                                                                                                                                                     <  1:626980/250‑1 (MQ=255)                                 gAACTTCGGCGATGTCGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGc                                                                                                                                                                                                                   >  2:78762/1‑250 (MQ=255)                                       cGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGGCAGCGAATTCTTTGGTGGTCGCAAACCa                                                                                                                                                                                                             <  1:150901/250‑1 (MQ=255)                                                  aTCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTcc                                                                                                                                                                                                  >  1:68699/1‑250 (MQ=255)                                                      gtgCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATgg                                                                                                                                                                                              <  2:555993/250‑1 (MQ=255)                                                          aGGAACGCCAGCGATGTCCGCTACGACCTTACGCCGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGtttt                                                                                                                                                                                          >  2:169661/1‑250 (MQ=255)                                                          aGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTa                                                                                                                                                                                         <  2:313519/251‑1 (MQ=255)                                                              aCGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAAt                                                                                                                                                                                     <  1:36737/251‑1 (MQ=255)                                                                    gCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATAccc                                                                                                                                                                                 >  1:495947/1‑249 (MQ=255)                                                                    gCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTCTACAATACCCAg                                                                                                                                                                               <  1:447520/251‑1 (MQ=255)                                                                                aCGACCTTACGCAGTTTTTCTTCGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCAt                                                                                                                                                                    <  2:595838/250‑1 (MQ=255)                                                                                aCGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCAtt                                                                                                                                                                   >  2:370590/1‑251 (MQ=255)                                                                                       tACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATg                                                                                                                                                             <  1:458468/250‑1 (MQ=255)                                                                                               gtttcatcttCTTCACCTTTTTTACCGAGAACGATACCCGGGCGAGCAATGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGGCGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCGTTGGTGTTCGCAAACCAGGTAGAGTTCCAGGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCt                                                                                                                                                     <  2:603830/241‑1 (MQ=255)                                                                                                     aCGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCAt                                                                                                                                              <  2:16300/251‑1 (MQ=255)                                                                                                       gTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTTTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTg                                                                                                                                            <  2:470614/251‑1 (MQ=255)                                                                                                       gTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTg                                                                                                                                            >  1:101278/1‑251 (MQ=255)                                                                                                       gTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTg                                                                                                                                            >  1:574928/1‑251 (MQ=255)                                                                                                          ttCACCTTTTTTACCGATAACGATACACGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCt                                                                                                                                          >  1:110978/1‑250 (MQ=255)                                                                                                           tCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTa                                                                                                                                         <  1:380792/250‑1 (MQ=255)                                                                                                            cACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGt                                                                                                                                       >  1:49218/1‑251 (MQ=255)                                                                                                            cACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGt                                                                                                                                       <  1:466564/251‑1 (MQ=255)                                                                                                             aCCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGtc                                                                                                                                      >  2:459809/1‑251 (MQ=255)                                                                                                              ccTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGgtgt                                                                                                                                                                                                                       <  2:421848/169‑1 (MQ=255)                                                                                                              ccTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGgtgt                                                                                                                                                                                                                       >  1:421848/1‑169 (MQ=255)                                                                                                              ccTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGCTGTCCTTTCTGACCCATTGCTAGtct                                                                                                                                     >  2:514698/1‑251 (MQ=255)                                                                                                                    ttACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCag                                                                                                                                 >  2:130298/1‑249 (MQ=255)                                                                                                                      aCCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGtc                                                                                                                             >  1:42178/1‑251 (MQ=255)                                                                                                                          aTAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGtctc                                                                                                                           >  1:220637/1‑249 (MQ=255)                                                                                                                              cGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGAt                                                                                                                     <  1:356626/251‑1 (MQ=255)                                                                                                                               gATACCCGGGCGAGCAGTGTGAATGGGTACACGGATGCTCTTAGCCGGAAGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTAAGAGAATTCTTAGGTGTTCGCAAACCAGGAAGAGATCCATGGTTTTACAATACCCAGGAGAATACTATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATc                                                                                                                    >  1:574182/1‑251 (MQ=255)                                                                                                                                   cccGGGCGAGCCGGGGGAAAGGTTACCCCGATTCTCTTTGCCCGACGCTCGATAACGATTCGAGATTCGGGCCCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGac                                                                                                                <  1:11260/251‑1 (MQ=255)                                                                                                                                          gAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAAcc                                                                                                          >  1:85111/1‑250 (MQ=255)                                                                                                                                          gAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACAAAAcc                                                                                                          >  2:58517/1‑250 (MQ=255)                                                                                                                                                 gtgAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGCTCGGACACAACCACAGTGAt                                                                                                  >  2:416748/1‑251 (MQ=255)                                                                                                                                                       ggTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCt                                                                                            >  1:302828/1‑251 (MQ=255)                                                                                                                                                           acacGGATGCTCTTAGCCGTACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTTCAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAACGATCGGACACAACCACAGTGATGTGGCTGGt                                                                                         <  2:606205/250‑1 (MQ=255)                                                                                                                                                                     tctTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACtt                                                                                                                                                        <  1:295597/177‑1 (MQ=255)                                                                                                                                                                     tctTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGGACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACtt                                                                                                                                                        >  2:295597/1‑177 (MQ=255)                                                                                                                                                                      ctTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAACTCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGGGCGCTTCAGGAt                                                                             >  1:474031/1‑251 (MQ=255)                                                                                                                                                                       ttAGCCGGACGCTCGATAACTATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGAt                                                                             <  2:626084/250‑1 (MQ=255)                                                                                                                                                                        tAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGAt                                                                             <  2:209806/249‑1 (MQ=255)                                                                                                                                                                            cGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGCGCGCTTCAGGATGCGATc                                                                       >  2:316541/1‑251 (MQ=255)                                                                                                                                                                             ggACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCt                                                                      <  2:593466/251‑1 (MQ=255)                                                                                                                                                                                cGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCa                                                                   <  1:78762/251‑1 (MQ=255)                                                                                                                                                                                 gCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCAc                                                                  <  2:319396/251‑1 (MQ=255)                                                                                                                                                                                     gATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGGTGCGATCTGCACGAcc                                                              >  2:89727/1‑251 (MQ=255)                                                                                                                                                                                      aTAACGATACGAGATACGGACGCTTTAGCCAATTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTTTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAATCCAGGTAGATTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTTCACGGCCt                                                             <  2:68699/251‑1 (MQ=255)                                                                                                                                                                                         aCGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCtttt                                                          >  1:322503/1‑251 (MQ=255)                                                                                                                                                                                             tACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCAc                                                      >  1:37497/1‑251 (MQ=255)                                                                                                                                                                                              aCGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCAc                                                      <  2:433666/250‑1 (MQ=255)                                                                                                                                                                                                gagaTACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCAcgcg                                                   >  1:581835/1‑251 (MQ=255)                                                                                                                                                                                                   aTACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGc                                                 >  1:373317/1‑250 (MQ=255)                                                                                                                                                                                                     aCGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCAt                                               <  2:298705/250‑1 (MQ=255)                                                                                                                                                                                                       ggACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAAt                                            <  2:98761/251‑1 (MQ=255)                                                                                                                                                                                                            cTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCCCGACCTTTTGCACGCGGCATAATGCGCt                                       >  2:105677/1‑251 (MQ=255)                                                                                                                                                                                                                           tAGTCAGGTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCt                         <  1:370590/250‑1 (MQ=255)                                                                                                                                                                                                                                 ggTACTGACGTACTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTa                  <  2:42178/251‑1 (MQ=255)                                                                                                                                                                                                                                            aCTTTAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTc       >  2:496933/1‑251 (MQ=255)                                                                                                                                                                                                                                              tttAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGACTACCATTCGGATGTACTTTCTGACCCATTGCTCGCCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGCTGTGGCTGGCGCGCTTCAGGAGTCGATCTGCCCGACCTTTTGCACGCGGCCTAATGCGCTTCATGCTCGGGCCTTCGTCTCCGTAAATTTTCGt     >  2:328079/1‑251 (MQ=255)                                                                                                                                                                                                                                              tttAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCtt                        <  2:569944/232‑1 (MQ=255)                                                                                                                                                                                                                                              tttAAAATCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCtt                        >  1:569944/1‑232 (MQ=255)                                                                                                                                                                                                                                                 aaaaTCACTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAAc  >  1:167036/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     AGCAACCAGTTTTGCGTCCAGTTCAGGCTTACGAACTTCGGCGATGTTGATCTGTGCAGGAACGCCAGCGATGTCCGCTACGACCTTACGCAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAAC  >  NC_000913/3449533‑3450023

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