Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	968,252	G→A	K544K (AAG→AAA)	msbA →	ATP‑binding lipopolysaccharide transport protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	968,252	0	G	A	100.0%	154.6 / NA	48	K544K (AAG→AAA)	msbA	ATP‑binding lipopolysaccharide transport protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/29);  total (19/29)

CAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTG  >  NC_000913/968009‑968495                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    caATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGAc                                                                                                                                                                                                                                               <  1:477480/250‑1 (MQ=255)       gATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCg                                                                                                                                                                                                                                        >  1:7620/1‑251 (MQ=255)              aaCGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGgtcgt                                                                                                                                                                                                                                 >  2:238727/1‑251 (MQ=255)               aCGGTCGCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTTGCGGTCAGCGTCAGCGTATTTCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCATACGAAATCGTGgtcggc                                                                                                                                                                                                                                <  2:446195/251‑3 (MQ=255)                  gtcgcGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCGGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAg                                                                                                                                                                                                                             <  2:193829/247‑1 (MQ=255)                      cgagaCAGTGATTGGTGAAAACGGCGTGCGGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGAt                                                                                                                                                                                                                          <  2:226825/246‑1 (MQ=255)                          aCAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGgtgt                                                                                                                                                                                                                     <  1:238727/251‑1 (MQ=255)                           cAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGCTTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTAGTCGAGGTTGGTGTc                                                                                                                                                                                                                    >  2:124642/1‑251 (MQ=255)                               gtttGGTGTAAACGCTGTGCTGCTCCTTGGCGGTCAGCGTCCTCGTAGTGCTATCGCTCGGGCCTTGTTGCGTGATAGCCTGATTCTGATTCTGTACGAAGCTCCCGCGTCTTTGGATACCGAATCCGAACGTTCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCACTCTGGTGATTGCCCACCGCGTGTCTTCCATTGAAAAAGCAGACGAAACCGTGGTCGTCGAGGATGGGGTCAtt                                                                                                                                                                                                                 <  2:137127/248‑1 (MQ=255)                                        aaaCGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCCTTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCCTTGTGGAAcgcg                                                                                                                                                                                                       >  1:165366/1‑251 (MQ=255)                                               gTGCTGCTCTCTGGCGGGCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGACCGTGCGGTTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACAAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGcc                                                                                                                                                                                                >  2:303205/1‑250 (MQ=255)                                                         ctGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTg                                                                                                                                                                                      >  2:313844/1‑251 (MQ=255)                                                             cGGTCCGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCGACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTg                                                                                                                                                                                  <  2:456300/251‑1 (MQ=255)                                                             cGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTg                                                                                                                                                                                  <  1:193995/251‑1 (MQ=255)                                                                  agcgtcatcgtGTTGCTATCGCTCTAGCCGTGTTGCGCGATAGCCCGATTCTGATTCTGGGCGAAGCTACCTCGGCTCTGGATTCCGAATCCGAACGTGCTATTCAGGCGGCACTGGATGTGTTGCAGAAACACCGTACCTCTCTGGTGGTTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCAc                                                                                                                                                                             <  1:304227/251‑1 (MQ=255)                                                                  agcgtcagcgtATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCAc                                                                                                                                                                             >  2:342403/1‑251 (MQ=255)                                                                                   aTCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGGCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACgcgc                                                                                                                                                             <  2:82779/250‑1 (MQ=255)                                                                                      gcaccaGCCTTGTTGCGTGATAGCCCGTGACTGATTCTGTACGAAGCTGCCGAGGCTCTCGATACCGAAACCGAACGTGCGAGTCTGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCGGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATTGTGTCATGGTGGAACGCGGTACGCAGAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAAc                                                                                                                                                          <  2:6860/245‑1 (MQ=255)                                                                                       cTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTgg                                                                                                                                                                                                             >  2:78958/1‑198 (MQ=255)                                                                                       cTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTgg                                                                                                                                                                                                             <  1:78958/198‑1 (MQ=255)                                                                                             ccTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTGCCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACaaa                                                                                                                                                  <  1:217885/251‑1 (MQ=255)                                                                                                         aTAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGc                                                                                                                                      <  2:99196/251‑1 (MQ=255)                                                                                                         aTAGCCCGATGCTGATTCTGGACGCAGCTACCTCGGCTCTGGATACCGAATCCTAACGTGCGATTTAGGCGGCACTGGATGAGTTGCAGAAAAACCGTCCCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCTTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATTACGATTTGCTTGGGCACCGCGGCGTTTACGCGCAACATCACAAAATGCAGTTTgg                                                                                                                                       <  1:112716/250‑1 (MQ=255)                                                                                                              ccGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATg                                                                                                                                 <  2:137668/251‑1 (MQ=255)                                                                                                                  ttcggattctgGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTAAGGCGGCACTGGATGAGTTGCAGAAAAACAGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGGCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTgg                                                                                                                                       <  1:368315/237‑1 (MQ=255)                                                                                                                  atcggattctgGACGAAGCTACCTCGGCTCTGTATACCGAATCCGAACGGGCGATTCAGGCTGCACTGGATGAGTTGCAGACAAACCGTACCTCTCTGGTGATTGCCCACCGCTTTTCTACCATTGAAAAAGCCGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGACCGCGTTAAGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCCCAACATGCATTTTggaaatctgtcg                                                                                                                             >  2:368315/5‑241 (MQ=255)                                                                                                                          ctgGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAAtc                                                                                                                     >  1:258293/1‑251 (MQ=255)                                                                                                                                               ctGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCTCCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCtt                                                                                                  <  1:119367/249‑1 (MQ=255)                                                                                                                                                      ccgaatccgaatGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTggcggc                                                                                         <  1:303205/251‑1 (MQ=255)                                                                                                                                                       cGAATCCGAACGTGCGATTCAGGAGGCACTGGATGAGGTGCAGAAACACCGTACCTCTCTGGTGATTGCCCACCGATTGTCTACCATTGGAAAAGCAGACGAAATCGTGGTCGTCGCGGATGGTGTCATTGTGGAACGAGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTCAGCGCAACTTCCCAAAATGCTGTTTTACCAATGATCGAAAAAAtctggtaaggtgatacacctttttggcgcgt                                                                                        >  2:206315/1‑226 (MQ=255)                                                                                                                                                             cGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATtgct                                                                                  <  1:342403/251‑1 (MQ=255)                                                                                                                                                             cGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGAGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTAtttct                                                                                  >  2:350300/1‑248 (MQ=255)                                                                                                                                                                 cGTGCGCTTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCtttttcttccc                                                                              >  2:60111/1‑241 (MQ=255)                                                                                                                                                                          cAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTgg                                                                     <  1:350300/251‑1 (MQ=255)                                                                                                                                                                                gcACTGGATGAGTTGCAGAAAAACCGTTCCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTAt                                                               <  1:211146/251‑1 (MQ=255)                                                                                                                                                                                            ttGCAGAAAAACCGTACCCCTCTGGTGATTGCCCACCGCTTTTCTACCATTGAAAAAGCAGACTAAACCGTGGTCGTCGAGTATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGTTGAATCCCCTTTGTGGCGTCTATTGCTGCCACTCTCCTGGTTGTATTGCCCGGttag                                                    <  2:211864/250‑4 (MQ=255)                                                                                                                                                                                              gCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTg                                                  <  1:482827/250‑1 (MQ=255)                                                                                                                                                                                                aGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGgcg                                               <  1:60111/251‑1 (MQ=255)                                                                                                                                                                                                 aaaaaaacgttccgctgtGGTGCTTGCCCAGCGCTTGTCTATCATTGAAAAAGCAGACTAAATCGTGGGCGTCTATGATGGTGTCATTGTGGACCGCGGTACGCATAACGCTTTGCGTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCAGGTGAATCCGCTTTGTGGCGGGTATTGCTGCCACTCTCCTGGTTGTATTGCCTGTTGAGTGgcgt                                              <  2:4141/234‑2 (MQ=255)                                                                                                                                                                                                            ccTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGGGGTCGGCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCttt                                   <  2:112115/251‑1 (MQ=255)                                                                                                                                                                                                             ctctctGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGCGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTTCTGCCACTAGCCTTGTTGTTTGGCCTgggggggggcggagtccggctttt                                  >  2:11525/1‑229 (MQ=255)                                                                                                                                                                                                             ctctctGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGCGTCATTGAGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGGGAGGGGCGCGATCCGACTTTg                                  >  1:463914/1‑251 (MQ=255)                                                                                                                                                                                                              tctctGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGATTGGCGCGATCCGTCTTTGc                                 >  1:430524/1‑251 (MQ=255)                                                                                                                                                                                                              tctctGGTGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGGTGTATGGCCTGGTGAGTGGGGCGATCCGGCTTTGc                                 >  1:107706/1‑251 (MQ=255)                                                                                                                                                                                                                     tGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAAt                                                                                                                                                >  1:315737/1‑133 (MQ=255)                                                                                                                                                                                                                     tGATTGCCCACCGCTTGTCTACCATTGAAAAAGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAAt                                                                                                                                                <  2:315737/133‑1 (MQ=255)                                                                                                                                                                                                                                   ttGTCTACCATTGAAAAAGCAGACGAAATCGTGGACGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGACCCTCTCCTGGTTGTATGGCCCGGGGAGTGGCCGGGTCCGCCTTTTTTTTAAACTAAAACTAAAgcgc            >  2:426464/1‑251 (MQ=255)                                                                                                                                                                                                                                              tGAAAAAGCAGACGAAACCGGGGGCGTCGAGGATGGTGTCGTTGTGGAACGCGGTACGCATAACGATTTGCTGGAGCCCCGCGGCGTTTACGCGCAACTTCACAAAATGCATTTTGGCCAATGATCGAAAAAATCTGGGCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCTCTCCTGGCGTg  <  2:395851/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    CAATAAGATGGATAACGGTCTCGATACAGTGATTGGTGAAAACGGCGTGCTGCTCTCTGGCGGTCAGCGTCAGCGTATTGCTATCGCTCGAGCCTTGTTGCGTGATAGCCCGATTCTGATTCTGGACGAAGCTACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTG  >  NC_000913/968009‑968495

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: ‑