Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,222,930	G→A	D100N (GAT→AAT)	ebgA →	evolved beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,222,930	0	G	A	100.0%	278.5 / NA	80	D100N (GAT→AAT)	ebgA	evolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (43/37);  total (43/37)

ATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCGATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACGCTGATTAAATTTGACGGCGTCGAAACCTATTTTGAAGTCTATGTTAAC  >  NC_000913/3222796‑3223064                                                                                                                                       |                                                                                                                                       aTCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGACAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAAtgt                                                                                                                                    >  1:1282371/1‑139 (MQ=255)   ccGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATATGc                                                                                                                                  <  1:854548/139‑1 (MQ=255)   ccGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGATGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGCCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGc                                                                                                                                  >  1:412464/1‑139 (MQ=255)              aCCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTccc                                                                                                                        <  2:907717/138‑1 (MQ=255)              aCCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTccc                                                                                                                        >  1:907717/1‑138 (MQ=255)                       cTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAAc                                                                                                              >  1:1293885/1‑139 (MQ=255)                        ttCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTcc                                                                                                                         <  2:1725632/127‑1 (MQ=255)                        ttCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTcc                                                                                                                         >  1:1725632/1‑127 (MQ=255)                        ttCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAAcc                                                                                                             >  2:1371271/1‑139 (MQ=255)                            ccTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAAc                                                                                                         <  2:1239823/139‑1 (MQ=255)                            ccTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAAc                                                                                                         <  2:1036751/139‑1 (MQ=255)                             ctctGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAAcc                                                                                                        >  1:1036821/1‑139 (MQ=255)                                tGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAAc                                                                                                         <  1:1360474/135‑1 (MQ=255)                                tGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAAc                                                                                                         >  2:1360474/1‑135 (MQ=255)                                tGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGt                                                                                                     <  2:1282371/139‑1 (MQ=255)                                tGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGt                                                                                                     >  2:1442210/1‑139 (MQ=255)                                         ggCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCaa                                                                                            <  2:122189/139‑1 (MQ=255)                                          gCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAAc                                                                                           >  2:1680053/1‑139 (MQ=255)                                           ccgACCGGGGGCATAATTCCGGCCCCGGCCAGGGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACg                                                                                          <  2:1168440/137‑1 (MQ=255)                                           cTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACg                                                                                          <  1:1427277/139‑1 (MQ=255)                                                        tatTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGt                                                                                                                               <  2:355384/89‑1 (MQ=255)                                                        tatTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGt                                                                                                                               >  1:355384/1‑89 (MQ=255)                                                        tatTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTc                                                                             >  2:25513/1‑139 (MQ=255)                                                                tCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGc                                                                                                                                  >  2:266519/1‑78 (MQ=255)                                                                tCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGc                                                                                                                                  <  1:266519/78‑1 (MQ=255)                                                                 ccccGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGc                                                                                                   >  2:1269963/1‑108 (MQ=255)                                                                 ccccGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGc                                                                                                   <  1:1269963/108‑1 (MQ=255)                                                                 ccccGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGc                                                                    >  1:1504516/1‑139 (MQ=255)                                                                    cGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTgg                                                                 >  2:1538155/1‑139 (MQ=255)                                                                     gCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGc                                                                >  2:113854/1‑139 (MQ=255)                                                                      ccATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCa                                                               <  1:1371271/139‑1 (MQ=255)                                                                       gccatgGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCt                                                                              >  1:25229/5‑123 (MQ=255)                                                                       gccatgGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCt                                                                              <  2:25229/119‑1 (MQ=255)                                                                       cATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACg                                                                                          <  2:994322/111‑1 (MQ=255)                                                                       cATGTGGCAAATGGAAGGACACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACg                                                                                          >  1:994322/1‑111 (MQ=255)                                                                         tgtgGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAggg                                                            <  2:1036821/139‑1 (MQ=255)                                                                                   ggAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCa                                                               >  2:709727/1‑126 (MQ=255)                                                                                   ggAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCa                                                               <  1:709727/126‑1 (MQ=255)                                                                                   ggAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACg                                                  <  1:1442210/139‑1 (MQ=255)                                                                                       ggTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTg                                                                  <  1:1223824/119‑1 (MQ=255)                                                                                       ggTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTg                                                                  >  2:1223824/1‑119 (MQ=255)                                                                                       ggTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACGCTGa                                              <  1:113854/139‑1 (MQ=255)                                                                                            cGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACg                                                  <  1:768399/130‑1 (MQ=255)                                                                                            cGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACg                                                  >  2:768399/1‑130 (MQ=255)                                                                                            cGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACGCTGATTaaa                                         >  1:671347/1‑139 (MQ=255)                                                                                             ggCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGAt                                                                                                                 >  2:1534904/1‑66 (MQ=255)                                                                                             ggCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGAt                                                                                                                 <  1:1534904/66‑1 (MQ=255)                                                                                             ggCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGAc                                                                       >  2:290136/1‑108 (MQ=255)                                                                                             ggCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGAc                                                                       <  1:290136/108‑1 (MQ=255)                                                                                              gCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACGCTGATTAAAtt                                       >  1:1045248/1‑139 (MQ=255)                                                                                                              cgacgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACg                                                                                          <  1:889531/72‑1 (MQ=255)                                                                                                              cgacgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACg                                                                                          >  2:889531/1‑72 (MQ=255)                                                                                                              cgacgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACGCTGATTAAATTTGACGGCGTCGAAAcc                       >  2:1188256/1‑139 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAAc                                                                                                         >  1:612534/1‑55 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAAc                                                                                                         <  2:612534/55‑1 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTg                                                                                                    <  1:200734/60‑1 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTg                                                                                                    >  2:200734/1‑60 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAAc                                                                                           >  1:1692000/1‑69 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAAc                                                                                           <  2:1692000/69‑1 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTAtttt                                                                                    >  1:202159/1‑76 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTAtttt                                                                                    <  2:202159/76‑1 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCAcc                                                                                <  2:1298523/80‑1 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCAcc                                                                                >  1:1298523/1‑80 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAgg                                                             >  1:418080/1‑99 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAgg                                                             <  2:418080/99‑1 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACGCTGATTAAATTTGACGGCGTCGAAACCTa                     >  1:1276824/1‑139 (MQ=255)                                                                                                                acgaAGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACGCTGATTAAATTTGACGGCGTCGAAACCTa                     >  2:65306/1‑139 (MQ=255)                                                                                                                   aaGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCa                                                                                  <  1:1142047/75‑1 (MQ=255)                                                                                                                   aaGGTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCa                                                                                  >  2:1142047/1‑75 (MQ=255)                                                                                                                     ggTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACGCTGATTAAATTTGACGGCGTCGAAACCTATTTTg                <  1:25513/139‑1 (MQ=255)                                                                                                                     ggTTTTCCGTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACGCTGATTAAATTTGACGGCGTCGAAACCTATTTTg                >  2:1123568/1‑139 (MQ=255)                                                                                                                             gTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACg                                                  <  1:1020176/97‑1 (MQ=255)                                                                                                                             gTTCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACg                                                  >  2:1020176/1‑97 (MQ=255)                                                                                                                               tCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATc                                                                                              >  2:201347/1‑51 (MQ=255)                                                                                                                               tCCCCATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATc                                                                                              <  1:201347/51‑1 (MQ=255)                                                                                                                                ccccATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTgg                                                                 <  1:9897/79‑1 (MQ=255)                                                                                                                                ccccATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTgg                                                                 >  2:9897/1‑79 (MQ=255)                                                                                                                                   cATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACGCTGATTAAATTTGACGGCGTCGAAACCTAtttt                 <  1:1349725/124‑1 (MQ=255)                                                                                                                                   cATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACGCTGATTAAATTTGACGGCGTCGAAACCTAtttt                 >  2:1349725/1‑124 (MQ=255)                                                                                                                                   cATCAATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACGCTGATTAAATTTGACGGCGTCGAAACCTATTTTGAAGTCTATGTTAAc  >  2:772016/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                       ATCCGCTGCAAGTACCAGAAGCCTTCACCTCTGAGTTAATGGCTGACTGGGGGCATATTACCGTCCCCGCCATGTGGCAAATGGAAGGTCACGGCAAACTGCAATATACCGACGAAGGTTTTCCGTTCCCCATCGATGTGCCGTTTGTCCCCAGCGATAACCCAACCGGTGCCTATCAACGTATTTTCACCCTCAGCGACGGCTGGCAGGGTAAACAGACGCTGATTAAATTTGACGGCGTCGAAACCTATTTTGAAGTCTATGTTAAC  >  NC_000913/3222796‑3223064

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