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_0009134,577,9580CT100.0% 317.1 / NA 98M1M (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 (53/45);  total (53/45)

ACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCAT  >  NC_000913/4577829‑4578091
                                                                                                                                 |                                                                                                                                     
acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCCTAGGTAAGCATGTAATAGATATTACGGCCAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                              >  2:396977/1‑139 (MQ=255)
acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                              >  1:23505/1‑139 (MQ=255)
acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                              >  1:363805/1‑139 (MQ=255)
acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                              >  1:1107201/1‑139 (MQ=255)
acCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTcc                                                                                                                              >  1:756148/1‑139 (MQ=255)
 ccccAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTccc                                                                                                                             <  1:1219646/139‑1 (MQ=255)
    cAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCta                                                                                                                          >  2:1004235/1‑139 (MQ=255)
            cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg                                                                                                                                 <  2:796125/124‑1 (MQ=255)
            cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttg                                                                                                                  >  1:1297511/1‑139 (MQ=255)
            cAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATCAg                                                                                                                                 >  1:796125/1‑124 (MQ=255)
                atgaTTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTgg                                                                                                              <  2:1004852/137‑1 (MQ=255)
                           cGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgt                                                                                                   <  1:803007/139‑1 (MQ=255)
                            gggTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                  <  1:1264072/139‑1 (MQ=255)
                                 tGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTta                                                                                             >  1:1297029/1‑139 (MQ=255)
                                       cAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                  >  1:899987/1‑128 (MQ=255)
                                       cAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtgtt                                                                                                  <  2:899987/128‑1 (MQ=255)
                                        aaGGTTTGCCTGCTTCATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGt                                                                                      >  2:866673/1‑139 (MQ=255)
                                        aaGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGt                                                                                      <  2:23505/139‑1 (MQ=255)
                                         aGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt                                                                                                     <  2:1258491/123‑1 (MQ=255)
                                         aGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCt                                                                                                     >  1:1258491/1‑123 (MQ=255)
                                            tttGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatc                                                                                  <  2:1297029/139‑1 (MQ=255)
                                             ttGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                                 <  1:1004235/139‑1 (MQ=255)
                                                    cTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg                                                                                                                                 >  2:639530/1‑84 (MQ=255)
                                                    cTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAg                                                                                                                                 <  1:639530/84‑1 (MQ=255)
                                                              tGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg                                                                                                                                   <  1:367349/72‑1 (MQ=255)
                                                              tGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg                                                                                                                                   >  2:367349/1‑72 (MQ=255)
                                                              tGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAg                                                                >  2:320165/1‑139 (MQ=255)
                                                              tGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAg                                                                >  2:649284/1‑139 (MQ=255)
                                                                 aaaTAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTaa                                                             <  1:182639/139‑1 (MQ=255)
                                                                 aaaTAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTaa                                                             <  1:649284/139‑1 (MQ=255)
                                                                       ccccATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg                                                                                                                                   >  2:418804/1‑63 (MQ=255)
                                                                       ccccATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATg                                                                                                                                   <  1:418804/63‑1 (MQ=255)
                                                                       ccccATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                    <  1:175355/94‑1 (MQ=255)
                                                                       ccccATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                    >  2:175355/1‑94 (MQ=255)
                                                                        cccATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCg                                                      <  2:363805/139‑1 (MQ=255)
                                                                          cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAgggg                                                        >  2:373453/1‑135 (MQ=255)
                                                                          cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAgggg                                                        <  1:373453/135‑1 (MQ=255)
                                                                          cATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGat                                                    >  1:422159/1‑139 (MQ=255)
                                                                            tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGatat                                                  >  1:377736/1‑139 (MQ=255)
                                                                            tGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGatat                                                  >  2:422117/1‑139 (MQ=255)
                                                                                      aGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                    <  2:1021156/79‑1 (MQ=255)
                                                                                      aGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                    >  1:1021156/1‑79 (MQ=255)
                                                                                       gCATGTAATAGATACTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATcaca                                       >  1:1168022/1‑139 (MQ=255)
                                                                                        cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGaat                                                                                                                                <  2:776825/49‑3 (MQ=255)
                                                                                        cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGt                                                                                                                                >  1:776825/1‑49 (MQ=255)
                                                                                        cATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACaa                                      >  2:506108/1‑139 (MQ=255)
                                                                                         aTGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                     <  2:1297511/139‑1 (MQ=255)
                                                                                         aTGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAAt                                     <  2:756148/139‑1 (MQ=255)
                                                                                               tAGATATTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt                               <  1:404753/139‑1 (MQ=255)
                                                                                                    atTCCGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGcc                          >  1:880439/1‑139 (MQ=255)
                                                                                                    atTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCtaataa                                                                                                                      >  2:1025802/1‑47 (MQ=255)
                                                                                                    atTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCtaataa                                                                                                                      <  1:1025802/47‑1 (MQ=255)
                                                                                                    atTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTta                                                                                             >  2:976599/1‑72 (MQ=255)
                                                                                                    atTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTta                                                                                             <  1:976599/72‑1 (MQ=255)
                                                                                                    atTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGcc                          >  1:820772/1‑139 (MQ=255)
                                                                                                    atTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGcc                          >  2:239875/1‑139 (MQ=255)
                                                                                                    atTACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGcc                          <  2:880439/139‑1 (MQ=255)
                                                                                                     ttACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCa                         >  1:645172/1‑139 (MQ=255)
                                                                                                      tACGGACAGGTATCACGGGCTGTTCCATTATGAGTccc                                                                                                                             >  1:624313/1‑38 (MQ=255)
                                                                                                      tACGGACAGGTATCACGGGCTGTTCCATTATGAGTccc                                                                                                                             <  2:624313/38‑1 (MQ=255)
                                                                                                      tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTta                                                                                             <  1:854573/70‑1 (MQ=255)
                                                                                                      tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt                               >  2:1060107/1‑132 (MQ=255)
                                                                                                      tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt                               <  1:1060107/132‑1 (MQ=255)
                                                                                                      tACGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTACTGTTGTTta                                                                                             >  2:854573/1‑70 (MQ=255)
                                                                                                       aCGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttgt                                                                                                                 <  2:1293431/49‑1 (MQ=255)
                                                                                                       aCGGACAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATttgt                                                                                                                 >  1:1293431/1‑49 (MQ=255)
                                                                                                           aCAGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGt                   >  1:231281/1‑139 (MQ=255)
                                                                                                             aGGTATCACGGGCTGTTCCATTATGAGTCCCTTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAATGGGTGATATCCGTCTTCACACTTTCATTAAGCCATTTCGTAt                 >  2:587419/1‑139 (MQ=255)
                                                                                                             aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                    <  1:777614/56‑1 (MQ=255)
                                                                                                             aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                    <  1:1245761/56‑1 (MQ=255)
                                                                                                             aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                    >  2:1245761/1‑56 (MQ=255)
                                                                                                             aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCtg                                                                                                    >  2:777614/1‑56 (MQ=255)
                                                                                                             aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                                 >  2:95890/1‑75 (MQ=255)
                                                                                                             aGGTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                                 <  1:95890/75‑1 (MQ=255)
                                                                                                               gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATc                <  1:512420/138‑1 (MQ=255)
                                                                                                               gTATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATc                >  2:512420/1‑138 (MQ=255)
                                                                                                                tATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                                 >  1:1034075/1‑72 (MQ=255)
                                                                                                                tATCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGtcatca                                                                                 <  2:1034075/72‑1 (MQ=255)
                                                                                                                  tCACGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGga            <  2:645172/139‑1 (MQ=255)
                                                                                                                    aCGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt                               <  1:1262550/118‑1 (MQ=255)
                                                                                                                    aCGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt                               >  2:1262550/1‑118 (MQ=255)
                                                                                                                    aCGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTg                      >  1:1059503/1‑127 (MQ=255)
                                                                                                                    aCGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTg                      <  2:1059503/127‑1 (MQ=255)
                                                                                                                     cGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt                               <  1:634389/117‑1 (MQ=255)
                                                                                                                     cGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCAtt                               >  2:634389/1‑117 (MQ=255)
                                                                                                                     cGGGCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGagag         <  1:320165/139‑1 (MQ=255)
                                                                                                                       ggCTGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGt       >  2:1235657/1‑139 (MQ=255)
                                                                                                                         cTGTTCCATTATGAGTCCCCTAATAATTTCTTGGTCCATTTCTGTAGTTTATGGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGc     >  2:747612/1‑139 (MQ=255)
                                                                                                                          tGTTCCATTATGAGTCCCCTAATAATTTGTTGGTc                                                                                                            <  1:1123325/35‑1 (MQ=255)
                                                                                                                          tGTTCCATTATGAGTCCCCTAATAATTTGTTGGTc                                                                                                            >  2:1123325/1‑35 (MQ=255)
                                                                                                                          tGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcat                                           >  1:779058/1‑100 (MQ=255)
                                                                                                                          tGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGtcat                                           <  2:779058/100‑1 (MQ=255)
                                                                                                                          tGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGcc    >  2:562201/1‑139 (MQ=255)
                                                                                                                          tGTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGcc    <  1:506108/139‑1 (MQ=255)
                                                                                                                           gTTCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCa   <  2:820772/139‑1 (MQ=255)
                                                                                                                            ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCAttt                                                                                                       >  1:1027654/1‑38 (MQ=255)
                                                                                                                            ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCAttt                                                                                                       <  2:1027654/38‑1 (MQ=255)
                                                                                                                            ttCCATTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTCGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGAAGAGGTGCCAt  >  1:815632/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                     
ACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCAT  >  NC_000913/4577829‑4578091

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