Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,581,809	G→A	P414S (CCG→TCG)	hsdM ←	type I restriction enzyme EcoKI methyltransferase component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,581,809	0	G	A	100.0%	339.2 / NA	95	P414S (CCG→TCG)	hsdM	type I restriction enzyme EcoKI methyltransferase component
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (42/53);  total (42/53)

CGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTC  >  NC_000913/4581677‑4581942                                                                                                                                     |                                                                                                                                      cggtgcggATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTaaa                                                                                                                                <  1:45142/140‑1 (MQ=255) cggtgcggATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTaaa                                                                                                                                <  1:1562716/140‑1 (MQ=255)     gcggATGCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGt                                                                                                                            >  1:911643/1‑140 (MQ=255)     gcggATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGt                                                                                                                            >  2:243979/1‑140 (MQ=255)         aTCCACTCACGGCTGAACTTGCGCCAGCGGCTAGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCgg                                                                                                                        >  1:759150/1‑140 (MQ=255)           ccACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTTCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGt                                                                                                                      <  2:1416570/140‑1 (MQ=255)            cACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTc                                                                                                                     <  1:1037141/140‑1 (MQ=255)              cTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGt                                                                                                                      <  2:150098/137‑1 (MQ=255)              cTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGt                                                                                                                      >  1:150098/1‑137 (MQ=255)              cTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCtt                                                                                                                   <  2:1466819/140‑1 (MQ=255)              cTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCtt                                                                                                                   <  1:243979/140‑1 (MQ=255)              cTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCtt                                                                                                                   >  1:889665/1‑140 (MQ=255)              cTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCtt                                                                                                                   >  1:794454/1‑140 (MQ=255)              cTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCtt                                                                                                                   <  2:794454/140‑1 (MQ=255)              cTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCtt                                                                                                                   >  2:954732/1‑140 (MQ=255)                  cGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGcc                                                                                                               >  2:825196/1‑140 (MQ=255)                   ggCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCa                                                                                                              >  1:861782/1‑140 (MQ=255)                   ggCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCa                                                                                                              >  1:1027933/1‑140 (MQ=255)                     cTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATa                                                                                                            >  2:292392/1‑140 (MQ=255)                     cTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATa                                                                                                            >  2:1184071/1‑140 (MQ=255)                       gAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATaca                                                                                                          <  1:706482/140‑1 (MQ=255)                       gAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATaca                                                                                                          <  2:889665/140‑1 (MQ=255)                       gAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATaca                                                                                                          <  2:1536176/140‑1 (MQ=255)                             gcgcCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCt                                                                                                    >  2:356955/1‑140 (MQ=255)                             gcgcCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCt                                                                                                    >  1:521362/1‑140 (MQ=255)                             gcgcCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCt                                                                                                    >  1:101758/1‑140 (MQ=255)                               gcCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCa                                                                                                  <  1:1308065/140‑1 (MQ=255)                                    cGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACgg                                                                                             <  1:1184071/140‑1 (MQ=255)                                    cGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACgg                                                                                             >  2:285556/1‑140 (MQ=255)                                     ggCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGc                                                                                            <  2:1027933/140‑1 (MQ=255)                                               aaGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCt                                                                                  <  2:779493/140‑1 (MQ=255)                                               aaGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCt                                                                                  <  1:431681/140‑1 (MQ=255)                                                 gATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCtcg                                                                                <  2:1392987/140‑1 (MQ=255)                                                 gATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCACTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCtcg                                                                                <  2:438959/140‑1 (MQ=255)                                                    gCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCtcgtcg                                                                             >  2:133037/1‑140 (MQ=255)                                                     cTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCgg                                                                            >  1:984373/1‑140 (MQ=255)                                                        gTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTa                                                                          >  1:1097749/1‑139 (MQ=255)                                                        gTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTa                                                                          <  2:1097749/139‑1 (MQ=255)                                                         tCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTaaa                                                                        <  2:703264/140‑1 (MQ=255)                                                         tCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTaaa                                                                        <  1:1276333/140‑1 (MQ=255)                                                           ggTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACg                                                                      >  1:1401794/1‑140 (MQ=255)                                                            gtgtTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACgg                                                                     >  1:1431516/1‑140 (MQ=255)                                                            gtgtTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACgg                                                                     >  1:1582676/1‑140 (MQ=255)                                                             tgtTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTTGGTAAACGgt                                                                    >  1:713234/1‑140 (MQ=255)                                                             tgtTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGgt                                                                    >  2:1176599/1‑140 (MQ=255)                                                                  ttGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTgcg                                                               <  2:861782/140‑1 (MQ=255)                                                                       ctctTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGc                                                          >  2:442389/1‑140 (MQ=255)                                                                       ctctTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGc                                                          >  2:1237193/1‑140 (MQ=255)                                                                         ctTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCg                                                        >  1:1087471/1‑140 (MQ=255)                                                                          ttCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAACCGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGa                                                       <  2:1591671/140‑1 (MQ=255)                                                                          ttCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGa                                                       >  1:1068285/1‑140 (MQ=255)                                                                          ttCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGa                                                       >  2:1092111/1‑140 (MQ=255)                                                                            cGCTGTCGGTAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGaaa                                                     >  2:257058/1‑140 (MQ=255)                                                                               tGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTc                                                  >  1:990922/1‑140 (MQ=255)                                                                                gTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCg                                                 <  1:954732/140‑1 (MQ=255)                                                                                 tCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCgg                                                >  1:411797/1‑140 (MQ=255)                                                                                   ggCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATCCACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCa                                              <  2:1431516/140‑1 (MQ=255)                                                                                   ggCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCa                                              <  1:285556/140‑1 (MQ=255)                                                                                   ggCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCa                                              <  2:679291/140‑1 (MQ=255)                                                                                   ggCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCa                                              <  2:18893/140‑1 (MQ=255)                                                                                   ggCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCa                                              <  1:356955/140‑1 (MQ=255)                                                                                      aaCTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCatat                                           <  2:521362/140‑1 (MQ=255)                                                                                      aaCTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCatat                                           <  2:1380875/140‑1 (MQ=255)                                                                                        cTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTg                                         <  1:8209/140‑1 (MQ=255)                                                                                           ccGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTa                                      >  1:1198103/1‑140 (MQ=255)                                                                                            cGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTTGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTAc                                     <  2:713234/140‑1 (MQ=255)                                                                                             gTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACg                                    <  1:292392/140‑1 (MQ=255)                                                                                             gTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACg                                    <  1:1176599/140‑1 (MQ=255)                                                                                              tctcTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGc                                   <  1:257058/140‑1 (MQ=255)                                                                                               ctctTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCa                                  <  1:133037/140‑1 (MQ=255)                                                                                                  ttCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGt                               <  2:645015/140‑1 (MQ=255)                                                                                                   tCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTAc                                     <  1:1416378/133‑1 (MQ=255)                                                                                                   tCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTc                              >  1:115649/1‑140 (MQ=255)                                                                                                   tCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGATTGCCGAAACTCGGCATATTGGTAc                                     >  2:1416378/1‑133 (MQ=255)                                                                                                      gCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATa                           <  1:825196/140‑1 (MQ=255)                                                                                                      gCGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATa                           <  1:1219290/140‑1 (MQ=255)                                                                                                       cGTTAAAACTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATac                          <  2:734376/140‑1 (MQ=255)                                                                                                            aaaCTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCca                     <  2:1401794/140‑1 (MQ=255)                                                                                                               cTCCATTCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCcacac                  <  2:798613/140‑1 (MQ=255)                                                                                                                    ttCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACAcatcat             <  1:1237193/140‑1 (MQ=255)                                                                                                                    ttCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACAcatcat             >  1:1505314/1‑140 (MQ=255)                                                                                                                    ttCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACAcatcat             >  2:548991/1‑140 (MQ=255)                                                                                                                    ttCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACAcatcat             >  1:373189/1‑140 (MQ=255)                                                                                                                    ttCACCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACAcatcat             <  1:1092111/140‑1 (MQ=255)                                                                                                                       aCCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCgg          <  2:411797/140‑1 (MQ=255)                                                                                                                       aCCTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCgg          <  2:1505314/140‑1 (MQ=255)                                                                                                                         cTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTa        <  1:546002/140‑1 (MQ=255)                                                                                                                         cTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTa        >  2:157286/1‑140 (MQ=255)                                                                                                                         cTTCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTa        <  1:442389/140‑1 (MQ=255)                                                                                                                          ttCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTTAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTAc       <  2:990922/140‑1 (MQ=255)                                                                                                                          ttCAGTGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTAc       <  2:1395368/140‑1 (MQ=255)                                                                                                                               tGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTc  >  2:386212/1‑140 (MQ=255)                                                                                                                               tGCGCGAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTc  >  1:335725/1‑140 (MQ=255)                                                                                                                                 cgcgAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTAc       >  2:692539/1‑133 (MQ=255)                                                                                                                                 cgcgAGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTAc       <  1:692539/133‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      CGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTC  >  NC_000913/4581677‑4581942

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