Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,862,518 Δ1 bp pseudogene (1086/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,5160C.100.0% 294.1 / NA 69pseudogene (1088/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 . (32/37);  total (32/37)

GAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGC  >  NC_000913/3862379‑3862647
                                                                                                                                         |                                                                                                                                   
gAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGGAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                                    >  2:828997/1‑139 (MQ=255)
       tgtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTca                                                                                                                            >  1:1124954‑M1/1‑137 (MQ=255)
        gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                                    <  2:929543/131‑1 (MQ=255)
        gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                                    >  1:929543/1‑131 (MQ=255)
        gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcac                                                                                                                           >  2:134176‑M1/1‑136 (MQ=255)
           cGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccac                                                                                                                        >  1:755755‑M1/1‑133 (MQ=255)
           cGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccac                                                                                                                        <  1:672765‑M1/139‑7 (MQ=255)
              aaaaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacacc                                                                                                                     <  1:134176‑M1/139‑10 (MQ=255)
               aaaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacacca                                                                                                                    <  2:259256‑M1/139‑11 (MQ=255)
               aaaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacacca                                                                                                                    >  1:91231‑M1/1‑129 (MQ=255)
                  ccACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaag                                                                                                                 <  2:91231‑M1/139‑14 (MQ=255)
                   cACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccaca                                                                                                                       <  1:954566‑M1/132‑8 (MQ=255)
                   cACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccaca                                                                                                                       >  2:954566‑M1/1‑125 (MQ=255)
                         gTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagataca                                                                                                          <  1:940754‑M1/139‑21 (MQ=255)
                               cagGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcacca                                                                                                                         <  2:1038616‑M1/118‑6 (MQ=255)
                               cagGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcacca                                                                                                                         >  1:1038616‑M1/1‑113 (MQ=255)
                                  ggTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcg                                                                                                 >  2:553170‑M1/1‑110 (MQ=255)
                                  ggTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcg                                                                                                 >  1:978526‑M1/1‑110 (MQ=255)
                                    tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgac                                                                                               >  2:684756‑M1/1‑108 (MQ=255)
                                    tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacacc                                                                                                                     >  2:789829‑M1/1‑108 (MQ=255)
                                    tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacacc                                                                                                                     <  1:789829‑M1/117‑10 (MQ=255)
                                    tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccac                                                                                                                        >  2:33540‑M1/1‑108 (MQ=255)
                                    tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccac                                                                                                                        <  1:33540‑M1/114‑7 (MQ=255)
                                         cacaACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagaga                                                                                                              >  2:995655‑M1/1‑103 (MQ=255)
                                         cacaACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagaga                                                                                                              <  1:995655‑M1/119‑17 (MQ=255)
                                          acaACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgag                                                                                         <  2:978526‑M1/139‑38 (MQ=255)
                                            aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTTCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaa                                                                                                                   <  2:135177‑M1/111‑12 (MQ=255)
                                            aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggc                                                                                       >  2:1049232‑M1/1‑100 (MQ=255)
                                            aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaa                                                                                                                   >  1:135177‑M1/1‑100 (MQ=255)
                                              cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTGCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccg                                                                                     >  1:718657‑M1/1‑98 (MQ=255)
                                                   tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                                                <  2:574284‑M1/139‑47 (MQ=255)
                                                   tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgag                                                                                         <  1:860131‑M1/130‑38 (MQ=255)
                                                   tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgag                                                                                         >  2:860131‑M1/1‑93 (MQ=255)
                                                        aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                                                                           >  1:520844‑M1/1‑88 (MQ=255)
                                                        aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                                                                           >  2:775054‑M1/1‑88 (MQ=255)
                                                                   catcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaa                                                                >  1:192097‑M1/1‑77 (MQ=255)
                                                                       atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacag                                                            <  1:775054‑M1/139‑67 (MQ=255)
                                                                       atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacag                                                            <  2:572186‑M1/139‑67 (MQ=255)
                                                                       atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacag                                                            <  1:553170‑M1/139‑67 (MQ=255)
                                                                              cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCGTGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggt                                                     <  2:718657‑M1/139‑74 (MQ=255)
                                                                               gTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtg                                                    <  2:767701‑M1/139‑75 (MQ=255)
                                                                                 ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa                                                  >  1:1077840‑M1/1‑63 (MQ=255)
                                                                                 ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa                                                  <  2:520844‑M1/139‑77 (MQ=255)
                                                                                 ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa                                                  <  1:1049232‑M1/139‑77 (MQ=255)
                                                                                 ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa                                                  <  1:266610‑M1/139‑77 (MQ=255)
                                                                                  gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgca                                                                 >  1:24732‑M1/1‑62 (MQ=255)
                                                                                  gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgca                                                                 <  2:24732‑M1/123‑62 (MQ=255)
                                                                                    tGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaat                                                <  2:396905‑M1/138‑79 (MQ=255)
                                                                                    tGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaat                                                >  1:396905‑M1/1‑60 (MQ=255)
                                                                                    tGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattac                                                                                                      >  2:229015‑M1/1‑60 (MQ=255)
                                                                                    tGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattac                                                                                                      <  1:229015‑M1/84‑25 (MQ=255)
                                                                                     gCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgag                                                                                         >  1:25228‑M1/1‑59 (MQ=255)
                                                                                     gCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgag                                                                                         <  2:25228‑M1/96‑38 (MQ=255)
                                                                                             tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGttcacc                                                                                                                          <  2:783041/55‑5 (MQ=255)
                                                                                             tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGttcacc                                                                                                                          >  1:783041/1‑51 (MQ=255)
                                                                                             tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGtt                                                                                                                              >  1:706467/1‑51 (MQ=255)
                                                                                             tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGtt                                                                                                                              <  2:706467/51‑1 (MQ=255)
                                                                                             tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                                                                             <  1:491344‑M1/100‑50 (MQ=255)
                                                                                             tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                                                                             >  2:491344‑M1/1‑51 (MQ=255)
                                                                                             tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTT‑‑‑CACCACACCAAAGAGATACAtt                                                                                                     <  2:863462/73‑1 (MQ=255)
                                                                                             tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTT‑‑‑CACCACACCAAAGAGATACAtt                                                                                                     >  1:863462/1‑73 (MQ=255)
                                                                                              cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgca                                                                 >  1:1044877‑M1/1‑50 (MQ=255)
                                                                                              cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgca                                                                 <  2:1044877‑M1/111‑62 (MQ=255)
                                                                                                ggATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattggggccgccagcacggcgtgtaccgcaaacag                                                            <  2:513781‑M1/114‑67 (MQ=255)
                                                                                                ggATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattggggccgccagcacggcgtgtaccgcaaacag                                                            >  1:513781‑M1/1‑48 (MQ=255)
                                                                                                      aGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccatg                                                                                                                                 >  2:502423/1‑37 (MQ=255)
                                                                                                      aGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccatg                                                                                                                                 <  1:502423/40‑4 (MQ=255)
                                                                                                       gTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggt                                                                                                   >  1:154633‑M1/1‑41 (MQ=255)
                                                                                                       gTTTTGCGGTAAAACTTAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggt                                                                                                   <  2:154633‑M1/68‑28 (MQ=255)
                                                                                                                       tAGTCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaattcc            <  2:930474‑M1/139‑115 (MQ=255)
                                                                                                                          tCAATCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaattcccac         <  1:828997‑M1/139‑118 (MQ=255)
                                                                                                                              tCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa                                                  >  2:393778‑M1/1‑18 (MQ=255)
                                                                                                                              tCAGACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa                                                  <  1:393778‑M1/94‑77 (MQ=255)
                                                                                                                                 gACCTCCG‑CCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaattcccaccagcatg  <  2:192097‑M1/139‑125 (MQ=255)
                                                                                                                                         |                                                                                                                                   
GAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGC  >  NC_000913/3862379‑3862647

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