Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,604,202	G→A	*242* (TAG→TAA)  E15K (GAA→AAA)	yjjQ → bglJ →	putative transcriptional regulator bgl operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,202	0	G	A	100.0%	256.6 / NA	78	*242* (TAG→TAA)  E15K (GAA→AAA)	yjjQ bglJ	putative transcriptional regulator bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/39);  total (39/39)

TCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA  >  NC_000913/4604070‑4604337                                                                                                                                     |                                                                                                                                        tCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGTGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                                                   <  1:1167419/139‑1 (MQ=255) tCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                                                   <  2:308647/139‑1 (MQ=255)      gCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                              <  2:1661255/139‑1 (MQ=255)      gCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                              <  2:1038231/139‑1 (MQ=255)          cGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                                                   <  1:391854/130‑1 (MQ=255)          cGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                                                   >  2:391854/1‑130 (MQ=255)          cGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGa                                                                                                                          >  1:1469890/1‑139 (MQ=255)           gCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                         <  1:1590413/139‑1 (MQ=255)           gCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                         <  1:1524279/139‑1 (MQ=255)            cTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGagt                                                                                                                        <  1:948316/139‑1 (MQ=255)                aaTGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                                    <  1:1692942/139‑1 (MQ=255)                         aTGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                                           <  2:1469890/139‑1 (MQ=255)                              tAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCAt                                                                                                                            <  2:1721032/117‑1 (MQ=255)                              tAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCAt                                                                                                                            >  1:1721032/1‑117 (MQ=255)                                     cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTAt                                                                                                                   <  2:799383/119‑1 (MQ=255)                                     cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTAt                                                                                                                   >  1:799383/1‑119 (MQ=255)                                     cgcgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCa                                                                                               >  1:1578551/1‑139 (MQ=255)                                       cgcAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGGTTATTCAGa                                                                                             >  1:462907/1‑139 (MQ=255)                                                        cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                                                   <  1:1069087/84‑1 (MQ=255)                                                        cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAAt                                                                                                                                   >  2:1069087/1‑84 (MQ=255)                                                        cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGa                                                                                                                          >  2:252669/1‑93 (MQ=255)                                                        cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGa                                                                                                                          <  1:252669/93‑1 (MQ=255)                                                        cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTAtt                                                                                                            >  2:227100/1‑107 (MQ=255)                                                        cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTAtt                                                                                                            <  1:227100/107‑1 (MQ=255)                                                        cttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAAccc                                                                            >  1:854099/1‑139 (MQ=255)                                                                gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCAt                                                                                                                            <  2:1490743/83‑1 (MQ=255)                                                                gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCAt                                                                                                                            >  1:1490743/1‑83 (MQ=255)                                                                gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                         >  1:614788/1‑86 (MQ=255)                                                                gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                         <  2:614788/86‑1 (MQ=255)                                                                 ccTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                                           >  1:249981/1‑99 (MQ=255)                                                                 ccTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTg                                                                                                           <  2:249981/99‑1 (MQ=255)                                                                 ccTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCt                                                                   >  2:1642364/1‑139 (MQ=255)                                                                   ttCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                                 <  2:854099/139‑1 (MQ=255)                                                                   ttCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                                 <  2:1706668/139‑1 (MQ=255)                                                                   ttCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                                 <  1:1642364/139‑1 (MQ=255)                                                                        tttATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                              >  1:1137201/1‑73 (MQ=255)                                                                        tttATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                              <  2:1137201/73‑1 (MQ=255)                                                                          tATGGAGATGCCGCAGAATGGAACACAGCCGAATTAATAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                              >  1:1203853/1‑71 (MQ=255)                                                                          tATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTc                                                                                                                              <  2:1203853/71‑1 (MQ=255)                                                                                gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAg                                                    >  2:607261/1‑139 (MQ=255)                                                                                gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAg                                                    >  1:1531600/1‑139 (MQ=255)                                                                                gaTGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAg                                                    >  2:1043553/1‑139 (MQ=255)                                                                                  tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTtata                                                                       <  2:200293/118‑1 (MQ=255)                                                                                  tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTtata                                                                       >  1:200293/1‑118 (MQ=255)                                                                                  tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                                  >  1:1601269/1‑123 (MQ=255)                                                                                  tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTc                                                                  <  2:1601269/123‑1 (MQ=255)                                                                                  tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTc                                                  >  2:895282/1‑139 (MQ=255)                                                                                  tGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTc                                                  >  1:802042/1‑139 (MQ=255)                                                                                   gCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCa                                                 <  1:607261/139‑1 (MQ=255)                                                                                    ccGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                         <  2:1348520/66‑1 (MQ=255)                                                                                    ccGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                         >  1:1348520/1‑66 (MQ=255)                                                                                      gCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGa                                              <  1:895282/139‑1 (MQ=255)                                                                                            tGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGc                                                                                                                                 >  2:186414/1‑50 (MQ=255)                                                                                            tGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGc                                                                                                                                 <  1:186414/50‑1 (MQ=255)                                                                                               aaCACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCa                                     >  2:17313/1‑139 (MQ=255)                                                                                                  acaGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGa                                  <  1:17313/139‑1 (MQ=255)                                                                                                    aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                         <  2:779592/50‑1 (MQ=255)                                                                                                    aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGag                                                                                                                         >  1:779592/1‑50 (MQ=255)                                                                                                    aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                                    >  2:196700/1‑55 (MQ=255)                                                                                                    aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                                    <  1:196700/55‑1 (MQ=255)                                                                                                    aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATg                                <  2:1531600/139‑1 (MQ=255)                                                                                                    aGCCGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATg                                >  1:733219/1‑139 (MQ=255)                                                                                                      ccGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGcc                              <  2:1578551/139‑1 (MQ=255)                                                                                                      ccGAATTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGcc                              >  1:1092608/1‑139 (MQ=255)                                                                                                         aaTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATg                           >  1:1550549/1‑139 (MQ=255)                                                                                                          aTTAAGAAGAGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGt                          <  2:733219/139‑1 (MQ=255)                                                                                                              agaagaGAAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTcgcg                      <  2:1550549/139‑1 (MQ=255)                                                                                                                   agaAATGTCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGAtct                 >  1:72886/1‑139 (MQ=255)                                                                                                                        tGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                                    <  1:800466/35‑1 (MQ=255)                                                                                                                        tGTCGCACTCATAAAAAAATGCGTCATGAgtagta                                                                                                                    >  2:800466/1‑35 (MQ=255)                                                                                                                          tCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                                 >  2:1080149/1‑84 (MQ=255)                                                                                                                          tCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTcc                                                                 <  1:1080149/84‑1 (MQ=255)                                                                                                                          tCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCAttt                                       >  2:952184/1‑110 (MQ=255)                                                                                                                          tCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCAttt                                       <  1:952184/110‑1 (MQ=255)                                                                                                                          tCGCACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTgc          >  1:84362/1‑139 (MQ=255)                                                                                                                             cACTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTgcggc       >  2:1748499/1‑139 (MQ=255)                                                                                                                               cTCATAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGcggcgg     >  2:438240/1‑139 (MQ=255)                                                                                                                                  aTAAAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCa  >  2:355509/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                        TCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCGCAGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCA  >  NC_000913/4604070‑4604337

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