Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,577,854	A→G	G35G (GGT→GGC)	mcrC ←	5‑methylcytosine‑specific restriction enzyme subunit McrC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,577,854	0	A	G	100.0%	276.0 / NA	75	G35G (GGT→GGC)	mcrC	5‑methylcytosine‑specific restriction enzyme subunit McrC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/39);  total (36/39)

AGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTT  >  NC_000913/4577722‑4577990                                                                                                                                     |                                                                                                                                         aGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTa                                                                                                                                   <  1:387912/140‑1 (MQ=255) aGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTa                                                                                                                                   <  2:625814/140‑1 (MQ=255)     aaCTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGc                                                                                                                               >  1:10484/1‑140 (MQ=255)       ctctATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCtt                                                                                                                             <  1:320108/140‑1 (MQ=255)                 ccTTTGATGCCAGGAATGATCTCGGTGTTAGGCTTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTACTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGc                                                                                                                   >  2:829424/1‑140 (MQ=255)                 ccTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCtt                                                                                                                             <  1:281329/130‑1 (MQ=255)                 ccTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCtt                                                                                                                             >  2:281329/1‑130 (MQ=255)                 ccTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGc                                                                                                                   <  1:823612/140‑1 (MQ=255)                  cTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGcc                                                                                                                  >  2:937052/1‑140 (MQ=255)                     tGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGc                                                                                                               <  2:146807/140‑1 (MQ=255)                          ccAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAAt                                                                                                          <  2:229634/140‑1 (MQ=255)                             ggAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTc                                                                                                       >  1:353959/1‑140 (MQ=255)                               aaTGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCt                                                                                                     >  2:878898/1‑140 (MQ=255)                               aaTGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCt                                                                                                     <  2:539355/140‑1 (MQ=255)                                   aTCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGACTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTaa                                                                                                 >  1:908569/1‑140 (MQ=255)                                    tctcGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGACTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTaaa                                                                                                <  2:908569/140‑1 (MQ=255)                                    tctcGGTGATAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTaaa                                                                                                <  1:43902/140‑1 (MQ=255)                                      tcGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATa                                                                                              <  2:387861/140‑1 (MQ=255)                                      tcGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATa                                                                                              <  2:366851/140‑1 (MQ=255)                                        ggTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAAc                                                                                            >  2:417190/1‑140 (MQ=255)                                         gtgtTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAAcc                                                                                           <  1:465055/140‑1 (MQ=255)                                         gtgtTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAAcc                                                                                           <  1:430529/140‑1 (MQ=255)                                         gtgtTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAAcc                                                                                           <  2:92996/140‑1 (MQ=255)                                         gtgtTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAAcc                                                                                           <  1:937052/140‑1 (MQ=255)                                           gtTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGACATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTACTTAATTTCCTGTAAATAAcccc                                                                                         <  1:829424/140‑1 (MQ=255)                                              aGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATg                                                                                      >  1:418445/1‑140 (MQ=255)                                                 aTTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCAt                                                                                   <  2:799119/140‑1 (MQ=255)                                                 aTTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCAt                                                                                   <  2:944806/140‑1 (MQ=255)                                                   tGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAg                                                                                 >  2:158751/1‑140 (MQ=255)                                                    gTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAgg                                                                                >  1:96764/1‑140 (MQ=255)                                                     tAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGt                                                                               <  1:801958/140‑1 (MQ=255)                                                       aTCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTaa                                                                             <  2:128885/140‑1 (MQ=255)                                                         caagctcaagcCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGc                                                                           <  2:325454/140‑1 (MQ=255)                                                          aagctcaagcCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGc                                                                           >  2:845506/1‑139 (MQ=255)                                                          aagctcaagcCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGc                                                                           <  1:845506/139‑1 (MQ=255)                                                          aagctcaagcCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCa                                                                          >  2:673730/1‑140 (MQ=255)                                                           agctcaagcCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGc                                                                                     >  2:78369/1‑128 (MQ=255)                                                           agctcaagcCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGc                                                                                     <  1:78369/128‑1 (MQ=255)                                                           agctcaagcCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCAt                                                                         <  1:381392/140‑1 (MQ=255)                                                           agctcaagcCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCAt                                                                         >  2:542234/1‑140 (MQ=255)                                                             ctcaagcCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCTTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGt                                                                       >  1:409418/1‑140 (MQ=255)                                                               caagcCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTaa                                                                     >  1:524132/1‑140 (MQ=255)                                                                  gcCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAg                                                                  <  1:673730/140‑1 (MQ=255)                                                                  gcCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCAGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAg                                                                  <  2:735199/140‑1 (MQ=255)                                                                        ggCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTa                                                            >  1:608708/1‑140 (MQ=255)                                                                          cGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACg                                                          <  1:417190/140‑1 (MQ=255)                                                                                         aCCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACggg                                           <  2:524132/140‑1 (MQ=255)                                                                                         aCCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACggg                                           <  2:608708/140‑1 (MQ=255)                                                                                         aCCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACggg                                           <  2:10484/140‑1 (MQ=255)                                                                                          ccccTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGc                                          >  1:580330/1‑140 (MQ=255)                                                                                           cccTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCt                                         <  2:580330/140‑1 (MQ=255)                                                                                             cTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGt                                       <  1:878898/140‑1 (MQ=255)                                                                                              tttatttaATACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCTTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGtt                                      <  2:409418/140‑1 (MQ=255)                                                                                                     aaTACATACCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTAt                               <  2:353959/140‑1 (MQ=255)                                                                                                        acatacCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAg                            <  1:339428/140‑1 (MQ=255)                                                                                                          atacCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTc                          >  1:498052/1‑140 (MQ=255)                                                                                                          atacCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTc                          >  2:251301/1‑140 (MQ=255)                                                                                                           tacCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTcc                         >  2:4598/1‑140 (MQ=255)                                                                                                            acCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                        >  2:835945/1‑140 (MQ=255)                                                                                                            acCCCAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTccc                        >  2:892153/1‑140 (MQ=255)                                                                                                                cAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAg                            <  1:101797/132‑1 (MQ=255)                                                                                                                cAGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAg                            >  2:101797/1‑132 (MQ=255)                                                                                                                 aGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCtaat                   >  2:385333/1‑140 (MQ=255)                                                                                                                 aGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCtaat                   >  1:652987/1‑140 (MQ=255)                                                                                                                 aGGATATCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGAATCACGGGCTGTTCCACTATGAGTCCCCtaat                   >  2:568037/1‑140 (MQ=255)                                                                                                                    atatCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTaataat                >  2:772405/1‑140 (MQ=255)                                                                                                                    atatCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTaataat                >  2:915745/1‑140 (MQ=255)                                                                                                                       tCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATttg             <  2:340627/140‑1 (MQ=255)                                                                                                                       tCAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATttg             >  1:339562/1‑140 (MQ=255)                                                                                                                        cAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATttgt            >  2:758321/1‑140 (MQ=255)                                                                                                                        cAAGAAGATTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATttgt            <  1:4598/140‑1 (MQ=255)                                                                                                                           gaagaTTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTgg         >  2:811652/1‑140 (MQ=255)                                                                                                                           gaagaTTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTgg         >  2:260267/1‑140 (MQ=255)                                                                                                                            aagaTTGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGt        >  1:854585/1‑140 (MQ=255)                                                                                                                                 tGTTGCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGT‑CCCTAATAATTTGTTGGTCCAttt  >  2:126727/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                         AGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTT  >  NC_000913/4577722‑4577990

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