Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	3,560,455	+G	intergenic (‑2/+1)	glpR ← / ← glpR	pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	3,560,455	1	.	G	98.5%	191.8 / ‑3.6	67	intergenic (‑2/+1)	glpR/glpR	pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (31/35);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTG  >  NC_000913_3_bme_pgi/3560307‑3560593                                                                                                                                                       |                                                                                                                                           gCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑c                                                                                                                                             >  2:925921/1‑149 (MQ=255)          gATAAACCGCGTCGAGCCATTGGGGATATGATCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGAGTGGTTGGCCTTGCCCTCGTGCCACGGCGTGTAACCCGCACTGGCAGGCAGCTCCGCACCGCCATGATGG‑CGCAGGAGCa                                                                                                                                    >  2:812582/1‑149 (MQ=255)          gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCa                                                                                                                                    >  1:1102066/1‑149 (MQ=255)          gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGACTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCa                                                                                                                                    >  2:634908/1‑149 (MQ=255)          gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTCCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGAGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCa                                                                                                                                    >  1:467576/1‑149 (MQ=255)           caaaaCAGCCCCGAGCCATTGGGGCTTGGCGCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGCCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAg                                                                                                                                   <  2:1476274/147‑1 (MQ=255)           aTAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAg                                                                                                                                   <  1:781481/149‑1 (MQ=255)           aTAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGACGCACCGCCATGATGG‑CGCAGGATCAg                                                                                                                                   >  1:333263/1‑149 (MQ=255)            tAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAgg                                                                                                                                  <  2:84809/149‑1 (MQ=255)             aaaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGt                                                                                                                                 <  1:687631/149‑1 (MQ=255)               aCAGCGTCGAGCCATTGGGGATTTGCTCCGACACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGACTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGttt                                                                                                                               >  1:1354314/1‑149 (MQ=255)                    gTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTGTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATGAGGTTTTGCTc                                                                                                                          >  2:163875/1‑149 (MQ=255)                     tCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTcc                                                                                                                         <  1:951168/149‑1 (MQ=255)                       gAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGc                                                                                                                       <  1:1116417/149‑1 (MQ=255)                          ccATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAg                                                                                                                    <  2:1214359/149‑1 (MQ=255)                            aTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGca                                                                                                                  >  1:189658/1‑148 (MQ=255)                             ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                 >  1:1196480/1‑149 (MQ=255)                                gggATTTGCTCCGCCACTTTGCGGGCGAGGCGCTCTTTTTCTTCGGTCTGGGTGCCCTTGCGTTCGTGCCACGGCGTGTTATCCGAACTGGAATGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                              <  1:1104807/149‑1 (MQ=255)                                 ggATGTGCTCCGCCTCTTTGCGGGCGATGCGCGCTTTTTCTTCGGTCTGGGTGGCCTTGCGAGCGTGCCACGGCGTGTTTACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                             <  2:1376600/149‑1 (MQ=255)                                   aTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                           <  1:1422221/149‑1 (MQ=255)                                              cACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                                                                <  2:880151/149‑1 (MQ=255)                                              cACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                                                                <  1:1023830/149‑1 (MQ=255)                                                  ttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt                                                                                            <  2:1063457/149‑1 (MQ=255)                                                   tGCGGGCGCAGCTCTCATTATATACGGTCTGGGCGGCCTTGCGATCGTCCCACGGCGTGCTCACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTc                                                                                           <  1:1370352/149‑1 (MQ=255)                                                        gCGATGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                      <  2:889481/149‑1 (MQ=255)                                                                 tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                                                                             >  1:194247/1‑149 (MQ=255)                                                                 tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                                                                             <  2:1102066/149‑1 (MQ=255)                                                                 tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTGATTGCGGTGGCGGCGAATAGTCTGCGGGcggacgga                                                                             >  2:699487/1‑149 (MQ=255)                                                                  ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                            >  2:978052/1‑149 (MQ=255)                                                                  ctTTTTCTTCGGGCTGGGGGGCCTTGCGAGCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCTCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                            <  2:1423137/149‑1 (MQ=255)                                                                      ttcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGTAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCCGCTCAGTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                        >  2:1243770/1‑149 (MQ=255)                                                                       tcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTTCCCGAACTGGAAGGCAGCGCCGCTCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                                       <  1:275227/149‑1 (MQ=255)                                                                       tcttcCGTCTGGGTGGCCTTGCGATCCTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGGCTGCGGGCTGACGGAGAAATg                                                                       <  1:820052/149‑1 (MQ=255)                                                                         ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAAATGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGTGAAATGct                                                                     >  2:938765/1‑149 (MQ=255)                                                                             gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                 <  1:925921/149‑1 (MQ=255)                                                                             gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGGTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                 <  1:615799/149‑1 (MQ=255)                                                                             gTCTGGGTGGCCATGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                 <  1:409582/149‑1 (MQ=255)                                                                              gcTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                <  1:446561/148‑1 (MQ=255)                                                                               cTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCGGCGGGCTGACGGAGAAATGCTCTcccc                                                               >  2:1247945/1‑148 (MQ=255)                                                                                  ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGct                                                            >  1:321442/1‑149 (MQ=255)                                                                                   gTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctc                                                           <  2:140102/149‑1 (MQ=255)                                                                                        cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTCCCAGCTCTTCgg                                                      >  1:269895/1‑149 (MQ=255)                                                                                        cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCAGTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAACTGCTCTACCAGCTCTTCgg                                                      >  2:178855/1‑149 (MQ=255)                                                                                        cTTGCGATCGTGCCACGGCGTGTGACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGAG‑GAGGATCAGGTTTTGCTCCGCCTGCTCATTGAGGTCGCGGCGGATAGTCTGGGGGCGGACGGAGAACTGCTCTACAAGCTCTTCgg                                                      >  2:1304394/1‑149 (MQ=255)                                                                                            cGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACt                                                  <  2:761735/149‑1 (MQ=255)                                                                                             gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAGAGTCTGCGGGCTGACGGAGAAATGCTCTACCCGCGCTGCGGGAcgg                                                 >  2:297727/1‑147 (MQ=255)                                                                                              aTCGTGCCACGGCGTGTTAACCGCACTGGAAGGCCGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTATTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCGCTACCAGCGCATCGGTACTGa                                                >  2:1008574/1‑149 (MQ=255)                                                                                              aTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTgc                                                >  1:720479/1‑148 (MQ=255)                                                                                                        ggCGGGTTAACCGATCTGGACGGCAGCTCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACTGATAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctg                                      <  1:812582/149‑1 (MQ=255)                                                                                                                aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTaa                              >  1:183414/1‑149 (MQ=255)                                                                                                                   cGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCa                           >  1:1256245/1‑149 (MQ=255)                                                                                                                     aaCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGt                         >  2:1067921/1‑149 (MQ=255)                                                                                                                     aaCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGt                         >  1:1192404/1‑149 (MQ=255)                                                                                                                      aCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                        <  2:194247/149‑1 (MQ=255)                                                                                                                        tGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGGTTTGCTCCGACAGCGCATTGAGGTCGCGCCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACAGACATAACCCTGCTGTGAACCCAGTTCg                      >  2:965052/1‑149 (MQ=255)                                                                                                                         ggAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGCAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                     >  2:1280377/1‑149 (MQ=255)                                                                                                                         ggAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                     <  2:321442/149‑1 (MQ=255)                                                                                                                         ggAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATAGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGa                     <  2:1196480/149‑1 (MQ=255)                                                                                                                          gAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCGTCGGGTCTGACATAACCCTGCTGTTTACCCAGTTCGat                    >  1:1269386/1‑149 (MQ=255)                                                                                                                          gAAGGCAGCGCAGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                    <  2:1354314/149‑1 (MQ=255)                                                                                                                              gCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATCCTCTACCAGCTCTTCGGTACTTACAAAACCCTGCTGATTAACCAGTACGataaaa                >  2:250704/1‑147 (MQ=255)                                                                                                                                   gccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt           <  2:333263/149‑1 (MQ=255)                                                                                                                                   gccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt           <  1:774754/149‑1 (MQ=255)                                                                                                                                   gccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGATCTTCGGTACTGACATAACCCTGCTGTGTAACCAGTTCGATAATACCGtt           <  1:978052/149‑1 (MQ=255)                                                                                                                                   gccgcACCGCCATGATGG‑CGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt           <  1:1225418/149‑1 (MQ=255)                                                                                                                                      gcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTTACCAGTTCGATAATACCGTtgtg        >  1:1404191/1‑149 (MQ=255)                                                                                                                                           gCCATGATGG‑CGCAGGTTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGtt   <  1:163875/149‑1 (MQ=255)                                                                                                                                            ccATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGATTAGTCGGCGGGCTGGCGGAGAAATGCTCTGCCAGCTCTTCGGTAGTGAGATAACCCTGCTGTTTCCCCAGTTCGATAATACCGTTGTGAGGtgg  >  2:717539/1‑147 (MQ=255)                                                                                                                                                       |                                                                                                                                           GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTG  >  NC_000913_3_bme_pgi/3560307‑3560593

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