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%	163.4 / NA	47	E115K (GAA→AAA)	fliC	flagellar filament structural protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/24);  total (23/24)

TGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGA  >  NC_000913/2003018‑2003507                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    tgtttgtGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTc                                                                                                                                                                                                                                                  <  1:513536/250‑1 (MQ=255)  gtttgtGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCt                                                                                                                                                                                                                                                >  1:167917/1‑251 (MQ=255)        ggtggtAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATa                                                                                                                                                                                                                                           <  1:510497/250‑1 (MQ=255)                cACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGAc                                                                                                                                                                                                                                     >  2:334584/1‑248 (MQ=255)                cACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGa                                                                                                                                                                                                                                  <  1:130709/251‑1 (MQ=255)                    aaaaGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGTGATTTAATTTTGTCCTGGATAGAAGACAGATCAg                                                                                                                                                                                                                              <  1:206021/251‑1 (MQ=255)                                 ggAGCGCTAGTGGAAACAATCTCGTTATTTTTAACACTCAAACCATAAAGGACGCGAGTTTTAGCATCAATATGCTACAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGt                                                                                                                                                                                                                 >  2:296922/1‑251 (MQ=255)                                        tAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGt                                                                                                                                                                                                           >  2:653999/1‑250 (MQ=255)                                            ggTNACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTgg                                                                                                                                                                                                      <  2:601839/251‑1 (MQ=255)                                                    tATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATCGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTcc                                                                                                                                                                                                                                                 >  1:77036/1‑200 (MQ=255)                                                    tATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTcc                                                                                                                                                                                                                                                 <  2:77036/200‑1 (MQ=255)                                                     aTCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGATTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTAcc                                                                                                                                                                                              <  2:72578/250‑1 (MQ=255)                                                     aTCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTAcc                                                                                                                                                                                              <  2:220162/250‑1 (MQ=255)                                                                      tAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGc                                                                                                                                                                            >  1:314551/1‑251 (MQ=255)                                                                               cAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATag                                                                                                                                                                                                                                          >  1:401072/1‑180 (MQ=255)                                                                               cAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTGGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATag                                                                                                                                                                                                                                          <  2:401072/180‑1 (MQ=255)                                                                                                         tGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCAAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTACCGCACACGCTGTAAGTTGATGTTGATTTCGGACAGCGCGCCTTc                                                                                                                                         >  2:195191/1‑251 (MQ=255)                                                                                                          gCTTCAGGTCGATATTGTTAGTCTGGTGATCATTGGCGCCATCCGGGATTTTCAGGGAGCCTTTTTTTGCCAGCCCGTTCACGCCGTTTAATTGGGTCTGACCAGATACGCTGCGAATTTCATCCAGACTGGATTTAATTTTTTCCTGGATAGAAGAAAGATCAGACTCAGCGTTAGTACCGGTAGTGGCCTGTTCCGTCAGTTCGCGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCTCGCCTTCg                                                                                                                                        <  2:376748/251‑1 (MQ=255)                                                                                                                    gatagtgatagtCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTgcg                                                                                                                              >  1:611164/1‑251 (MQ=255)                                                                                                                    gatagtgatagtCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTgcg                                                                                                                              <  1:297507/251‑1 (MQ=255)                                                                                                                                                ccAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCgg                                                                                                                                       <  1:605914/214‑1 (MQ=255)                                                                                                                                                ccAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCgg                                                                                                                                       >  2:605914/1‑214 (MQ=255)                                                                                                                                                 cAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTAACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACgg                                                                                                 >  1:192117/1‑251 (MQ=255)                                                                                                                                                    ccTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGgcc                                                                                              <  2:377815/251‑1 (MQ=255)                                                                                                                                                       ggATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGgccg                                                                                             <  1:553593/249‑1 (MQ=255)                                                                                                                                                         aTTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTg                                                                                         >  1:141303/1‑251 (MQ=255)                                                                                                                                                                 ggAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAgg                                                                                  <  2:464156/250‑1 (MQ=255)                                                                                                                                                                 ggAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGc                                                                                 >  2:260653/1‑251 (MQ=255)                                                                                                                                                                   aGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCt                                                                               <  1:334584/251‑1 (MQ=255)                                                                                                                                                                    gCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCtt                                                                              <  1:613320/251‑1 (MQ=255)                                                                                                                                                                         tttttGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTACATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATg                                                                         >  2:646455/1‑251 (MQ=255)                                                                                                                                                                                        tCACGACGGTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTTTCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGTCCTTTAATGTTAGATGTGAAACg                                                           <  2:139883/250‑1 (MQ=255)                                                                                                                                                                                                 tGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAAt                                                  <  2:431670/250‑1 (MQ=255)                                                                                                                                                                                                   aaCTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAgg                                                                                  <  1:445911/216‑1 (MQ=255)                                                                                                                                                                                                   aaCTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAgg                                                                                  >  2:445911/1‑216 (MQ=255)                                                                                                                                                                                                         gTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCATTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTACTGTTAGAGGTGAAACGGTTAGCAATCGCCTGAcc                                         >  2:158796/1‑251 (MQ=255)                                                                                                                                                                                                           cTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGGCTGCGCAACGGAGATACCGTCGTTGGCGGTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCg                                       >  2:406240/1‑251 (MQ=255)                                                                                                                                                                                                                  gATACGCGGTCAATTTCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTc                                <  1:442806/251‑1 (MQ=255)                                                                                                                                                                                                                                ttCATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGt                  >  1:314157/1‑251 (MQ=255)                                                                                                                                                                                                                                  cATCCAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATAAGACTCAGAGCTAGTACCGGTAGTGGCCTGTAACGTCAATTCAACCACACGCTGTAAGATGTTGTTGATTTCGGACAGCGCGCTTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCACCCTAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAAACGCCTGACAAGCTGCGTCCCCCTTCGCGCTNtt                 >  2:593411/1‑250 (MQ=255)                                                                                                                                                                                                                                     ccAGACGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTGATa             <  1:653999/251‑1 (MQ=255)                                                                                                                                                                                                                                         aCGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATcc                            <  2:205045/232‑1 (MQ=255)                                                                                                                                                                                                                                         aCGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATcc                            >  1:205045/1‑232 (MQ=255)                                                                                                                                                                                                                                          cGGGATTTAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCAAACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACNCa         >  2:606869/1‑250 (MQ=255)                                                                                                                                                                                                                                               tttAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTc                                                                                     <  1:556220/169‑1 (MQ=255)                                                                                                                                                                                                                                               tttAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTc                                                                                     >  2:556220/1‑169 (MQ=255)                                                                                                                                                                                                                                                 tAATTTTGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCaga  >  1:38301/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    TGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTTCATGGAGCCATTTTTTGCCAGCACGTTCACGCCGTTGAACTGGGTCTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGA  >  NC_000913/2003018‑2003507

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