Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,503,689	A→G	I266T (ATT→ACT)	fryA ←	putative PTS enzyme: Hpr, enzyme I and IIA components

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,503,689	0	A	G	100.0%	230.4 / NA	59	I266T (ATT→ACT)	fryA	putative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/29);  total (30/29)

CAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAATGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCT  >  NC_000913/2503544‑2503836                                                                                                                                                  |                                                                                                                                                    cAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTc                                                                                                                                                <  1:365876/151‑1 (MQ=255) cAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTc                                                                                                                                                <  1:65776/151‑1 (MQ=255)    gCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGt                                                                                                                                             <  1:350731/151‑1 (MQ=255)      cACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGttt                                                                                                                                           >  1:281932/1‑151 (MQ=255)          cAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCTGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCa                                                                                                                                       <  2:680917/151‑1 (MQ=255)           aGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCa                                                                                                                                       <  2:319037/150‑1 (MQ=255)             tACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGc                                                                                                                                    >  1:624900/1‑151 (MQ=255)                        ccGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCttgtt                                                                                                                         >  2:176766/1‑151 (MQ=255)                          gAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCttgtt                                                                                                                         >  2:638648/1‑149 (MQ=255)                          gAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTgc                                                                                                                       >  2:230976/1‑151 (MQ=255)                           aGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTgcg                                                                                                                      >  2:445888/1‑151 (MQ=255)                             gCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGcgcg                                                                                                                    <  1:411395/151‑1 (MQ=255)                              cACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCgg                                                                                                                   <  2:554241/151‑1 (MQ=255)                                  ggCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTc                                                                                                               <  1:217070/151‑1 (MQ=255)                                     cAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGt                                                                                                            >  1:113468/1‑151 (MQ=255)                                        aaTCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTcc                                                                                                         <  2:280360/151‑1 (MQ=255)                                                      aGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATg                                                                                           >  1:277959/1‑151 (MQ=255)                                                      aGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATg                                                                                           <  1:557516/151‑1 (MQ=255)                                                      aGGTGCGCCCGGTTGTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATg                                                                                           <  2:597139/151‑1 (MQ=255)                                                            gcCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATgcagca                                                                                      >  2:746175/1‑150 (MQ=255)                                                             cccGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAAc                                                                                    <  2:612576/151‑1 (MQ=255)                                                               cGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTg                                                                                  <  2:632668/151‑1 (MQ=255)                                                                   tttCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTc                                                                               <  2:294713/150‑1 (MQ=255)                                                                    ttCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGc                                                                             >  1:410895/1‑151 (MQ=255)                                                                     tctcCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGc                                                                             <  1:121333/150‑1 (MQ=255)                                                                      ctcCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTg                                                                           >  2:51526/1‑151 (MQ=255)                                                                        ccAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTg                                                                           <  1:386099/149‑1 (MQ=255)                                                                        ccAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGat                                                                         >  2:353574/1‑151 (MQ=255)                                                                          aGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGGTATGCAGCAACTGTTCGCTGatat                                                                       <  2:26473/151‑1 (MQ=255)                                                                          aGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGatat                                                                       >  2:346837/1‑151 (MQ=255)                                                                          aGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGCGATATGCAGCAACTGTTCGCTGatat                                                                       >  2:28275/1‑151 (MQ=255)                                                                                aTGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGaa                                                                 <  1:80458/151‑1 (MQ=255)                                                                                  gCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAAt                                                                <  1:561644/150‑1 (MQ=255)                                                                                  gCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATg                                                               >  2:306854/1‑151 (MQ=255)                                                                                  gCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATg                                                               <  1:351122/151‑1 (MQ=255)                                                                                  gCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATg                                                               >  1:480279/1‑151 (MQ=255)                                                                                    ccGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTc                                                             >  1:477980/1‑151 (MQ=255)                                                                                        aaGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCt                                                          >  1:754092/1‑150 (MQ=255)                                                                                           aTTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAaca                                                      <  1:248590/151‑1 (MQ=255)                                                                                           aTTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAac                                                       >  2:283035/1‑150 (MQ=255)                                                                                             ttttCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACg                                                    >  2:551101/1‑151 (MQ=255)                                                                                               ttCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGt                                                   <  1:758680/150‑1 (MQ=255)                                                                                                      cAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAg                                           <  1:562812/151‑1 (MQ=255)                                                                                                       aaGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCgcgc                                             >  1:188715/1‑148 (MQ=255)                                                                                                           cTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAAt                                      >  1:310179/1‑151 (MQ=255)                                                                                                             cAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAAt                                      <  2:283122/149‑1 (MQ=255)                                                                                                             cAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCg                                    <  1:502117/151‑1 (MQ=255)                                                                                                             cAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCg                                    >  2:388665/1‑151 (MQ=255)                                                                                                              aaaaaCTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGc                                   >  2:347022/1‑151 (MQ=255)                                                                                                               aaaaCTGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCt                                  >  2:381433/1‑151 (MQ=255)                                                                                                                    tGGCTTGGGGTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTggc                               >  2:501500/1‑149 (MQ=255)                                                                                                                          ggggTTAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTggcagaggcag                       <  2:187516/151‑1 (MQ=255)                                                                                                                               tAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAg                   <  1:176773/150‑1 (MQ=255)                                                                                                                               tAAATCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGt                  >  1:118589/1‑151 (MQ=255)                                                                                                                                aaaTCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGtt                 <  2:666675/151‑1 (MQ=255)                                                                                                                                aaaTCTTCAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGtt                 >  1:119536/1‑151 (MQ=255)                                                                                                                                       cAGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGgcgc           >  1:354723/1‑150 (MQ=255)                                                                                                                                        aGCCACCAGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCaa         <  1:48293/151‑1 (MQ=255)                                                                                                                                           caccaGAGTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAAcc      >  2:294694/1‑151 (MQ=255)                                                                                                                                               agagTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGc   >  2:39511/1‑150 (MQ=255)                                                                                                                                               agagTGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCt  <  1:143822/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CAAGCCACTCAGTACCGGGATCGCCGAGGCACGGGCCAGAATCAGTGTATGCGAGGTGCGCCCGGTTTTCTCCAGAATCATGCCCGCAAGATTTTTCAAATCAAGGCTCAAAAACTGGCTTGGGGTTAAATCTTCAGCCACCAGAATGGTCGGTTTTTCAAGCACCAGCTTGTTGCGCGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCT  >  NC_000913/2503544‑2503836

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