Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,577,958	C→T	M1M (GTG→ATG) † *460* (TAG→TAA)	mcrC ← mcrB ←	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,577,958	0	C	T	100.0%	262.5 / NA	80	M1M (GTG→ATG) † *460* (TAG→TAA)	mcrC mcrB	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (39/41);  total (39/41)

ATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGT  >  NC_000913/4577823‑4578086                                                                                                                                        |                                                                                                                                 atacatacCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTAt                                                                                                                               <  2:1408402/139‑1 (MQ=255)       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGa                                                                                                                             >  1:238773/1‑135 (MQ=255)       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGa                                                                                                                             <  2:238773/135‑1 (MQ=255)       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                         >  1:68207/1‑139 (MQ=255)       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                         >  2:707105/1‑139 (MQ=255)       acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                         >  2:445521/1‑139 (MQ=255)            aGGATATCAAGAAGCTTGTTACCGGGTATGGCTTCAAGGTATGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCtaa                                                                                                                    >  2:980121/1‑139 (MQ=255)                   cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttg                                                                                                             >  1:1489126/1‑139 (MQ=255)                   cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCACTAATAATttg                                                                                                             >  2:1046929/1‑139 (MQ=255)                         gaTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTc                                                                                                       <  1:707105/139‑1 (MQ=255)                            tGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCAt                                                                                                    >  1:576875/1‑139 (MQ=255)                            tGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCAt                                                                                                    <  2:576875/139‑1 (MQ=255)                              ttACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCAttt                                                                                                  <  2:1114030/139‑1 (MQ=255)                                aCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt                                                                                                <  2:126979/139‑1 (MQ=255)                                aCCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt                                                                                                <  1:802642/139‑1 (MQ=255)                                 ccGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                               <  2:68207/139‑1 (MQ=255)                                 ccGGGTATCGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                               <  1:980121/139‑1 (MQ=255)                                  cGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgt                                                                                              <  1:445521/139‑1 (MQ=255)                                   gggTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                             <  1:285210/139‑1 (MQ=255)                                   gggTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                             <  1:1046929/139‑1 (MQ=255)                                        tGGCTTTAAGGTTTGCCTGCTGAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGt                                                                                                                           <  1:11102/104‑1 (MQ=255)                                        tGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGt                                                                                                                           >  2:11102/1‑104 (MQ=255)                                          gCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTtata                                                                                      <  1:1300552/139‑1 (MQ=255)                                               aaGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGt                                                                                 >  1:700966/1‑139 (MQ=255)                                                  ttttGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcat                                                                              <  1:102370/138‑1 (MQ=255)                                                  gTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcat                                                                              <  2:1489126/139‑1 (MQ=255)                                                    ttGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                            <  1:1585258/139‑1 (MQ=255)                                                     tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg                                                                                                                            >  2:77380/1‑90 (MQ=255)                                                     tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg                                                                                                                            <  1:77380/90‑1 (MQ=255)                                                     tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTccc                                                                                                                        >  2:906482/1‑94 (MQ=255)                                                     tGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTccc                                                                                                                        <  1:906482/94‑1 (MQ=255)                                                      gCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg                                                                                                                            <  2:922751/89‑1 (MQ=255)                                                      gCCTGCTTAATTTCCTGAAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACAGGCTGTTCCATTATGAg                                                                                                                            >  1:922751/1‑89 (MQ=255)                                                              aaTTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgt                                                                                              >  1:1073906/1‑111 (MQ=255)                                                              aaTTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgt                                                                                              <  2:1073906/111‑1 (MQ=255)                                                                     tGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAg                                                           <  1:448134/139‑1 (MQ=255)                                                                                 cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGTGTCATCAAAGAAATAttct                                                                <  2:1096277/122‑1 (MQ=255)                                                                                 cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATAttct                                                                >  1:1096277/1‑122 (MQ=255)                                                                                  aTGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGata                                              >  2:810746/1‑139 (MQ=255)                                                                                   tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttgt                                                                                                            >  1:933314/1‑76 (MQ=255)                                                                                   tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttgt                                                                                                            <  2:933314/76‑1 (MQ=255)                                                                                         ggTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAgggg                                                   >  2:216649/1‑127 (MQ=255)                                                                                         ggTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAgggg                                                   <  1:216649/127‑1 (MQ=255)                                                                                         ggTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtca                                       >  2:805985/1‑139 (MQ=255)                                                                                               cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc                                                                                                                          <  2:661458/50‑1 (MQ=255)                                                                                               cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc                                                                                                                          >  1:661458/1‑50 (MQ=255)                                                                                               cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAatat                                                                   <  2:991274/105‑1 (MQ=255)                                                                                               cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAatat                                                                   >  1:991274/1‑105 (MQ=255)                                                                                               cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATAttctt                                                               >  2:357267/1‑109 (MQ=255)                                                                                               cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATAttctt                                                               <  1:357267/109‑1 (MQ=255)                                                                                               cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTaa                                                        >  1:466111/1‑116 (MQ=255)                                                                                               cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTaa                                                        <  2:466111/116‑1 (MQ=255)                                                                                                aTGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                >  2:762372/1‑139 (MQ=255)                                                                                                  gTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttgt                                                                                                            <  1:863676/61‑1 (MQ=255)                                                                                                  gTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttgt                                                                                                            >  2:863676/1‑61 (MQ=255)                                                                                                    aaTAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg                                                                                                                              <  2:83248/41‑1 (MQ=255)                                                                                                    aaTAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg                                                                                                                              >  1:83248/1‑41 (MQ=255)                                                                                                    aaTAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt                                                                                                >  2:1178963/1‑71 (MQ=255)                                                                                                    aaTAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt                                                                                                <  1:1178963/71‑1 (MQ=255)                                                                                                    aaTAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCa                            >  1:303366/1‑139 (MQ=255)                                                                                                         atatTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                             >  1:1009401/1‑69 (MQ=255)                                                                                                         atatTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                             <  2:1009401/69‑1 (MQ=255)                                                                                                           atTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc                                                                                                                          >  1:686006/1‑38 (MQ=255)                                                                                                           atTACGGACAGGTATCACGGGCTGTTCCATTATGAGTc                                                                                                                          <  2:686006/38‑1 (MQ=255)                                                                                                           atTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGcc                     <  2:700966/139‑1 (MQ=255)                                                                                                             tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTc                                                                                                       <  2:681924/55‑1 (MQ=255)                                                                                                             tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTc                                                                                                       >  1:681924/1‑55 (MQ=255)                                                                                                              aCGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCAtt                  >  2:1454735/1‑139 (MQ=255)                                                                                                                   cAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttgt                                                                                                            <  2:447514/44‑1 (MQ=255)                                                                                                                   cAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttgt                                                                                                            >  1:447514/1‑44 (MQ=255)                                                                                                                     ggTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATAttcttc                                                              <  1:360396/88‑1 (MQ=255)                                                                                                                     ggTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATAttcttc                                                              >  2:360396/1‑88 (MQ=255)                                                                                                                      gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATcc          >  2:1309164/1‑139 (MQ=255)                                                                                                                              ggCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGt  >  1:705000/1‑139 (MQ=255)                                                                                                                                cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                             <  2:30018/46‑1 (MQ=255)                                                                                                                                cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                             <  2:850196/46‑1 (MQ=255)                                                                                                                                cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                             >  1:850196/1‑46 (MQ=255)                                                                                                                                cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                             >  1:30018/1‑46 (MQ=255)                                                                                                                                cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                               >  1:1153500/1‑44 (MQ=255)                                                                                                                                cTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                               <  2:1153500/44‑1 (MQ=255)                                                                                                                                        |                                                                                                                                 ATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGT  >  NC_000913/4577823‑4578086

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