Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,862,497 C→T pseudogene (1107/1617 nt) glvC ← pseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,4970CT100.0% 258.4 / NA 79pseudogene (1107/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/39);  total (40/39)

TTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGA  >  NC_000913/3862369‑3862601
                                                                                                                                |                                                                                                        
ttGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCa                                                                                                >  1:259256/1‑139 (MQ=255)
   aaCGGCAGAATGAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGAc                                                                                             <  1:734632/139‑1 (MQ=255)
        cAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                        >  1:572186/1‑139 (MQ=255)
        cAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                        >  1:1104907/1‑139 (MQ=255)
          gAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGGAAAACTTAGTCAATCAGACCTCCGcc                                                                                      >  2:828997/1‑139 (MQ=255)
                 tgtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTca                                                                              >  1:1124954‑M1/1‑137 (MQ=255)
                  gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                      <  2:929543/131‑1 (MQ=255)
                  gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                      >  1:929543/1‑131 (MQ=255)
                  gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcac                                                                             >  2:134176‑M1/1‑136 (MQ=255)
                     cGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccac                                                                          >  1:755755‑M1/1‑133 (MQ=255)
                     cGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccac                                                                          <  1:672765‑M1/139‑7 (MQ=255)
                        aaaaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacc                                                                       <  1:134176‑M1/139‑10 (MQ=255)
                         aaaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacca                                                                      >  1:91231‑M1/1‑129 (MQ=255)
                         aaaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacca                                                                      <  2:259256‑M1/139‑11 (MQ=255)
                            ccACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaag                                                                   <  2:91231‑M1/139‑14 (MQ=255)
                             cACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccaca                                                                         >  2:954566‑M1/1‑125 (MQ=255)
                             cACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccaca                                                                         <  1:954566‑M1/132‑8 (MQ=255)
                                   gTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagataca                                                            <  1:940754‑M1/139‑21 (MQ=255)
                                         cagGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcacca                                                                           >  1:1038616‑M1/1‑113 (MQ=255)
                                         cagGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcacca                                                                           <  2:1038616‑M1/118‑6 (MQ=255)
                                            ggTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcg                                                   >  1:978526‑M1/1‑110 (MQ=255)
                                            ggTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcg                                                   >  2:553170‑M1/1‑110 (MQ=255)
                                              tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgac                                                 >  2:684756‑M1/1‑108 (MQ=255)
                                              tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacc                                                                       >  2:789829‑M1/1‑108 (MQ=255)
                                              tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacc                                                                       <  1:789829‑M1/117‑10 (MQ=255)
                                              tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccac                                                                          <  1:33540‑M1/114‑7 (MQ=255)
                                              tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccac                                                                          >  2:33540‑M1/1‑108 (MQ=255)
                                                aaGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                        >  1:96516/1‑99 (MQ=255)
                                                aaGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                        <  2:96516/99‑1 (MQ=255)
                                                   cacaACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagaga                                                                <  1:995655‑M1/119‑17 (MQ=255)
                                                   cacaACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagaga                                                                >  2:995655‑M1/1‑103 (MQ=255)
                                                    acaACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgag                                           <  2:978526‑M1/139‑38 (MQ=255)
                                                      aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTTCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaa                                                                     <  2:135177‑M1/111‑12 (MQ=255)
                                                      aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggc                                         >  2:1049232‑M1/1‑100 (MQ=255)
                                                      aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaa                                                                     >  1:135177‑M1/1‑100 (MQ=255)
                                                        cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTGCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccg                                       >  1:718657‑M1/1‑98 (MQ=255)
                                                           aaTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                       <  2:295808/73‑1 (MQ=255)
                                                           aaTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                       >  1:295808/1‑73 (MQ=255)
                                                             tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                  <  2:574284‑M1/139‑47 (MQ=255)
                                                             tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgag                                           >  2:860131‑M1/1‑93 (MQ=255)
                                                             tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgag                                           <  1:860131‑M1/130‑38 (MQ=255)
                                                                  aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                             >  1:520844‑M1/1‑88 (MQ=255)
                                                                  aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                             >  2:775054‑M1/1‑88 (MQ=255)
                                                                             catcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaa                  >  1:192097‑M1/1‑77 (MQ=255)
                                                                                 atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacag              <  1:553170‑M1/139‑67 (MQ=255)
                                                                                 atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacag              <  2:572186‑M1/139‑67 (MQ=255)
                                                                                 atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacag              <  1:775054‑M1/139‑67 (MQ=255)
                                                                                        cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCGTGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggt       <  2:718657‑M1/139‑74 (MQ=255)
                                                                                         gTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtg      <  2:767701‑M1/139‑75 (MQ=255)
                                                                                           ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa    <  1:1049232‑M1/139‑77 (MQ=255)
                                                                                           ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa    >  1:1077840‑M1/1‑63 (MQ=255)
                                                                                           ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa    <  2:520844‑M1/139‑77 (MQ=255)
                                                                                           ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa    <  1:266610‑M1/139‑77 (MQ=255)
                                                                                            gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgca                   >  1:24732‑M1/1‑62 (MQ=255)
                                                                                            gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgca                   <  2:24732‑M1/123‑62 (MQ=255)
                                                                                              tGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaat  <  2:396905‑M1/138‑79 (MQ=255)
                                                                                              tGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaat  >  1:396905‑M1/1‑60 (MQ=255)
                                                                                              tGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattac                                                        <  1:229015‑M1/84‑25 (MQ=255)
                                                                                              tGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattac                                                        >  2:229015‑M1/1‑60 (MQ=255)
                                                                                               gCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgag                                           >  1:25228‑M1/1‑59 (MQ=255)
                                                                                               gCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgag                                           <  2:25228‑M1/96‑38 (MQ=255)
                                                                                                       tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGttcacc                                                                            >  1:783041/1‑51 (MQ=255)
                                                                                                       tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGttcacc                                                                            <  2:783041/55‑5 (MQ=255)
                                                                                                       tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGtt                                                                                <  2:706467/51‑1 (MQ=255)
                                                                                                       tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGtt                                                                                >  1:706467/1‑51 (MQ=255)
                                                                                                       tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                               >  2:491344‑M1/1‑51 (MQ=255)
                                                                                                       tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                               <  1:491344‑M1/100‑50 (MQ=255)
                                                                                                       tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTT‑‑‑CACCACACCAAAGAGATACAtt                                                       <  2:863462/73‑1 (MQ=255)
                                                                                                       tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTT‑‑‑CACCACACCAAAGAGATACAtt                                                       >  1:863462/1‑73 (MQ=255)
                                                                                                        cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgca                   >  1:1044877‑M1/1‑50 (MQ=255)
                                                                                                        cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgca                   <  2:1044877‑M1/111‑62 (MQ=255)
                                                                                                          ggATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattggggccgccagcacggcgtgtaccgcaaacag              <  2:513781‑M1/114‑67 (MQ=255)
                                                                                                          ggATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattggggccgccagcacggcgtgtaccgcaaacag              >  1:513781‑M1/1‑48 (MQ=255)
                                                                                                             tCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                        <  2:465753/38‑1 (MQ=255)
                                                                                                             tCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                        >  1:465753/1‑38 (MQ=255)
                                                                                                                aGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccatg                                                                                   >  2:502423/1‑37 (MQ=255)
                                                                                                                aGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccatg                                                                                   <  1:502423/40‑4 (MQ=255)
                                                                                                                 gTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggt                                                     >  1:154633‑M1/1‑41 (MQ=255)
                                                                                                                 gTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggt                                                     <  2:154633‑M1/68‑28 (MQ=255)
                                                                                                                                |                                                                                                        
TTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGA  >  NC_000913/3862369‑3862601

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: