Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,181,995	G→A	L279L (CTG→TTG)	potD ←	spermidine preferential ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,181,995	0	G	A	100.0%	131.6 / NA	45	L279L (CTG→TTG)	potD	spermidine preferential ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/16);  total (29/16)

CGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAGCAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCAGTTTTTTCAGCTCGT  >  NC_000913/1181752‑1182240                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      cGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCTCCGCAAGGTTTGGCGTTGGATAACCGTTAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGa                                                                                                                                                                                                                                                  >  2:237900/1‑249 (MQ=255)            accaccGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAAtt                                                                                                                                                                                                                                      <  1:565951/250‑1 (MQ=255)                    ccgcTTGCTGAATTAACGTACTGCTTACAGCTTCTGCTAATACTCTTCATAAATACTGCTGGCTGCGCCAACGTCATTCTGACATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGTTAAACGATAGTTTCAGCAAACTGTT‑TTGCCACTACCGGGCGAAACAGGATGGTGATCACTTTCAgcgcgc                                                                                                                                                                                                                             >  2:101307/1‑251 (MQ=255)                      gcTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCct                                                                                                                                                                                                                           <  1:333749/251‑1 (MQ=255)                        ttGCTTAATTAACGGCCTGCTTTCAGCTTCTGATAATACTCTGCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTCCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCcta                                                                                                                                                                                                                          <  2:233446/250‑2 (MQ=255)                          gCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGTACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCcttctt                                                                                                                                                                                                                       <  1:390700/251‑1 (MQ=255)                                ttAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGtttt                                                                                                                                                                                                                 >  1:451145/1‑251 (MQ=255)                                 tAACGTCCTGCTTTCAGCTTCTGATAATTCTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGCGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGtttt                                                                                                                                                                                                                 <  1:120728/250‑1 (MQ=255)                                  aaCGTCCTGCTTTCATCTTCTGATCATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCCTTCTGCCCTTCGCCCTTTTTAATGGTTTCAGCGTCCGGGTAGAGCGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGAAAGGATTGGCGTTGCATAACCGATAGTTTCCGCAACCTGTT‑TTGCCCCATCCGGGCGCAACAGGAAGTTTATCAAATTCAGCGCGCCTTCTTTGATTTTg                                                                                                                                                                                                               >  2:322625/1‑251 (MQ=255)                                                 gCTTCTGATAATACTCTTCATAACTGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGTTAGTTTCAGAAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATAAATTTCAGCGCGCCTTCTTTGTTTTTGTCAtttacagggttc                                                                                                                                                                                                >  2:467168/1‑242 (MQ=255)                                                    tCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGcc                                                                                                                                                                                             >  1:657107/1‑251 (MQ=255)                                                                ctTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGtccatc                                                                                                                                                                                 >  1:426644/1‑251 (MQ=255)                                                                              ctggctgCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGATGTTGATCAATTTCAGCGCTCCTTCTTTGTTTTTGGCCTTTTCCGGGATCGCCATGCTGGCCATCCAGacagtgcggcc                                                                                                                                                                   >  2:5115/1‑241 (MQ=255)                                                                                     cgcCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTc                                                                                                                                                            >  1:352050/1‑251 (MQ=255)                                                                                                           gCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTGCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCg                                                                                                                                      <  1:467168/251‑1 (MQ=255)                                                                                                              aTTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTAGCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACC‑GTTGTTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTGTGTCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTGCGGCCACACCACGTCAATTGGCGTAc                                                                                                                                   <  1:617886/251‑1 (MQ=255)                                                                                                                   tAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCTCCGAACGGTTTAGAGTAGGATAACCGATAGTTTAAGCAACCTGTT‑TTGCCACATCGGAGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTTTCCTTTGCCGGCATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGcc                                                                                                                              >  1:431955/1‑251 (MQ=255)                                                                                                                        gTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCg                                                                                                                         >  1:276497/1‑251 (MQ=255)                                                                                                                           tCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGTCTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCCGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAAGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACCCCACGTCATTTGGCGTAACCGCCTGGCGTGc                                                                                                                      >  2:677376/1‑251 (MQ=255)                                                                                                                             aGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCa                                                                                                                     >  1:636356/1‑250 (MQ=255)                                                                                                                                aTCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCCATTGGCGTACCCGCCTGGCGTGCaccg                                                                                                                  >  1:40680/1‑247 (MQ=255)                                                                                                                                        ggagTGTTTTATCGTTCGCCACTTCTGGAGTTAACAGCTTACGCACCGAACGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCCCATCCGCGCGCAACAGGCCATTGATCAATTTCAGCGCGCCGTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGTACATGCCGCCTTCTTTCGGCCACACCACGTCAATTGTCGTACCCGCCTCACGTGCCACGacaacata                                                                                                          >  2:578842/2‑243 (MQ=255)                                                                                                                                        ggagTGTTTTATCGTTCGCCACTTCGGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCTGACCCGCCTGGCGTGCAACGAAAGCAGaa                                                                                                         <  1:542456/250‑1 (MQ=255)                                                                                                                                         gagTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGACAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGaa                                                                                                         >  2:542456/1‑250 (MQ=255)                                                                                                                                                 tATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGATTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACCCCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGACCCGTTcc                                                                                                  >  2:233341/1‑249 (MQ=255)                                                                                                                                                  aTCGTACGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGCTAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCTGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAAAGAAAGCAGAACCGATcccg                                                                                                >  2:325504/1‑248 (MQ=255)                                                                                                                                                                 tGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCAGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTa                                                                                 >  1:59687/1‑250 (MQ=255)                                                                                                                                                                       tAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCg                                                                           <  1:325504/250‑1 (MQ=255)                                                                                                                                                                           aGCTTACGCGCCGCAAGGTTGGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTc                                                                      <  2:138830/251‑1 (MQ=255)                                                                                                                                                                              ttACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCAttt                                                                                                                                                                                                         >  1:137902/1‑117 (MQ=255)                                                                                                                                                                              ttACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCAttt                                                                                                                                                                                                         <  2:137902/117‑1 (MQ=255)                                                                                                                                                                                      cgcAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGt                                                           >  2:522066/1‑251 (MQ=255)                                                                                                                                                                                          aGGTTTAGCGTTGGATAACCGATAGGTTCAGCAACCTGAT‑TTGCCACTTCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGCCCATCCCGAACATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCAGGCGTGCAACGAAAGCAGAACCGTTCCAGCTCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGc                                                       >  2:480059/1‑251 (MQ=255)                                                                                                                                                                                                       gATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAg                                          <  1:419214/251‑1 (MQ=255)                                                                                                                                                                                                           aCCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTa                                       >  1:218307/1‑250 (MQ=255)                                                                                                                                                                                                              gATAGTTTCAGCAACCAGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCCAGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAATCAGAACCGTTCCATATCATGACGAGGTTAACTTCGCCTTCCATGTACAGGTTCGCCGATTTATCGGAGTTAAAgg                                   >  2:232902/1‑251 (MQ=255)                                                                                                                                                                                                                  gttctaGCAACCTGTT‑TGGCCACAGACGGGCGTAACAGGAAGTGGATAAATTTCAGCGGGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGACTTCTTTCGGCCACACCACGTCAAGTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTGCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTg                                <  2:636356/245‑1 (MQ=255)                                                                                                                                                                                                                         cAACCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTACAGGCTGCCACGttt                        >  2:335408/1‑251 (MQ=255)                                                                                                                                                                                                                          aaCCTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTg                       >  1:98608/1‑251 (MQ=255)                                                                                                                                                                                                                             cTGTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGAAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTAGCGGAGTTAAAGGCTGCCCCGTTTGGc                     >  2:474797/1‑250 (MQ=255)                                                                                                                                                                                                                               gTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTCTCGGAGTTAAAGGCTGCCACGTTTGGCATc                  >  2:6168/1‑251 (MQ=255)                                                                                                                                                                                                                               gTT‑TTGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGAACGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCAAGACAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGCTCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGGGTTACAGGCTGCCACGTTTGGCATc                  >  1:617849/1‑251 (MQ=255)                                                                                                                                                                                                                                 t‑ttGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCa                 <  1:6168/250‑1 (MQ=255)                                                                                                                                                                                                                                   ttGCCACATCCGGGCGCAACAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCAg                <  1:216946/250‑1 (MQ=255)                                                                                                                                                                                                                                                gcgcAACAGGAATTTGATCAATTGCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGAGCGCCAGGCGGTCCACCCAGATAATGCCGCCTTCTTTCGGCCACACCGCTGCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCAGGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCAGTTTTTTCAGCTCGt  <  1:297984/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      CGGGCGGTAATACCACCGCCCGCTTGCTGAATTAACGTCCTGCTTTCAGCTTCTGATAATACTCTTCATAAATGCTGCTGGCTGCGCCAACGTCATTCTGCCATTCGCCATTTTTAATGGTTTCAGCATCCGGGTAGAGTGTTTTATCGTTCGCCACTTCTGGACTTAACAGCTTACGCGCCGCAAGGTTTGGCGTTGGATAACCGATAGTTTCAGCAACCTGTT‑TTGCCACATCCGGGCGCAGCAGGAAGTTGATCAATTTCAGCGCGCCTTCTTTGTTTTTGGCATTTGCCGGGATCGCCAGGCTGTCCATCCAGAAAATGCCGCCTTCTTTCGGCCACACCACGTCAATTGGCGTACCCGCCTGGCGTGCAACGAAAGCAGAACCGTTCCAGATCATGCCGAGGTTAACTTCGCCTTCCATGTACGGGTTCGCCGGGTTATCGGAGTTAAAGGCTGCCACGTTTGGCATCAGTTTTTTCAGCTCGT  >  NC_000913/1181752‑1182240

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