Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,013,246	G→A	E416K (GAG→AAG)	pqiA →	intermembrane transport protein PqiA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,013,246	0	G	A	100.0%	145.0 / NA	42	E416K (GAG→AAG)	pqiA	intermembrane transport protein PqiA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/21);  total (21/21)

AAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGGAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAACCACCATTAAAAGCCGTAGCGTTGACGTCGGCGTGGTTGAAAGCGCCAC  >  NC_000913/1013003‑1013485                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                aaGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTGATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTg                                                                                                                                                                                                                                          <  1:21177/251‑1 (MQ=255)    cgcgACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGAc                                                                                                                                                                                                                                        >  2:6595/1‑250 (MQ=255)        aCAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGGGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTAt                                                                                                                                                                                                                                   >  1:51545/1‑251 (MQ=255)              aaaGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATc                                                                                                                                                                                                                             <  1:387767/251‑1 (MQ=255)               aaGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCt                                                                                                                                                                                                                            <  1:444430/251‑1 (MQ=255)                      cATTTGATTTATGCAGTTGTTGAGCTTGTCGGCCGCTGGCCGATGATTGACGCTTTCGTTATCGCGGTGCTCTCGGAGCTGGTCCGTATGGGTGGTTTAATGAGTATTTATCCGGCAATGGATGCATTCATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACTTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATg                                                                                                                                                                                                                     >  2:27173/1‑251 (MQ=255)                           gATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGaa                                                                                                                                                                                                                <  1:173085/251‑1 (MQ=255)                           gATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGa                                                                                                                                                                                                                 <  1:283305/250‑1 (MQ=255)                               tATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGcc                                                                                                                                                                                                             <  1:396740/250‑1 (MQ=255)                                                tGTAGGCCGCTGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTga                                                                                                                                                                                             >  1:219321/1‑249 (MQ=255)                                                 gTAGGCCGCTGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTgaaga                                                                                                                                                                                          <  1:118058/251‑1 (MQ=255)                                                   aGGCCGCTGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAAc                                                                                                                                                                                        <  1:208387/251‑1 (MQ=255)                                                          tGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGtctc                                                                                                                                                                                 >  1:383020/1‑251 (MQ=255)                                                               gatgatTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCgtg                                                                                                                                                                            >  1:2997/1‑251 (MQ=255)                                                                          ttttCGTTATCGCGGTGCTCTCGGCGCTGGTGCGAATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAACTCCCGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTcc                                                                                                                                                                 >  1:255381/1‑251 (MQ=255)                                                                             tCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTa                                                                                                                                                               <  1:216942/250‑1 (MQ=255)                                                                              cGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGCATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTAt                                                                                                                                                              <  2:46690/250‑1 (MQ=255)                                                                                  aTCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTc                                                                                                                                                          <  1:6595/250‑1 (MQ=255)                                                                                    cgcgGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACg                                                                                                                                                       >  1:117453/1‑251 (MQ=255)                                                                                       ggTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGgcg                                                                                                                                                    <  2:219321/251‑1 (MQ=255)                                                                                                                ggAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCtttttt                                                                                                                           <  1:86219/251‑1 (MQ=255)                                                                                                                      ttAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCAtt                                                                                                                     >  1:233148/1‑251 (MQ=255)                                                                                                                               aTTNATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATc                                                                                                            <  2:411525/251‑1 (MQ=255)                                                                                                                                           aTGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGa                                                                                                  >  1:16288/1‑249 (MQ=255)                                                                                                                                             gggTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTg                                                                                              <  2:216971/251‑1 (MQ=255)                                                                                                                                                         gTTTGCTTTAGTCTTCATAATGACAATGTTTTCTGCTATGACTTTGGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTCTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCAcc                                                                                  <  2:255381/251‑1 (MQ=255)                                                                                                                                                                      gtcATAATGACAATAGTATCTGCTAAGACGTTTGATCCGCGTTTGTCGAGGAATCGACAACCTGAATGAGAGCATGTGAAGTACTGACCGTATGGAATATAATGATGGGGAAGCCAAAATCCAGAAAGTGACGAACTGGTCTCCCGTGTGGATATTTCCTATCGTAACGGCTCTAGTTGGGGCCTGGGTTCTTTTTTATCATTACATCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATACgga                                                                       >  2:183292/1‑249 (MQ=255)                                                                                                                                                                        cATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCggaagg                                                                    >  1:443201/1‑250 (MQ=255)                                                                                                                                                                               aCAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGa                                                              >  1:231204/1‑249 (MQ=255)                                                                                                                                                                                  aTGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTg                                                         <  2:383020/251‑1 (MQ=255)                                                                                                                                                                                          tGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGAGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGCAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAacc                                                 >  2:305803/1‑251 (MQ=255)                                                                                                                                                                                              aTGACGTTTGCCCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCCTTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAaccaca                                              >  2:300429/1‑249 (MQ=255)                                                                                                                                                                                                gACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGTGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAACCACCAtt                                           >  2:43864/1‑251 (MQ=255)                                                                                                                                                                                                     ttGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAACCACCATTaaa                                        >  1:98927/1‑249 (MQ=255)                                                                                                                                                                                                       gACCCGCTTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAACCACCATTAAAAGcc                                    <  1:300429/251‑1 (MQ=255)                                                                                                                                                                                                         cccGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTGTGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAACCACCATTAAAAGCCGt                                  >  2:174190/1‑251 (MQ=255)                                                                                                                                                                                                          ccGCGTTTGTCGTGGGATCGTCATCCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATGACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAACCACCATTAAAAGCCGTa                                 <  1:43864/251‑1 (MQ=255)                                                                                                                                                                                                             cgTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAACCACCATTAAAAGCCGTAGCg                              <  1:251577/251‑1 (MQ=255)                                                                                                                                                                                                                                aaCCTGAATCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAACCACCATTAAAAGCCGTAGCGTTGACGTCGGCGTGGTTg            >  1:239486/1‑250 (MQ=255)                                                                                                                                                                                                                                aaCCTGAATCAGAGCATGAGAAGACCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAACCAACATTAAAAGCCGTAGCGGTGACGTCGGCGTGGNTg            >  2:417214/1‑250 (MQ=255)                                                                                                                                                                                                                                      aaTCAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAACCACCATTAAAAGCCGTAGCGTTGACGTCGGCGTGGTTGAAAgcgc     <  2:350400/251‑1 (MQ=255)                                                                                                                                                                                                                                         cAGAGCATGAGAAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAACCACCATTAAAAGCCGTAGCGTTGACGTCGGCGTGGTTGAAAGCGCcac  >  2:229260/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                AAGCGCGACAGTGAAAGAATGCATTTGATTTATGAAGTTGTTGAGTTTGTAGGCCGCTGGTCGATGATTGACGTTTTCGTTATCGCGGTGCTCTCGGCGCTGGTGCGTATGGGAGGTTTAATGAGTATTTATCCGGCAATGGGTGCATTAATGTTTGCTTTAGTCGTCATAATGACAATGTTTTCTGCTATGACGTTTGACCCGCGTTTGTCGTGGGATCGTCAACCTGAATCAGAGCATGAGGAGTCCTGACAGTATGGAATCTAATAATGGGGAAGCCAAAATCCAGAAAGTGAAGAACTGGTCTCCCGTGTGGATATTTCCTATCGTCACGGCGCTCATTGGGGCCTGGGTTCTTTTTTATCATTACAGCCATCAGGGACCGGAAGTGACCCTGATCACCGCGAATGCGGAAGGAATTGAAGGTGGCAAAACCACCATTAAAAGCCGTAGCGTTGACGTCGGCGTGGTTGAAAGCGCCAC  >  NC_000913/1013003‑1013485

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