Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,249,916	G→A	P78P (CCG→CCA)	yqjE →	inner membrane protein YqjE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,249,916	0	G	A	100.0%	221.9 / NA	70	P78P (CCG→CCA)	yqjE	inner membrane protein YqjE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/30);  total (40/30)

CGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCGCAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTC  >  NC_000913/3249670‑3250161                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      cgcgTCGATGATTGTGGCGGAATCTCACCCCGCACAACGGCCCGGCAAAATCTTTCTGGGCATCGGGCAGCGAATTGATTTTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGGTGGAAGACGAAAAAGCGAATTTCTTTCAACTTTTACTGAGGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGGCTTATGAGGCTGATGGTGCTAATTATTTGTGCGTTTGACCCACAat                                                                                                                                                                                                                                                   <  2:330989/251‑1 (MQ=255)  gcgTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAata                                                                                                                                                                                                                                                  <  1:2684/251‑1 (MQ=255)    gtcgtTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATc                                                                                                                                                                                                                                                >  1:93858/1‑251 (MQ=255)              aTGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATg                                                                                                                                                                                                                                       <  2:228031/250‑1 (MQ=255)              aTGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATg                                                                                                                                                                                                                                       <  1:221299/250‑1 (MQ=255)              aTGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGc                                                                                                                                                                                                                                      >  1:310163/1‑251 (MQ=255)                  cGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCgatg                                                                                                                                                                                                                                  <  1:641/251‑1 (MQ=255)                  cGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCgatg                                                                                                                                                                                                                                  >  1:191682/1‑251 (MQ=255)                          catcaCGCACAAGGGCCCGGTAAAAGCGTTCTGGGCATCGGGCCGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGGGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGCTTTccac                                                                                                                                                                                                                          >  2:95537/1‑251 (MQ=255)                          catcaCGCACAAGGGCCCGGTAAAAGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGcca                                                                                                                                                                                                                           <  1:301256/250‑1 (MQ=255)                           atcaCGCACAAGGGCCCGGTCAAAGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGCACTGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGCGCTGGAAGAGGAAAAAGCGCATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGACGGTGCTAATTATTTGGGCGGTTGACCTACAATACCGCCTGAATGCGCAGATTGccacc                                                                                                                                                                                                                         >  2:72184/1‑251 (MQ=255)                             caCGAACAAGGGCCCGGTAAAAGTGTTCTGGGCATCGGGCAGCGAATTGTTTCTGTCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGAGTGGCTGGAAGAGGAAAAAGCGAATCTCTTTTAACTTTTACTGATGCTGGGCCTGGCGATGCTTTTCGCTGCATTTGGTATTATGAGCCTGATGGTGCTAATTATTTGGGAGGTTGACCCACAATATCGCCTGAATGCAATCATTGCcaccac                                                                                                                                                                                                                       >  2:10096/1‑251 (MQ=255)                                           cGGTAAAAGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTtgctgc                                                                                                                                                                                                         <  1:212797/251‑1 (MQ=255)                                                aaaGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGAGTTGCTGCTACTg                                                                                                                                                                                                    >  1:97397/1‑251 (MQ=255)                                                  aGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTCCTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGGATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGc                                                                                                                                                                                                  >  2:81882/1‑251 (MQ=255)                                                   gCGTTCTGGGCTGCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCa                                                                                                                                                                                                 <  2:156929/251‑1 (MQ=255)                                                    cGTTCTGGGAATCGGGCAGCGAATTGTTTCTATCCTGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGGTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTACTTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCAc                                                                                                                                                                                                >  2:60087/1‑251 (MQ=255)                                                     gTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGAATTTGGTCTTATGCGCCTGATGGTGATAATTATTTGGGCGGTTGACCCACAATATAGCTTGAATGCGATGATTGCAACCGCCGTGGTGTTGCTGCTACTTGCACt                                                                                                                                                                                               >  2:1566/1‑251 (MQ=255)                                                            gCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGc                                                                                                                                                                                        <  1:64580/251‑1 (MQ=255)                                                            gCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGc                                                                                                                                                                                        <  2:95816/251‑1 (MQ=255)                                                                     aGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGACAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTgg                                                                                                                                                                               >  1:377154/1‑251 (MQ=255)                                                                      gCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGa                                                                                                                                                                              <  1:223123/251‑1 (MQ=255)                                                                      gCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGa                                                                                                                                                                              <  1:142492/251‑1 (MQ=255)                                                                       cGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGAc                                                                                                                                                                             >  1:106245/1‑251 (MQ=255)                                                                         aaTTGTTTCTATCATGGTTGAACTGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGATGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCATGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATATGGaagc                                                                                                                                                                           >  1:255590/1‑248 (MQ=255)                                                                            tgtctCTATCATGTTTGAAATGGGAGAGACACGTCTGCGGCTGGCGGTGTTGGAGCTGGAAGATGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGCCCCACAATAACGCCTGAATGCGATGATTGCCACCACCGTTGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTa                                                                                                                                                                         <  1:358776/246‑1 (MQ=255)                                                                                  tATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGCTTGGCGGTATCTGGCCGCTACGTAAAt                                                                                                                                                                  >  2:150562/1‑251 (MQ=255)                                                                                  tATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCAACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAAt                                                                                                                                                                  >  1:232523/1‑251 (MQ=255)                                                                                      aTGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAAATTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTCCGTAAATcgcg                                                                                                                                                              >  1:256690/1‑251 (MQ=255)                                                                                       tGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGCGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGt                                                                                                                                                             >  1:333042/1‑251 (MQ=255)                                                                                             aaaTGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCt                                                                                                                                                       >  1:407916/1‑251 (MQ=255)                                                                                             aaaTGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCt                                                                                                                                                       >  1:118305/1‑251 (MQ=255)                                                                                               aTGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTAc                                                                                                                                                     >  1:21464/1‑251 (MQ=255)                                                                                                  gTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTTGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGtt                                                                                                                                                  <  2:97397/251‑1 (MQ=255)                                                                                                   tAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGttt                                                                                                                                                 >  1:409172/1‑250 (MQ=255)                                                                                                            gTCTGCGACTGGCGGTGGTTGTGCTGGGAGCGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCCTTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGCATGCGATGCTTGCCACCACCGTGGTGTTGCTGCCACTAGCACTGATTGGCGGTATCTGGGCTCTACATAAATCCCGTAAGGCTCCGGTGCTGCGCCAt                                                                                                                                        >  2:271652/1‑251 (MQ=255)                                                                                                                      ggCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATg                                                                                                                              >  1:358864/1‑251 (MQ=255)                                                                                                                        gggtggtGGAGCTGGAAGAGGAAAAAGCGAAGCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGGAAAGCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAg                                                                                                                            <  2:153247/250‑1 (MQ=255)                                                                                                                           tggtggAGCTGGACGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGtt                                                                                                                          >  1:305010/1‑250 (MQ=255)                                                                                                                                gAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCaaa                                                                                                                    <  1:150562/251‑1 (MQ=255)                                                                                                                                   cTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACg                                                                                                                  >  1:317882/1‑250 (MQ=255)                                                                                                                                        agagGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTACTTATTTGGGCGGTTGACCCACAATATCGCCTGACTGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCggg                                                                                                            >  2:340597/1‑250 (MQ=255)                                                                                                                                         gagGAAAACGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACCCGCCCTGAGTTAGCAAACGATCGGc                                                                                                            >  2:202653/1‑250 (MQ=255)                                                                                                                                          cggAAAAAGCGCATCGCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGCTGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTACGTCTACGTTGCTGCGCCATACACGCCCTGAGTTAGCAACCGATCgcc                                                                                                            >  2:338100/2‑247 (MQ=255)                                                                                                                                                aaGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTc                                                                                                    >  1:245846/1‑251 (MQ=255)                                                                                                                                                 aGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCGGACGATGCTTTTCGCGGCATTTGGTCTTATGAGCCTGGTGGTGCTAATTATTTGGGCGGGTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACTTAAATCGCGTAAGTCTACGGTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCg                                                                                                   <  2:403893/251‑1 (MQ=255)                                                                                                                                                    gAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCg                                                                                                              >  1:87385/1‑237 (MQ=255)                                                                                                                                                    gAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCg                                                                                                              <  2:87385/237‑1 (MQ=255)                                                                                                                                                    gAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCgagg                                                                                                >  1:197979/1‑251 (MQ=255)                                                                                                                                                              cAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGCCCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGCGTCCCg                                                                                      >  1:398265/1‑251 (MQ=255)                                                                                                                                                                 cTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATAAACGCCCTGAGTTAGCAAAAGATCGGCAGCTGCTCGGGGAGGGGTCCCGTgg                                                                                   >  2:78331/1‑250 (MQ=255)                                                                                                                                                                     tACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGc                                                                                 >  1:137741/1‑249 (MQ=255)                                                                                                                                                                      aCTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCGTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAAGATCGCCTGAATGCGATGATGGCCACCACCGTGGTGGTGCTGCTACTGGCACTGAGTGGCGGTATCTGGACTCTACTTAAATCGCGTAAGTCTACGTGGCTGCGCCATACACGCCATGATTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAg                                                                               <  1:211253/250‑1 (MQ=255)                                                                                                                                                                         gATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTaa                                                                            <  1:95537/250‑1 (MQ=255)                                                                                                                                                                                             ttttCGCTGCAATTGGGATTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGCATGCGATGATTGCCACCACCGTGGTGTTGCTGCTCCTGGCACTGCTTGGCGGTATATGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGAGCCATACCCGCCATGAGTTAGCAAACGATCGGCAGCTGCCCGGGGAGGAGTCCCGTGAGCAGTAAAGTCGAACgggacaggagga                                                        >  2:414187/1‑239 (MQ=255)                                                                                                                                                                                                 cGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGGCGAACGTGAACGACGTAAGgc                                                    <  1:81882/250‑1 (MQ=255)                                                                                                                                                                                                           ggTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTa                                         <  1:67086/251‑1 (MQ=255)                                                                                                                                                                                                                 aTGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCTACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTAGCCAg                                    <  1:72184/250‑1 (MQ=255)                                                                                                                                                                                                                    aGCCTGATGGTGCTAATTTTTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCGGCTACTGGCACTGATTGGCGGTATCTGGACGCTACTTAAATCGCGTAAGGCTACGTTGCTGCGCTATAGACGCCATGAGTTAGCGAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTAGCCAGATcc                                <  1:178641/251‑1 (MQ=255)                                                                                                                                                                                                                     gCCTGATGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTACCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAAAGTGAACGACGTAAGGCGCAACTGCTTAGCCAGATcc                                >  1:147051/1‑250 (MQ=255)                                                                                                                                                                                                                          aTGGTGCTAATTATTTGGGCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTCCGTGGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGcaac                          <  2:224703/251‑1 (MQ=255)                                                                                                                                                                                                                                      attggggCGGTTGACCCACAATATCGCCTGAATGCGATGATTGCACCCGCCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGAt               <  2:147974/246‑1 (MQ=255)                                                                                                                                                                                                                                         tGGGCGGTTGACCCACATTATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCtt            <  1:305057/250‑1 (MQ=255)                                                                                                                                                                                                                                           ggCGGTGGGCCCACAATATCGCCTGAATGCGTTGTTTGCCACCGCCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGTAACAACGGCGGGATCTTTcc         <  2:256146/251‑1 (MQ=255)                                                                                                                                                                                                                                             cGGTTGTCCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTTCGTAAAGCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGGAAGGCGCAACTGCTTAGc                                       <  1:318393/219‑1 (MQ=255)                                                                                                                                                                                                                                             cGGTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAAGCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGGACGGGAACGACGGAAGGCGCAACTGATTAGc                                       >  2:318393/1‑219 (MQ=255)                                                                                                                                                                                                                                              ggTTGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGcc      >  1:288599/1‑251 (MQ=255)                                                                                                                                                                                                                                                 tGACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGGCGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGCGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAgt   >  1:165730/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gACCCACAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAgtc  >  1:322724/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gACCCACAATATCGCCAGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCTCTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGGAGCGGCTCGAGGAAGAGTCCCGTGAGCAGTAAAGTTGAACGTGAACAACGTAAGGCGCAACCGCTTAGCCAGATCCAGCAACAGCGGCTGGATCTTTCCGCCAgtc  <  1:38971/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      CGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGTTCTGGGCATCGGGCAGCGAATTGTTTCTATCATGGTTGAAATGGTAGAGACACGTCTGCGGCTGGCGGTGGTGGAGCTGGAAGAGGAAAAAGCGAATCTCTTTCAACTTTTACTGATGCTGGGCCTGACGATGCTTTTCGCTGCATTTGGTCTTATGAGCCTGATGGTGCTAATTATTTGGGCGGTTGACCCGCAATATCGCCTGAATGCGATGATTGCCACCACCGTGGTGTTGCTGCTACTGGCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTC  >  NC_000913/3249670‑3250161

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