Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	4,235,408	Δ1 bp	intergenic (‑2175/‑499)	lysC ← / → yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	4,235,407	0	A	.	100.0%	274.1 / NA	60	intergenic (‑2174/‑500)	lysC/yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (24/36);  total (24/36)

CAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATA  >  NC_000913_3_bme_pgi/4235260‑4235542                                                                                                                                                    |                                                                                                                                        cAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGat                                                                                                                                        <  1:290536/149‑2 (MQ=255)      ccTTAGATTAAAGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcg                                                                                                                                  >  2:8464/1‑149 (MQ=255)        ttCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgtc                                                                                                                                >  1:436171/1‑149 (MQ=255)                    gggttGGAACTGGGTAAGCAACGGGCGGCGGAACGTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGTTCGCCACGCTCCCACGCAAGGGCGTGTCGCGCATCTGCATGTGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGcc                                                                                                                    <  2:429496/145‑1 (MQ=255)                    ggCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGcc                                                                                                                    <  1:8464/149‑1 (MQ=255)                        tGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGAt                                                                                                                >  1:671535/1‑149 (MQ=255)                          gaaaggggTAAGGACCGGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAGGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTTCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATaa                                                                                                              <  2:153836/144‑1 (MQ=255)                           aaCTGGGTAAGGCACTGGCGACGGAACTTTTGCCGGGTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                             <  1:657327/149‑1 (MQ=255)                           aaCTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                             <  2:731232/149‑1 (MQ=255)                                 gTAAGGAACTGGCGACGGAACTTTTGCCGGTTGGATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTc                                                                                                       <  1:648761/149‑1 (MQ=255)                                      gAACTGGCGACGGAACTTTTGCCGGTTGTCTCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                  <  1:696737/149‑1 (MQ=255)                                           ggggACGGAACTTTTGCCGGTTGGATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                             <  2:521683/146‑1 (MQ=255)                                             cggcGGAACTTTTGCCGGTTGTATGTGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                           <  2:26617/146‑1 (MQ=255)                                              gACGGAACTTTTGGCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                          <  1:291124/149‑1 (MQ=255)                                                cGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACCAGGCCGCTCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                        >  1:158101/1‑149 (MQ=38)                                                   acagtttGCGGGTTGGATCGGGAAAGGAAGGGGTGAGTGTTCGCGACGCGTCCACGCAAGGGCTTGTCGCGCATCCGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                                                                     <  2:84331/145‑1 (MQ=38)                                                        ttGCTGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGAGACGCTTCCACGCAAGGGCATGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ACAACCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                <  1:27533/149‑1 (MQ=38)                                                        ttGCCGGTTGTATCGGGCAGGGAAGGGGCGAGGGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                <  1:71799/149‑1 (MQ=255)                                                           ccGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCCGCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCCGCCGGCAACCGATGCCTg                                                                             >  1:476392/1‑149 (MQ=255)                                                             ggTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCACCGATGCCTGAt                                                                           >  1:167872/1‑149 (MQ=255)                                                              gTTGTATCGGCCAAGGAAGGGGCGATTGGTCTCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCACGCGCGAAGGAAAGCTTG‑ACCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                                                          <  2:669862/149‑1 (MQ=38)                                                               ttGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                           >  2:572780/1‑131 (MQ=255)                                                               ttGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                           <  1:572780/131‑1 (MQ=255)                                                               ttGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGc                                                                         >  2:66677/1‑149 (MQ=255)                                                                   aTCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcg                                                                     >  1:152075/1‑149 (MQ=255)                                                                    tCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                    >  2:757428/1‑149 (MQ=255)                                                                       ggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                                                                 >  2:426458/1‑149 (MQ=255)                                                                        gCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                >  1:131423/1‑149 (MQ=255)                                                                        gCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTc                                                                >  1:506721/1‑149 (MQ=255)                                                                            ggaaggGGCGAGTGGTTGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGt                                                            <  2:184210/149‑1 (MQ=255)                                                                            ggaaggGGCGAGTCGTAGCGGCGTTGCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAACGCTTG‑ACCATCGTCGATCTGTAGGCCGCATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCCGATGCCACGCGGTCGCGt                                                            <  2:417108/149‑1 (MQ=37)                                                                            ggaaggGGCGACTGGTCGCGACGCTTTCACGCAAGGGCTTGGCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATACGGCGTTCACGCCGCATCCGGCAACCGAGTCCTGATGCGACGCGGTCGCGt                                                            <  2:176890/149‑1 (MQ=38)                                                                                ggggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                        >  2:647896/1‑149 (MQ=255)                                                                                 gggCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt                                                       <  2:225634/149‑1 (MQ=255)                                                                                   gCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                    >  1:133016/1‑134 (MQ=255)                                                                                   gCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                    <  2:133016/134‑1 (MQ=255)                                                                                    cgggTGGTCGCGACGCTTCCACGAAAGGGCTTGGCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                    <  2:68380/146‑1 (MQ=255)                                                                                    cgggTGGGCGCGACGGTTCCAAGGAAGGGGTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                    <  2:232949/146‑1 (MQ=37)                                                                                    cGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                    <  1:647896/149‑1 (MQ=255)                                                                                     gAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg                                                   <  1:757428/149‑1 (MQ=255)                                                                                        tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCt                                                >  1:413287/1‑149 (MQ=255)                                                                                           tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCGAGCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGACTACa                                             >  1:556610/1‑149 (MQ=38)                                                                                            cgcgACGCTTCCACGCAAGGGCTCGTCGCGCATCTGCATGCGCGACGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGAATAAGGCGCTCACGCGGCACCCGGCTACCGATGCCTGATGTGACGCGGTCGCGTCTTATCAGGCCTACAg                                            <  1:683888/149‑1 (MQ=37)                                                                                             gcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTCGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGCTGCCTGATGCGACGCGGTCGCGTCTTATCAGGCATACAgg                                           >  1:571973/1‑149 (MQ=255)                                                                                                    ttCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGAt                                                                                                                >  1:675566/1‑73 (MQ=255)                                                                                                    ttCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGAt                                                                                                                <  2:675566/73‑1 (MQ=255)                                                                                                      ccACCCAAGGGCTTGTCGCGCACCTGCATGCGCGAAGGAAAGCGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCg                                  <  1:710150/149‑1 (MQ=255)                                                                                                       taggaaaGGGCTTGTCGGGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                 <  2:78387/144‑1 (MQ=255)                                                                                                        aCGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                >  1:126907/1‑149 (MQ=255)                                                                                                           cAAGGGCTGGTCGCGCATCTGCACGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGTCTACAGGTCGATGCCGATATg                             <  2:278813/149‑1 (MQ=255)                                                                                                              gggCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCAGCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTAc                          <  1:94061/149‑1 (MQ=255)                                                                                                                  ttGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCg                      >  2:302655/1‑149 (MQ=255)                                                                                                                       gcgcGTCGGGATGCGGGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                 <  2:588923/149‑1 (MQ=25)                                                                                                                        cgcATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCg                <  1:426458/149‑1 (MQ=255)                                                                                                                                taTGCGCGAAGGAAAGCTTG‑ATCACCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa        <  1:488023/148‑1 (MQ=255)                                                                                                                                   gcgcGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTGCAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATAc     <  2:506721/149‑1 (MQ=255)                                                                                                                                     gcgAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACAt   >  1:327667/1‑149 (MQ=255)                                                                                                                                      ggAAGGAAAGCTTG‑ATCATCGTAAATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGAAACGCAATGGTGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa  <  1:560085/148‑1 (MQ=25)                                                                                                                                      cgAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa  >  2:323326/1‑149 (MQ=255)                                                                                                                                        aaGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                  <  1:659175/67‑1 (MQ=255)                                                                                                                                        aaGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                  >  2:659175/1‑67 (MQ=255)                                                                                                                                                    |                                                                                                                                        CAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATA  >  NC_000913_3_bme_pgi/4235260‑4235542

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