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% 280.7 / NA 85pseudogene (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/45);  total (40/45)

TGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCG  >  NC_000913/3862367‑3862627
                                                                                                                                  |                                                                                                                                  
tGTTGAACTGCGGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGACCCAGTTTTGCGGTAAAACTTAGTCAAt                                                                                                                            <  1:1096219/139‑1 (MQ=255)
           aGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                                                  <  2:1439747/138‑1 (MQ=255)
           aGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                                                  >  1:1439747/1‑138 (MQ=255)
           aGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGc                                                                                                                 >  1:694445/1‑139 (MQ=255)
           aGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGc                                                                                                                 >  1:1592946/1‑139 (MQ=255)
            gAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                >  2:334215/1‑139 (MQ=255)
                  gtgtGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGCCAATCAGACCt                                                                                                                     <  1:37747/128‑1 (MQ=255)
                  gtgtGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGCCAATCAGACCt                                                                                                                     >  2:37747/1‑128 (MQ=255)
                   tgtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTca                                                                                                        <  2:1592946‑M1/139‑3 (MQ=255)
                   tgtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTca                                                                                                        <  2:747762‑M1/139‑3 (MQ=255)
                   tgtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTca                                                                                                        <  1:334215‑M1/139‑3 (MQ=255)
                    gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccat                                                                                                              >  1:258179/1‑131 (MQ=255)
                    gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccat                                                                                                              <  2:258179/133‑3 (MQ=255)
                       cGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccac                                                                                                    >  1:359479‑M1/1‑133 (MQ=255)
                         gAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                                                  <  2:1040508/124‑1 (MQ=255)
                         aaaaaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTGCg                                                                                                                  >  1:1040508/2‑124 (MQ=255)
                          aaaaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacc                                                                                                 <  2:791901‑M1/139‑10 (MQ=255)
                            aaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcac                                                                                                       >  2:1189681‑M1/1‑128 (MQ=255)
                            aaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcac                                                                                                       <  1:1189681‑M1/131‑4 (MQ=255)
                             aCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaa                                                                                              >  1:426890‑M1/1‑127 (MQ=255)
                             aCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGCCAATCAGACCTCCGCCATGTTcaccacaccaaa                                                                                              <  1:184978‑M1/139‑13 (MQ=255)
                                   aaGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagata                                                                                        <  2:275333‑M1/139‑19 (MQ=255)
                                    aGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatac                                                                                       >  2:817082‑M1/1‑120 (MQ=255)
                                              ggTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcg                                                                             >  2:1142557‑M1/1‑110 (MQ=255)
                                                tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgac                                                                           >  2:942130‑M1/1‑108 (MQ=255)
                                                   aGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagaga                                                                                          <  2:769255‑M1/121‑17 (MQ=255)
                                                   aGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagaga                                                                                          >  1:769255‑M1/1‑105 (MQ=255)
                                                    gCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattg                                                                       <  2:1680142‑M1/139‑36 (MQ=255)
                                                        aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCa                                                                                                                          >  2:1475696/1‑85 (MQ=255)
                                                        aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCa                                                                                                                          <  1:1475696/85‑1 (MQ=255)
                                                        aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggc                                                                   >  1:535187‑M1/1‑100 (MQ=255)
                                                         aCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccaca                                                                                                   >  2:1630671‑M1/1‑99 (MQ=255)
                                                         aCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccaca                                                                                                   <  1:1630671‑M1/106‑8 (MQ=255)
                                                                    aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAAt                                                                                                                            >  2:1374576/1‑71 (MQ=255)
                                                                    aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAAt                                                                                                                            <  1:1374576/71‑1 (MQ=255)
                                                                    aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCt                                                                                                                     <  1:734952/78‑1 (MQ=255)
                                                                    aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCt                                                                                                                     >  2:734952/1‑78 (MQ=255)
                                                                    aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                                                       >  2:929491‑M1/1‑88 (MQ=255)
                                                                    aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                                                       <  2:1359113‑M1/139‑52 (MQ=255)
                                                                    aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccaca                                                                                                   <  2:73323‑M1/95‑8 (MQ=255)
                                                                    aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccaca                                                                                                   >  1:73323‑M1/1‑88 (MQ=255)
                                                                               catcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaa                                            >  2:688852‑M1/1‑77 (MQ=255)
                                                                                 tcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaac                                          >  1:622184‑M1/1‑75 (MQ=255)
                                                                                   atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacag                                        <  1:1142557‑M1/139‑67 (MQ=255)
                                                                                     atCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagca                                      <  1:680528‑M1/139‑69 (MQ=255)
                                                                                      tCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaa                                     <  2:694445‑M1/139‑70 (MQ=255)
                                                                                        cgcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacg                                   <  1:1687827‑M1/139‑72 (MQ=255)
                                                                                         gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaag                                                                                             <  2:1189283‑M1/80‑14 (MQ=255)
                                                                                         gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaag                                                                                             >  1:1189283‑M1/1‑67 (MQ=255)
                                                                                          cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagca                                      >  2:1454846‑M1/1‑66 (MQ=255)
                                                                                          cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagca                                      <  1:1454846‑M1/134‑69 (MQ=255)
                                                                                          cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagca                                                           <  2:1401390‑M1/113‑48 (MQ=255)
                                                                                          cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagca                                                           >  1:1401390‑M1/1‑66 (MQ=255)
                                                                                          cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                            <  2:1421251‑M1/112‑47 (MQ=255)
                                                                                          cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                            >  1:1421251‑M1/1‑66 (MQ=255)
                                                                                          cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagaga                                                                                          <  1:982221‑M1/82‑17 (MQ=255)
                                                                                          cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagaga                                                                                          >  2:982221‑M1/1‑66 (MQ=255)
                                                                                          cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTca                                                                                                        >  2:511996‑M1/1‑66 (MQ=255)
                                                                                          cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTca                                                                                                        <  1:511996‑M1/68‑3 (MQ=255)
                                                                                           gTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtg                                <  2:426890‑M1/139‑75 (MQ=255)
                                                                                            tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                                                 >  1:335159/1‑42 (MQ=255)
                                                                                            tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                                                 <  2:335159/42‑1 (MQ=255)
                                                                                             ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccatg                                                                                                             >  2:340758/1‑58 (MQ=255)
                                                                                             ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccatg                                                                                                             <  1:340758/61‑4 (MQ=255)
                                                                                                  cTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaa                         <  1:942130‑M1/139‑82 (MQ=255)
                                                                                                  cTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaa                         <  2:535187‑M1/139‑82 (MQ=255)
                                                                                                      aCATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacagg                     <  2:598081‑M1/139‑86 (MQ=255)
                                                                                                          cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagca                                      >  2:135419‑M1/1‑50 (MQ=255)
                                                                                                          cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagca                                      <  1:135419‑M1/118‑69 (MQ=255)
                                                                                                            ggATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTT‑‑‑CACCACACCaa                                                                                            >  2:1368252/1‑59 (MQ=255)
                                                                                                            ggATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTT‑‑‑CACCACACCaa                                                                                            <  1:1368252/59‑1 (MQ=255)
                                                                                                               tCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTT‑‑‑CACCACACCAAag                                                                                          >  1:656608/1‑58 (MQ=255)
                                                                                                               tCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTT‑‑‑CACCACACCAAag                                                                                          <  2:656608/58‑1 (MQ=255)
                                                                                                                 cAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacagg                     <  1:22811‑M1/128‑86 (MQ=255)
                                                                                                                 cAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacagg                     >  2:22811‑M1/1‑43 (MQ=255)
                                                                                                                    ttttGCGGCAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaa                                     >  2:111684‑M1/1‑40 (MQ=255)
                                                                                                                    ttttGCGGCAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaa                                     <  1:111684‑M1/109‑70 (MQ=255)
                                                                                                                       tGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcgg    >  2:933804‑M1/1‑37 (MQ=255)
                                                                                                                         cGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggtt  <  1:817082‑M1/139‑105 (MQ=255)
                                                                                                                         cGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagc                                       >  2:1077292‑M1/1‑35 (MQ=255)
                                                                                                                         cGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagc                                       <  1:1077292‑M1/102‑68 (MQ=255)
                                                                                                                         cGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                                                         <  2:844380‑M1/84‑50 (MQ=255)
                                                                                                                         cGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                                                         >  1:844380‑M1/1‑35 (MQ=255)
                                                                                                                         cGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                            >  2:48933‑M1/1‑35 (MQ=255)
                                                                                                                         cGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                            <  1:48933‑M1/81‑47 (MQ=255)
                                                                                                                                  |                                                                                                                                  
TGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCG  >  NC_000913/3862367‑3862627

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