Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,584,837	A→G	V642A (GTT→GCT)	hsdR ←	type I restriction enzyme EcoKI endonuclease component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,584,837	0	A	G	100.0%	258.1 / NA	70	V642A (GTT→GCT)	hsdR	type I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/34);  total (36/34)

AGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATC  >  NC_000913/4584702‑4584970                                                                                                                                        |                                                                                                                                      aGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCcgc                                                                                                                                   >  1:963072/1‑140 (MQ=255) aGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCcgc                                                                                                                                   <  1:815024/140‑1 (MQ=255) aGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCcgc                                                                                                                                   <  1:312453/140‑1 (MQ=255)  gggTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCcgcg                                                                                                                                  <  1:603009/140‑1 (MQ=255)     tgtCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGcggta                                                                                                                               <  2:171351/140‑1 (MQ=255)     tgtCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGcggta                                                                                                                               >  1:171351/1‑140 (MQ=255)         aTTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGcggtacggt                                                                                                                           <  1:373391/140‑1 (MQ=255)         aTTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGcggtacggt                                                                                                                           <  2:618634/140‑1 (MQ=255)                aCTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGtata                                                                                                                    <  1:827587/140‑1 (MQ=255)                aCTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGtata                                                                                                                    >  1:458198/1‑140 (MQ=255)                 cTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATaa                                                                                                                   >  2:379410/1‑140 (MQ=255)                 cTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATaa                                                                                                                   >  1:525016/1‑140 (MQ=255)                 cTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATaa                                                                                                                   >  1:511894/1‑140 (MQ=255)                         tGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAAc                                                                                                           >  2:313701/1‑140 (MQ=255)                            gggCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCgg                                                                                                        >  1:358729/1‑140 (MQ=255)                            gggCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCgg                                                                                                        <  1:150093/140‑1 (MQ=255)                               cTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTc                                                                                                     >  1:949897/1‑140 (MQ=255)                                   tGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCg                                                                                                 <  1:365103/140‑1 (MQ=255)                                   tGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCg                                                                                                 <  1:317474/140‑1 (MQ=255)                                        tctACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAAt                                                                                            >  1:461451/1‑140 (MQ=255)                                          tACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTTCGGTAGGTATAACGGTAAACCGGCTCGCCGAAAAACt                                                                                          >  2:690622/1‑140 (MQ=255)                                             cTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGca                                                                                       <  1:785745/140‑1 (MQ=255)                                             cTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGca                                                                                       <  2:458198/140‑1 (MQ=255)                                             cTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGca                                                                                       <  2:525016/140‑1 (MQ=255)                                               gCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGcaca                                                                                     >  1:142639/1‑140 (MQ=255)                                                cTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAg                                                                                    <  1:313701/140‑1 (MQ=255)                                                    ccTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCt                                                                                          <  1:215319/130‑1 (MQ=255)                                                    ccTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCt                                                                                          >  2:215319/1‑130 (MQ=255)                                                    ccTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGca                                                                                       >  1:161389/1‑133 (MQ=255)                                                    ccTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGca                                                                                       <  2:161389/133‑1 (MQ=255)                                                      tttGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTa                                                                              >  1:866317/1‑140 (MQ=255)                                                       ttGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAg                                                                             <  1:944429/140‑1 (MQ=255)                                                              gaTAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCggg                                                                      <  1:721614/140‑1 (MQ=255)                                                                tAAGCCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGGCGAAAATCTGCACAGTATGTAGCGCCGGGGt                                                                    <  2:253079/140‑1 (MQ=255)                                                                      cccTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGt                                                              <  1:379410/140‑1 (MQ=255)                                                                              cgcgTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGa                                                      >  1:429738/1‑140 (MQ=255)                                                                                          gatgatCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGTCGGGGTGGCGGTGAGAGCGATTTTTACCGCAt                                          <  2:606134/140‑1 (MQ=255)                                                                                          gatgatCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCAt                                          <  1:773693/140‑1 (MQ=255)                                                                                                 tGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTa                                                >  2:273651/1‑127 (MQ=255)                                                                                                 tGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTa                                                <  1:273651/127‑1 (MQ=255)                                                                                                 tGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTg                                   >  1:40367/1‑140 (MQ=255)                                                                                                   aaTAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTa                                                <  1:496702/125‑1 (MQ=255)                                                                                                   aaTAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTa                                                >  2:496702/1‑125 (MQ=255)                                                                                                   aaTAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGaa                                      <  2:816584/135‑1 (MQ=255)                                                                                                   aaTAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGaa                                      >  1:816584/1‑135 (MQ=255)                                                                                                        gCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCgaga                            >  1:298710/1‑140 (MQ=255)                                                                                                        gCGGATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCgaga                            >  2:481524/1‑140 (MQ=255)                                                                                                           gATCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATg                         <  2:358729/140‑1 (MQ=255)                                                                                                             tCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCg                       >  2:17657/1‑140 (MQ=255)                                                                                                             tCCTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCg                       >  1:432961/1‑140 (MQ=255)                                                                                                              ccTGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGa                      <  2:40367/140‑1 (MQ=255)                                                                                                                tGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACg                    <  2:429738/140‑1 (MQ=255)                                                                                                                tGGTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACg                    <  1:481524/140‑1 (MQ=255)                                                                                                                 ggTCGATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACgg                   >  2:904202/1‑140 (MQ=255)                                                                                                                     gATCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAgg               <  2:298710/140‑1 (MQ=255)                                                                                                                      aTCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGc              >  1:1028579/1‑140 (MQ=255)                                                                                                                      aTCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGc              >  2:243290/1‑140 (MQ=255)                                                                                                                      aTCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGc              >  2:165095/1‑140 (MQ=255)                                                                                                                      aTCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGc              >  1:363982/1‑140 (MQ=255)                                                                                                                      aTCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGACGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATg                         <  1:895710/129‑1 (MQ=255)                                                                                                                      aTCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGACGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATg                         >  2:895710/1‑129 (MQ=255)                                                                                                                      aTCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCACCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGc              >  1:140083/1‑140 (MQ=255)                                                                                                                       tCAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCa             <  2:511894/140‑1 (MQ=255)                                                                                                                        cAGAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCag            <  2:949897/140‑1 (MQ=255)                                                                                                                          gAAAACCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCagag          <  2:1028579/140‑1 (MQ=255)                                                                                                                           aaaaCCGTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAgaga         >  1:657109/1‑140 (MQ=255)                                                                                                                                 gTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAAt   >  1:256364/1‑140 (MQ=255)                                                                                                                                 gTCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAAt   >  1:1019362/1‑140 (MQ=255)                                                                                                                                  tCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACg                    <  1:753644/122‑1 (MQ=255)                                                                                                                                  tCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACg                    >  2:753644/1‑122 (MQ=255)                                                                                                                                  tCGATAGCCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATc  >  1:773434/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                      AGGGTGTCATTGATCACTTCTCCCTGCGGGCTGATGCGCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGATCCTGGTCGATCAGAAAACCGTCGATAACCGCGGTACGGTAGGTATAACGGTAAACCGGCTCGCCGAAAATCTGCACAGTATGTAGCGCCGGGGTGGCGGTGAGAGCGATTTTTACCGCATCGAAGTGATCGAGAATGCGACGGTAGGCAGAGACGTAATC  >  NC_000913/4584702‑4584970

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