Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,210,436	T→C	W31R (TGG→CGG)	mcrA →	e14 prophage; 5‑methylcytosine‑specific restriction enzyme McrA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,210,436	0	T	C	100.0%	230.6 / NA	65	W31R (TGG→CGG)	mcrA	e14 prophage; 5‑methylcytosine‑specific restriction enzyme McrA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/31);  total (34/31)

TAAAAGTAGTATTGTCGTGAAAAATTGATTAAAGATTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCAT  >  NC_000913/1210303‑1210564                                                                                                                                      |                                                                                                                                 tAAAAGTAGTATTGTCGTGAAAAATTGATTAAAGATTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGc                                                                                                                            <  2:39606/140‑1 (MQ=255)               tCGTGAAAAATTGATTAAAGATTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAAc                                                                                                              >  2:383603/1‑140 (MQ=255)               tCGTGAAAAATTGATTAAAGATTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAAc                                                                                                              <  1:436071/140‑1 (MQ=255)                 gTGAAAAATTGATTAAAGATTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACaa                                                                                                            >  1:471929/1‑140 (MQ=255)                    aaaaaTTGATTAAAGATTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGa                                                                                                         >  2:299195/1‑140 (MQ=255)                     aaaaTTGATTAAAGATTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGAt                                                                                                        <  2:471929/140‑1 (MQ=255)                          tGATTAAAGATTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACaa                                                                                                   <  2:689161/140‑1 (MQ=255)                               aaaGATTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCg                                                                                              <  1:631316/140‑1 (MQ=255)                                aaGATTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGc                                                                                             <  1:41356/140‑1 (MQ=255)                                 aGATTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGct                                                                                            >  2:703699/1‑140 (MQ=255)                                  gATTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGctc                                                                                           <  2:496000/140‑1 (MQ=255)                                   aTTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGctct                                                                                          <  1:564333/140‑1 (MQ=255)                                      aaTATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGa                                                                                       <  1:703699/140‑1 (MQ=255)                                          ttatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGct                                                                                            <  1:39465/131‑1 (MQ=255)                                          ttatGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGct                                                                                            >  2:39465/1‑131 (MQ=255)                                             tGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATc                                                                                >  2:476457/1‑140 (MQ=255)                                             tGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATc                                                                                >  2:85458/1‑140 (MQ=255)                                             tGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATc                                                                                >  2:16527/1‑140 (MQ=255)                                                  gTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtat                                                                                   >  1:39398/1‑132 (MQ=255)                                                  gTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATtat                                                                                   <  2:39398/132‑1 (MQ=255)                                                  gTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGa                                                                           <  1:16527/140‑1 (MQ=255)                                                  gTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGa                                                                           >  2:361691/1‑140 (MQ=255)                                                    ttttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGctc                                                                                           <  2:489221/122‑1 (MQ=255)                                                    ttttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGctc                                                                                           >  1:489221/1‑122 (MQ=255)                                                     tttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACg                                                                        >  2:593907/1‑140 (MQ=255)                                                      ttGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACgg                                                                       <  1:559114/140‑1 (MQ=255)                                                       tGATAATAATGGAATTGGACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACATAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggt                                                                      >  2:73809/1‑140 (MQ=255)                                                       tGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggt                                                                      >  2:138022/1‑140 (MQ=255)                                                           aataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACATAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggt                                                                  <  1:73809/140‑1 (MQ=255)                                                           aataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggt                                                                  <  1:631186/140‑1 (MQ=255)                                                           aataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggt                                                                  <  1:138022/140‑1 (MQ=255)                                                           aataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggt                                                                  <  1:383603/140‑1 (MQ=255)                                                           aataatGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACggttggt                                                                  <  1:361691/140‑1 (MQ=255)                                                                tGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGAtt                                                             <  1:476457/140‑1 (MQ=255)                                                                   aaTTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTg                                                          >  1:1071/1‑140 (MQ=255)                                                                      tGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGgtg                                                       >  1:585817/1‑140 (MQ=255)                                                                        aaCTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTa                                                     >  1:1713/1‑140 (MQ=255)                                                                        aaCTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTa                                                     >  1:561371/1‑140 (MQ=255)                                                                            gAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATccc                                                 <  2:1071/140‑1 (MQ=255)                                                                               aGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAg                                              <  2:561371/140‑1 (MQ=255)                                                                               aGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAg                                              <  1:49722/140‑1 (MQ=255)                                                                                gCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGt                                             <  2:606064/140‑1 (MQ=255)                                                                                  tGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCg                                           >  2:703126/1‑140 (MQ=255)                                                                                       gtTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCgtagta                                      >  1:136482/1‑140 (MQ=255)                                                                                              aGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTgg                               >  2:236134/1‑140 (MQ=255)                                                                                                  gAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGtc                           >  2:384915/1‑140 (MQ=255)                                                                                                   aaGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGtca                          >  1:189490/1‑140 (MQ=255)                                                                                                          tGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGtcatcatc                     >  1:509644/1‑138 (MQ=255)                                                                                                          tGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGtcatcatc                     <  2:509644/138‑1 (MQ=255)                                                                                                          tGGTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCAtcag                   >  1:727201/1‑140 (MQ=255)                                                                                                           ggTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCAtcagt                  >  1:530125/1‑140 (MQ=255)                                                                                                           ggTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCAtcagt                  >  1:637930/1‑140 (MQ=255)                                                                                                           ggTGTTTATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCAtcagt                  <  2:637930/140‑1 (MQ=255)                                                                                                               tttATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGa              <  2:727201/140‑1 (MQ=255)                                                                                                                ttATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGaa             >  2:171535/1‑140 (MQ=255)                                                                                                                ttATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGaa             <  2:1713/140‑1 (MQ=255)                                                                                                                ttATGGTCTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGaa             >  2:443970/1‑140 (MQ=255)                                                                                                                       cTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACatat      >  2:553962/1‑140 (MQ=255)                                                                                                                       cTAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACatat      >  2:592610/1‑140 (MQ=255)                                                                                                                        tAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCAtcagt                  <  2:242094/127‑1 (MQ=255)                                                                                                                        tAATCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCAtcagt                  >  1:242094/1‑127 (MQ=255)                                                                                                                          aTCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCa   <  2:189490/140‑1 (MQ=255)                                                                                                                           tCCTTGAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCAt  >  2:512991/1‑140 (MQ=255)                                                                                                                                gAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCa   <  2:520842/134‑1 (MQ=255)                                                                                                                                gAGTCGCGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCa   >  1:520842/1‑134 (MQ=255)                                                                                                                                      |                                                                                                                                 TAAAAGTAGTATTGTCGTGAAAAATTGATTAAAGATTAATATTATGCATGTTTTTGATAATAATGGAATTGAACTGAAAGCTGAGTGTTCGATAGGTGAAGAGGATGGTGTTTATGGTCTAATCCTTGAGTCGTGGGGGCCGGGTGACAGAAACAAAGATTACAATATCGCTCTTGATTATATCATTGAACGGTTGGTTGATTCTGGTGTATCCCAAGTCGTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCAT  >  NC_000913/1210303‑1210564

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