Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	3,486,574	C→A	A152E (GCA→GAA)	crp →	cAMP‑activated global transcription factor, mediator of catabolite repression

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	3,486,574	0	C	A	100.0%	185.2 / NA	63	A152E (GCA→GAA)	crp	cAMP‑activated global transcription factor, mediator of catabolite repression
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (28/35);  total (28/35)

CGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACAC  >  NC_000913_3_bme_pgi/3486429‑3486719                                                                                                                                                  |                                                                                                                                                  cGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGaaaa                                                                                                                                                <  1:962583/149‑1 (MQ=255)   ccAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACTTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaac                                                                                                                                              <  2:1278696/149‑1 (MQ=255)    cAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCCCAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGGGGGCACCCTGGCGTTCCTCGACGTGCCGGGCCGCATTGCACAGACTCTGCTGACTCTGGAAAaacc                                                                                                                                             >  1:526265/1‑148 (MQ=255)     aaTTGATTCAGGTAACCCCGGACATTCTGATGCGTTTGTCTGCACAGAAGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaacaa                                                                                                                                            <  1:938070/149‑1 (MQ=255)     aaTTGATTCAGGTAAACCCGGACATTCTGATGCGTATGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaacaa                                                                                                                                            <  1:506768/149‑1 (MQ=255)     aaTTGATTCAGCTAAACCCGGACACGCTGAGGCGTTTGTCCGCACAGATGGCGCGTCGTCTGCAAGGCACTTCCGAGAAAGTGGTCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaacca                                                                                                                                            <  2:872439/149‑3 (MQ=255)     aaTTGATTAAGGTAAACCCGGACATTATGATTCGTTTGTCTGCACATATGGCGCTTAGTCTTCAAGTCACATCAGAGAAAGTGGGAAACATGCCGTTCCTCGACGTGACGGGCCGCATTGCAAAGAGTCTGGTGAATCTGTAAAaacaa                                                                                                                                            >  2:1317029/1‑149 (MQ=255)      aTTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaacaac                                                                                                                                           <  2:550114/149‑1 (MQ=255)      aTTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAaacaac                                                                                                                                           >  2:33333/1‑149 (MQ=255)        tGATTCAGGTAAACCCGGACCTTCTGATGCGTGTGGCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCa                                                                                                                                         <  1:425997/149‑1 (MQ=255)           ttCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGAc                                                                                                                                      <  1:856878/149‑1 (MQ=255)                gTAAACTCGGACATTCTGTTGCGTTTGTCTGCACAGATGGCGCGGCGGCTGCTAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGTCGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTAt                                                                                                                                 <  1:73227/149‑1 (MQ=255)                gTAAACCCGGACATTCTGATGCTTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAATTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTAt                                                                                                                                 <  1:1287654/149‑1 (MQ=255)                gTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTAt                                                                                                                                 <  1:579921/149‑1 (MQ=255)                                tGATGGGTTTGTCTGCACAGATGGCGCGTCTTCTGCAAGTGACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACTTGACGGGCAGCATTGAACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGt                                                                                                                 <  1:1358245/149‑1 (MQ=255)                                        ttGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAAt                                                                                                         <  1:1406861/149‑1 (MQ=255)                                        ttGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAAt                                                                                                         <  1:1357052/149‑1 (MQ=255)                                          gTCTGCACAGATGGCGAGTCGTCTGCAAGTCACTTCAGAGAAAGTGGTCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCAGGACGGTATGCAAATCa                                                                                                       >  1:825314/1‑149 (MQ=255)                                           tCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCACGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCaa                                                                                                      >  1:236870/1‑149 (MQ=255)                                           gcTGCACAGATGGCGCGGCGTCTGCAACTCACTTCAGAGAAAGTGGGCAACCGGGCGATCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCaa                                                                                                      <  2:265793/148‑1 (MQ=255)                                            cTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCaaa                                                                                                     <  2:912371/149‑1 (MQ=255)                                               cacaGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAATCATCCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAAtt                                                                                                  >  2:1298773/1‑149 (MQ=255)                                               cacaGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAAtt                                                                                                  >  1:893716/1‑149 (MQ=255)                                               cacaGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTAAGACTCACCCGGCCGGTATGCAAATCAAAAtt                                                                                                  >  2:657587/1‑149 (MQ=255)                                                acaGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGAAACCAGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTa                                                                                                 >  1:1085869/1‑149 (MQ=255)                                                  aGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTAcc                                                                                               <  2:193746/149‑1 (MQ=255)                                                   gATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTAccc                                                                                              <  1:481845/149‑1 (MQ=255)                                                       gcgcGGCGCCTGCAAGTCACTTCAGAGAAAGTGGGCTCCCTGGCGTTCCTCGACGGGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCa                                                                                          <  2:477577/149‑1 (MQ=255)                                                         gcgTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGTCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGATATTACTAACACGGACGGTATGCAAATCAAAATTAAACGTAAgg                                                                                        >  1:285336/1‑149 (MQ=255)                                                                 tGCAAGTCACTTCAGAGAAAGTGGGCAAGCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGt                                                                                <  1:902199/149‑1 (MQ=255)                                                                  gCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTc                                                                               <  1:1159202/149‑1 (MQ=255)                                                                  gCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCAGTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACACGTCAGGAAATTGGTc                                                                               >  1:809623/1‑149 (MQ=255)                                                                    aaGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATGGCACAGACTCTGCTGAATCTGGAAAAAAAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACAAGTCAGGAAATTGGTCAg                                                                             >  2:1150144/1‑149 (MQ=255)                                                                     aGTTACTTCAGAGAATGTGGGCAACCTTGCGTTCCTAGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAATACAACCAGACGCTATGACTCTCCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGa                                                                            <  1:203046/149‑1 (MQ=255)                                                                           ttCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGCCGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGAGTGTCg                                                                      >  1:1058049/1‑149 (MQ=255)                                                                           ttCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCg                                                                      >  2:918565/1‑149 (MQ=255)                                                                           ttCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCg                                                                      >  1:1134917/1‑149 (MQ=255)                                                                           tgcAGAGCAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCAGTGCACAGACTCTGCTGAATCTGGAAAAACAACCCGCCGCTATGCCTCCCCCGGACGGTATGCACATCAAAATTACCCGTCTGGAAATTGGTCAGATTGTCg                                                                      >  2:1113286/3‑149 (MQ=255)                                                                              agagAAAGTGGGCAACCGGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCt                                                                   <  2:1222045/149‑1 (MQ=255)                                                                               gagaAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTg                                                                  >  1:1470801/1‑149 (MQ=255)                                                                                agaAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCCCAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTAAGACTCTCCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGt                                                                 >  1:481656/1‑149 (MQ=255)                                                                                agaAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGt                                                                                            >  2:1266440/1‑122 (MQ=255)                                                                                agaAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGt                                                                                            <  1:1266440/122‑1 (MQ=255)                                                                                          cAACCTGGCGATCCTCGACGACACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAAACAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCTTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGaaa                                                       <  2:865546/149‑1 (MQ=255)                                                                                              cTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGt                                                   <  1:33333/149‑1 (MQ=255)                                                                                              cTGGCGGGCCTCGACGTGACGGGCCGCATTGAAGAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTGACCGGGACGGTATGCATATCAAAATTCCCCGTCAGGAAATTGGTCAGCTTGTCGGCTGTTCTCGTGAAACCGt                                                   <  2:526265/149‑1 (MQ=255)                                                                                                       cTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCAt                                          >  1:541552/1‑149 (MQ=255)                                                                                                         cGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTc                                        <  2:1085869/149‑1 (MQ=255)                                                                                                               gACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCCGCCGGTATGACTCACCCGGCCGGTTTCCAATGCAAATTTACCCGGATGGAATTTGTGCAGATTGTCGGCGGTTCACGTGAACCCGTGGGCCGCATTCTGAAgt                                  >  2:720712/1‑148 (MQ=255)                                                                                                               gACGGGCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCATACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGACAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGa                                  >  2:497724/1‑149 (MQ=255)                                                                                                                    gCCGCATTGCACAGACTCTGCTGAATCTGGAAAAACAACCAGCCGCTATGACTCACCCGGCCGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTg                             >  1:1221747/1‑149 (MQ=255)                                                                                                                     gcGCATGGCCCAGCCTCGGCTGAATCTGGAAAAACACCCAGCCGCTATGACTCACCCTGACGGTAGGCAACTCAAAATTACCCGTCAGGAAATTGGTCAGCTTGTCGGCTGTTATCGTGAAACCGTGGGACGCATTCTGAAGATGCTgg                            <  2:877440/148‑1 (MQ=255)                                                                                                                      ctcATTGCACATATTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCATATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGa                           <  1:1150144/147‑1 (MQ=255)                                                                                                                          ttGCACAGACTCTGCTGAATGTTGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGAt                       <  1:561065/149‑1 (MQ=255)                                                                                                                           tccacaGACTCTGCTGAATCTGGAAAAACAACCAGACCCTAAGACTCACCCGGACGGTAAGCAACTCAAAATTACCCGCCAGGAAATTGCTCAGATAGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATc                      <  2:220166/147‑1 (MQ=255)                                                                                                                             cacaGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCAAGTGAAACCGTGGTACGCATTCTGAAGATGCTGGCAGATCAg                    >  1:978341/1‑149 (MQ=255)                                                                                                                              acaGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGACGATCAGa                   >  2:1218705/1‑149 (MQ=255)                                                                                                                              acaGACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCAGTCTGAAGATGCTGGAAGATCAGa                   >  1:1281960/1‑149 (MQ=255)                                                                                                                               caGCCTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGGACATCAAAATTACCCGTCAGGAAATTGGGCAGATTGTCGGCGGTTCTCGTGAAACCGTGGTACGCTTTCTGAAGATGATGGCAGagccgac                  >  2:90564/1‑143 (MQ=255)                                                                                                                                 gACTCTGCTGAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAAcc                <  2:1134917/149‑1 (MQ=255)                                                                                                                                     atgctgATTGTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGCCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGAt            <  2:596347/148‑1 (MQ=255)                                                                                                                                      tgctgAATCTGGAAAAACAACCAGACGCTATGACTCACCCGGACGGTATCCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATGCTGACGATTCTGGAAGATCAGAAACTGAtc           <  2:280686/149‑1 (MQ=255)                                                                                                                                            aTCAGGAAAAACAACCAGACGATATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTGCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGca     <  1:520384/149‑1 (MQ=255)                                                                                                                                               tGGAAAAACAACCAGACGCTATGACTCCCCCGGACGGTATGCAAATCAAAATTACCGGTCAGGTAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGCACCTGATCTCCGcactc  >  1:831545/1‑147 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACAC  >  NC_000913_3_bme_pgi/3486429‑3486719

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