Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	793,050	C→T	E247K (GAA→AAA)	modF ←	ABC family protein ModF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	793,050	0	C	T	100.0%	236.1 / NA	69	E247K (GAA→AAA)	modF	ABC family protein ModF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/33);  total (36/33)

CCGCTGCCGCGACGTCGTCCGAAAAGCGTCAAATCGTTGCTGTAACCTTGCGGATGATCGCCAGTAACCAGGCTTAATAACGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTCATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCC  >  NC_000913/792806‑793296                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       ccgcTGCCGCGACGTCGTCCGAAAAGCGTCAAATCGTTGCTGTAACCTTGCGGATGATCGCCAGTAACCAGGCTTAATAACGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTAtccgg                                                                                                                                                                                                                                                  >  1:520633/1‑251 (MQ=255)                   ccGAAAAGCGTCAAATCGTTGCTGTAACCTTGCGGATGATCGCCAGTAACCAGGCTTAATAACGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACAcc                                                                                                                                                                                                                                <  1:198323/251‑1 (MQ=255)                    ctaaaaGCGTCAAATCGTTGCTGTAACCTTGCGGATGATCGCCAGTAACCAGGCTTAATAACGTCGCTTTTCCTGCCCCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTCAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTACGGCTCCACGTGCTGCGGGTTTATCCCGCTCCGTCAGTTGCTCAccc                                                                                                                                                                                                                               >  2:15978/3‑250 (MQ=255)                     gAAAAGCGTAAAATCGTTGCTTTAACCTTGCGGATGATCGCCAGTAACCAGGCTTAATAACGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCtt                                                                                                                                                                                                                              <  2:19417/251‑1 (MQ=255)                        aaGCGTCAAATCGTTGCTGTAACCTTGCGGATGATCGCCAGTAACCAGGCTTAATAACGTCGATTTTCCTGCACCAGTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCaa                                                                                                                                                                                                                           <  1:473540/251‑1 (MQ=255)                          gCGTCAAATCGTTGCTGTAACCTTGCGGATGATCGCCAGTAACCAGGCTTAATAACGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGc                                                                                                                                                                                                                         >  1:115656/1‑251 (MQ=255)                                aaTCGTTGCTGTAACCTTGCGGATGATCGCCAGTAACCAGGCTTAATAACGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCa                                                                                                                                                                                                                   <  2:222127/251‑1 (MQ=255)                                   cGTTGCTGTAACCTTGCGGATGATCGCCAGTAACCAGGCTTAATAACGTCGATTTTCCCGCACCAGTTGGGCCGACAATTTGCCCGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGGTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTa                                                                                                                                                                                                                <  2:133711/251‑1 (MQ=255)                                      tGCTGTAACCTTGCGGATGATCGCCAGTAACCAGGCTTAATAACGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATg                                                                                                                                                                                                              <  1:495814/250‑1 (MQ=255)                                          gTAACCTTGCGGATGGTCGCCAGTAACCAGGCTTAATAACGGCTATTGTCCTGCCCCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTGCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGGTTAGCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCa                                                                                                                                                                                                         <  1:265642/251‑1 (MQ=255)                                          gTAACCTTGCGGATGATCGCCAGTAACCAGGCTTAATAACGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGATCGTTGGCGGTTAAGGCGTGACGTGCTGATGGTTTATACGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCAATATGCGccc                                                                                                                                                                                                         >  1:129400/1‑250 (MQ=255)                                                    ggATGATCGCCAGTAACCAGGCTTAATAACGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGcgac                                                                                                                                                                                               >  1:231705/1‑251 (MQ=255)                                                            gCCAGTAACCAGGCTTAATAACGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTTTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCtt                                                                                                                                                                                       <  2:96157/251‑1 (MQ=255)                                                             ccAGTAACCAGGCTTAATAACGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCttg                                                                                                                                                                                      >  1:531500/1‑251 (MQ=255)                                                                           taataaCGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTAGCCGGCTCCTGCAGTTGCACACCTTCAAGCTGTTCACGATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTcc                                                                                                                                                                         <  2:36772/250‑1 (MQ=255)                                                                           taataaCGTCGATTTTACTGCACCATTTGGCCCGACAATTTGCCAGGGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAATCCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCGTCCAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCGGTTcc                                                                                                                                                                         <  2:418732/250‑1 (MQ=255)                                                                             ataaCGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCctct                                                                                                                                                                      <  1:338429/251‑1 (MQ=255)                                                                               taCGTCGATTTTCCTGCACCATTTGGCCCGACATTTTGAAAGTGTTCGCCTGGATTCACCTGCCAGCTCAGGTTCTTAAGAATGGGGCGATCGTTATAAGCAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCAGGCAGTTGCACACCTTCACGCTGTTCACTAGTCGCCAGTTCCGCTACGAGTGCTTGTTTGAGCAGTTCCTCttt                                                                                                                                                                    >  2:53508/2‑251 (MQ=255)                                                                                     gATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGcc                                                                                                                                                              <  1:137556/251‑1 (MQ=255)                                                                                      aTTTTCCTGCACCATGTGGCCCGACAATGTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATGGTTCAGCATAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCa                                                                                                                                                             <  2:235175/251‑1 (MQ=255)                                                                                        tttCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTAACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGt                                                                                                                                                           >  1:194592/1‑251 (MQ=255)                                                                                                         cccGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAAt                                                                                                                                          <  2:49792/251‑1 (MQ=255)                                                                                                            gACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCCTTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGGTGCCAACCTTCAAGCTGTTCACTATGCGCCAGTTGCCCGACGAGGGCTTGTTGGAGAAGTTCCTCTTTCGCGCCATTTTCCGCTAACGTGCAATccc                                                                                                                                       >  2:437370/1‑250 (MQ=255)                                                                                                            gACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTAtccggctc                                                                                                                                                                                                                                               >  1:315368/1‑147 (MQ=255)                                                                                                            gACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCcg                                                                                                                                       <  2:92801/251‑1 (MQ=255)                                                                                                            gACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTAGAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTAttcggctc                                                                                                                                                                                                                                               <  2:315368/147‑1 (MQ=255)                                                                                                                  ttGCCAGTGATCGCCTGGATTCACCCGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGAAGTGCTGAAGGTTTATCCGGCTCCCGCAGTTGCACACCTTCCAGCTGTTCACTATTCGCCAGTTGCGCGGCGGGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCCGTTTCCGCTAACGTGCCATCcgccagcc                                                                                                                                 >  2:487236/1‑250 (MQ=255)                                                                                                                       aGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCcgccagcacgcca                                                                                                                            <  1:42841/251‑1 (MQ=255)                                                                                                                       aGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCcgccagcacgcc                                                                                                                             >  1:477587/1‑250 (MQ=255)                                                                                                                                cTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGCAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCCGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGCCGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTgg                                                                                                                   >  2:214938/1‑251 (MQ=255)                                                                                                                                   gATTCACCTGCCAGCTAAGGTTATTAAGAATGGGACGATCGCTATAAGAAACCACGCCATTGTTCCGCACAATGCGCGGTTCGTTGGCGGGTACGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGAc                                                                                                                 >  1:232030/1‑250 (MQ=255)                                                                                                                                         ccTGCTAGCTAAGGTTATTAAGAATGGGGCGATCGTTAGAAGAAACCACGCCATTGTTCAGCGCAATGCGCGTTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTGTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGATGAGTGCGTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAATTGGACAAACTcc                                                                                                          <  2:172839/251‑1 (MQ=255)                                                                                                                                         ccTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACGATGCGCCAGTTGCGCGACGAGTGCGTGGTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTcc                                                                                                          <  2:228499/251‑1 (MQ=255)                                                                                                                                                   aaGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCCCGCCATTGTTCCGCACAATGCGCGGTTAGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGCGCAGTTACTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCAAGCAAACTGGACAAACTCCGGGATCTCAt                                                                                                >  2:272832/1‑251 (MQ=255)                                                                                                                                                         attaAGAATGGGGCGATAGTTATAAGAAACAACGCCATTGTTCAGTAAAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGATACGGCAGTTGCACAACTTAAAGGTGTTCAATATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCAAGTAAACTGGACAAACTCAGGGATCTCAACGAAgg                                                                                          >  2:213990/1‑250 (MQ=255)                                                                                                                                                           gaaGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAAGGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTTCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACGATGCGCCGGTGGCGCCACGAGTGCTTGTTGGAGCATTTCCTCTTTAGCGCCAGTTTCCGCTAACGTTCAAACCGCTAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCg                                                                                         <  1:516437/249‑1 (MQ=255)                                                                                                                                                             aGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTGGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGTATCTCATCGAAGCGAtt                                                                                      <  2:64674/251‑1 (MQ=255)                                                                                                                                                             aGAATAGGGCGATCGTTATCAGAAACCACGCCATTGTTCAGCACACTGCGCGGATCGTTGGTGGGTAAAGCGAGACGTGCTGAAGGTTTATCCGGCACCGGAAGGTGAACACCTTCACGCTGTTTCCTATGCGCCAGTTGCGCGACGGGTGCTTGTTGGAGCAGTTCCTCTTTAGCGACCGCTTCCGCTAACGAGCCATCCGCCAACAAGCACGCAAACTGGACCACCTCCAGGAtatccacgaaccgatt                                                                                      >  2:26515/1‑236 (MQ=255)                                                                                                                                                               aaTGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGACGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGCGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTg                                                                                     >  2:221716/1‑250 (MQ=255)                                                                                                                                                                   gggCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCCCACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCAc                                                                                >  1:248190/1‑251 (MQ=255)                                                                                                                                                                    ggcgttctttATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTGTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCAcc                                                                               <  2:392288/243‑1 (MQ=255)                                                                                                                                                                                aaGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTa                                                                   <  2:8185/251‑1 (MQ=255)                                                                                                                                                                                 aGAAACCACGCCATTGTTCAGCAAAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGGGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCTCTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTGTCCGCTAACGTGCAATCCGCGAGCTCGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGTTATTGATCACCAGTACCAGAGTaa                                                                  <  2:265340/251‑1 (MQ=255)                                                                                                                                                                                  gAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTTACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAAt                                                                 >  1:532496/1‑251 (MQ=255)                                                                                                                                                                                   aaaCCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGCAAGGCGTGACGTGATGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAAt                                                                 >  1:527499/1‑250 (MQ=255)                                                                                                                                                                                     cccACGCCATTGTTCAGCACAAGGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGTACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATAcc                                                              <  2:520633/250‑1 (MQ=255)                                                                                                                                                                                          gCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACt                                                         >  1:486128/1‑251 (MQ=255)                                                                                                                                                                                           ccATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTg                                                        <  1:468829/251‑1 (MQ=255)                                                                                                                                                                                                   cAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCCTCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACg                                                >  1:214811/1‑251 (MQ=255)                                                                                                                                                                                                    aGCACAAGGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGCCGGTTTATGCGGCTGCGGCAGTTGCACCTCTTCACGCTGTTCGCTATGCGGCAGTTGCGCGAAGAGTGATTGTTGGAGCCGTTCCTCTTTAGCGCTAGTTTCCCCTAACGTGCCATCCGCCAGCACGCCAGCAAACTGGGCAAACTCAGGGCTCTCATCGGAGCGATTGAGCACCAGTACCAGCGTAATACCGGACGGATGTCACg                                                >  1:52708/1‑250 (MQ=255)                                                                                                                                                                                                      cacaATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCa                                                                                                                                                             >  1:225293/1‑139 (MQ=255)                                                                                                                                                                                                      cacaATGCGCGGTTCGTTGGCGGGGAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGTAGCAGTTCCGCTTTAGCGCCa                                                                                                                                                             <  2:225293/139‑1 (MQ=255)                                                                                                                                                                                                         aaTGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGACGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGAAATACCGGACTGATGTAACGAGGCg                                           >  1:195618/1‑250 (MQ=255)                                                                                                                                                                                                              gcgGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGc                                     >  1:466745/1‑251 (MQ=255)                                                                                                                                                                                                               cgGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCt                                    >  1:79264/1‑251 (MQ=255)                                                                                                                                                                                                               cgGTTCGTTGGCGGGTAAGGCGCGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCCAACGTGCACTCCGCCAGCACGCCAGCCAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAACACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGACGcc                                    >  2:201816/1‑250 (MQ=255)                                                                                                                                                                                                                 gTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTca                                  <  1:520559/251‑1 (MQ=255)                                                                                                                                                                                                                  ttCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTc                                   >  1:491344/1‑249 (MQ=255)                                                                                                                                                                                                                       tGGCGGGTAAGGCGAGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTcagccagc                            >  1:359516/1‑251 (MQ=255)                                                                                                                                                                                                                          cGGGTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGTGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCACACTGGACAAACTCCGGGATCTGATCGAAGCGATTGAGCACCAGGACCAGAGTAATACCGGGCTGATGGAAAGAAGCGAGGCGCTCCGCGAgct                           >  2:162631/1‑249 (MQ=255)                                                                                                                                                                                                                             gTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTACTACCGGACTGATGTa                                                   >  2:430053/1‑222 (MQ=255)                                                                                                                                                                                                                             gTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTa                                                   <  1:430053/222‑1 (MQ=255)                                                                                                                                                                                                                             gTAAGGCGTGACGTGCTGAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGa                      <  2:477587/251‑1 (MQ=255)                                                                                                                                                                                                                                      gACGTGCTGAAGGTTTATCCGGCTCCGGCAGTGGCACACCTTAAAGCTGTGCACTACGCGCCAGGTTCGCGACGAGTGCTTGTTGGAGCAGTTCCTCttt                                                                                                                                                                    <  1:265638/100‑1 (MQ=255)                                                                                                                                                                                                                                              gAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATcc       >  2:12976/1‑249 (MQ=255)                                                                                                                                                                                                                                              gAAGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAg     <  1:214938/251‑1 (MQ=255)                                                                                                                                                                                                                                                aGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTca                                  >  1:46458/1‑220 (MQ=255)                                                                                                                                                                                                                                                aGGTTTATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCGTGTGGGAGCAGTTCCTCTGTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCGCAGCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGGCGCTca                                  <  2:46458/220‑1 (MQ=255)                                                                                                                                                                                                                                                 ggTTTATCCGGCTCCGGCAGTTGCACACCTTCACGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGcc  >  1:348766/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       CCGCTGCCGCGACGTCGTCCGAAAAGCGTCAAATCGTTGCTGTAACCTTGCGGATGATCGCCAGTAACCAGGCTTAATAACGTCGATTTTCCTGCACCATTTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACGTGCTGAAGGTTCATCCGGCTCCGGCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCC  >  NC_000913/792806‑793296

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