Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,003,264	C→T	E115K (GAA→AAA)	fliC ←	flagellar filament structural protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,003,264	0	C	T	100.0%	126.4 / NA	38	E115K (GAA→AAA)	fliC	flagellar filament structural protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/17);  total (21/17)

TGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAG  >  NC_000913/2003018‑2003506                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   tgtttgtGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTc                                                                                                                                                                                                                                                 <  2:148862/250‑1 (MQ=255)                   cAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCCCGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCa                                                                                                                                                                                                                              <  1:310914/251‑1 (MQ=255)                              aCTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGtt                                                                                                                                                                                                                   <  1:221465/251‑1 (MQ=255)                                            ggTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAAGCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTGTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTgg                                                                                                                                                                                                     <  1:28789/251‑1 (MQ=255)                                               aaCTGTATCGTTAGTTTTAACGCTAAAACCAGCAAGGCCAAGTGTTTTAGCTTAAATCTGCGCTAGATCGATAGTGATAGTCTGGTTATCATTTGCGCGATCCTGGATTTGCATGGAGCCATTGTTTGCAAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATGTCATCCAGACGGGATTTAATTTTGTCCTTGATAGAAGCCAGATCAGATTCAGAGTTATTTCCGGTAGGGGCCt                                                                                                                                                                                                  <  2:413928/251‑1 (MQ=255)                                                   ttATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTa                                                                                                                                                                                               <  2:217949/249‑1 (MQ=255)                                                   gTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTTCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCTTTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTTTAc                                                                                                                                                                                              <  2:189294/251‑1 (MQ=255)                                                       cGTTATTTTTAAGGCTAAAACCATCAAGGCCAAGAGTTTTCGCATCACTCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTCCACCAACCTGGATATTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCAAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGCGGCCTGTACCgac                                                                                                                                                                                          >  2:45137/1‑249 (MQ=255)                                                               ttAACGCTCACACCATCAAGGCCAAGAGTTTTAGCATCGATCTGCTTCAGATCGCTAGTGATAGTCTGGTTATCATTTGCGCCACCCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCAAGACGGGATTTAATTTTGTCCTGGATAGAAGACCGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTGAGTTAACg                                                                                                                                                                                  >  2:399680/1‑251 (MQ=255)                                                                tAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGc                                                                                                                                                                                 >  2:10002/1‑251 (MQ=255)                                                                                          gTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTc                                                                                                                                                       <  2:299522/251‑1 (MQ=255)                                                                                                      aTCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGcc                                                                                                                                           >  2:96303/1‑251 (MQ=255)                                                                                                                 aTCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCt                                                                                                                                >  1:193783/1‑251 (MQ=255)                                                                                                                           atagtCTTGTTATCATGTGCGCCACCCTTGATTTTCATGGTGCCATTTTTTGCCCGAACTTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTTAACTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCAGTAGTGGCCTCTCCCATCATTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACGGCGCGCCTTCGGTGGTCTGCGCTACGga                                                                                                                      <  2:145729/251‑1 (MQ=255)                                                                                                                             agtCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACgg                                                                                                                       >  2:48840/1‑248 (MQ=255)                                                                                                                                  tgtttccATTTGCGCCCCCCTGTATTTTCATGGAGCCATTTTTTCCCAGCCCGTTCCCGCCGTTGAACTGTGTCTGACCAGATACGCAGTCAATTTCATCCAGACTGTATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGTCTTTACCGTCAGTTCTCTCACACGCTGTCAGTTGTTGTTGACTTCGGCCAGCGCGCCTTCGGTGGTCTGCGCAACCGCGATAcc                                                                                                                <  2:1650/244‑1 (MQ=255)                                                                                                                                  ggTTATCATTTGCGCCAACCTGGATTTTCATGGCGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACcg                                                                                                               >  1:342131/1‑251 (MQ=255)                                                                                                                                              cgcCAACCTGGATTTTCATGTAGCCATTTGTTGCTAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCATACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCTCGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCTTTa                                                                                                   <  2:199657/251‑1 (MQ=255)                                                                                                                                                    ccTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGgcc                                                                                             >  1:410523/1‑251 (MQ=255)                                                                                                                                                           tttCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCCGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCTTCGTTGGCGTTACGGGCCGCCTg                                                                                        >  2:206200/1‑249 (MQ=255)                                                                                                                                                           tttCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTg                                                                                        >  1:160338/1‑249 (MQ=255)                                                                                                                                                            ttNATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAg                                                                                      <  2:382588/250‑1 (MQ=255)                                                                                                                                                                            ttGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGCCAGCGCGCCTTCGGTGGTCTGCGCAACGGGGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGtt                                                                      >  2:138715/1‑250 (MQ=255)                                                                                                                                                                             tGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGtt                                                                      >  2:9399/1‑249 (MQ=255)                                                                                                                                                                               ccAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTag                                                                    >  1:12115/1‑249 (MQ=255)                                                                                                                                                                                    aCGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGGTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCCGCGTCCGGCCTTTAATGTTAGAGGTGa                                                              >  2:218798/1‑250 (MQ=255)                                                                                                                                                                                    aCGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCAAGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGa                                                              >  1:43965/1‑250 (MQ=255)                                                                                                                                                                                      gTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAAc                                                           >  1:277257/1‑251 (MQ=255)                                                                                                                                                                                              cGTTGAACTGGGTCTGCCCAGATACGCGGTCAATTTACTCCAGAAGGGATTTAATTTTGTCCTGGATCGAAGACAGCTAAGACTCAGAGTTAGTACCGGTAGTGGCCTGTAACGTCAATTAACACACAAGCTGAAAGTTGTTTTTGATTTCGGACAGCGCGCCATCGGTGGTCTGCGCAACGGAGGTACCGTCGTTGGCGATACGGGCCGCCTGAGTCAGGACTTTAATGTTAGAGGTGAAACGGTTAg                                                     >  2:227964/1‑249 (MQ=255)                                                                                                                                                                                              cGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATGTTGTCTTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCGGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGAGTCCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAAAGTTAGAGGTGAAACGGTTAg                                                     <  1:227964/249‑1 (MQ=255)                                                                                                                                                                                                             gACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCt                                    <  2:131181/251‑1 (MQ=255)                                                                                                                                                                                                             gACCAGATACGCGGTCAATTTCACCCAGACGGGATTTAATGGTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTGCCGGTAGTGGTCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCTCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTGGGCGTTACGGGCCGCCTGAGTCGGGCCATTAATGTTAGAGGTGAAACGATTAGCAACCGCCTGACCCGcg                                    <  2:219883/251‑2 (MQ=255)                                                                                                                                                                                                                     aCGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATc                            >  2:402405/1‑251 (MQ=255)                                                                                                                                                                                                                     aCGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATc                            >  1:330860/1‑251 (MQ=255)                                                                                                                                                                                                                            cAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTcgcg                     <  2:43965/251‑1 (MQ=255)                                                                                                                                                                                                                                          cGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTAGCCTGTACCGTCAGTTCACGCACACGCTGTACGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGCGATACCGTCGTTGGCGTTACGGGCCGCCTGCGTCAg                                                                                  >  1:378820/1‑176 (MQ=255)                                                                                                                                                                                                                                              aTTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTAACGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCCACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGGAAGGCCTTTAATGTTAGAGGAGAAACGGGTAGCAATCGCCTGACCCGCTGCGTCATCCTTcgagccggtaataaggcagcc    >  2:310955/1‑231 (MQ=255)                                                                                                                                                                                                                                                ttAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCag  <  2:342131/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   TGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAG  >  NC_000913/2003018‑2003506

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