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% 239.5 / NA 75pseudogene (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 (32/43);  total (32/43)

GCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTT  >  NC_000913/3862362‑3862630
                                                                                                                                       |                                                                                                                                     
gCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                                                    <  1:1149391/139‑1 (MQ=255)
gCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAg                                                                                                                                    >  2:268848/1‑139 (MQ=255)
     tGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAAt                                                                                                                               <  2:1387713/139‑1 (MQ=255)
       ttGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCa                                                                                                                             <  1:1094904/139‑1 (MQ=255)
          aaCTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGTATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGAc                                                                                                                          <  2:16561/139‑1 (MQ=255)
          aaCTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGAc                                                                                                                          >  2:889282/1‑139 (MQ=255)
          aaCTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGAc                                                                                                                          <  1:634694/139‑1 (MQ=255)
              gCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCGATCAGACCTcc                                                                                                                      >  1:1083362/1‑139 (MQ=255)
               cAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                                                     >  2:370860/1‑139 (MQ=255)
               cAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                                                     <  1:712294/139‑1 (MQ=255)
                 gAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                   >  2:1599212/1‑139 (MQ=255)
                      aGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGtt                                                                                                             <  1:268848/139‑1 (MQ=255)
                         gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcac                                                                                                          <  1:1386920‑M1/139‑4 (MQ=255)
                           gCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcacca                                                                                                        <  1:370860‑M1/139‑6 (MQ=255)
                            cGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCATTCAGACCTCCGCCATGTTcaccac                                                                                                       >  1:924823‑M1/1‑133 (MQ=255)
                              gAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacac                                                                                                     <  2:967643‑M1/139‑9 (MQ=255)
                              gAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacac                                                                                                     <  1:1639876‑M1/139‑9 (MQ=255)
                              gAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATACCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacac                                                                                                     <  1:1187087‑M1/139‑9 (MQ=255)
                               aaaaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacc                                                                                                    <  2:924823‑M1/139‑10 (MQ=255)
                                 aaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCCGAACATCGGGATCCAGTTCTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaa                                                                                                  <  2:967435‑M1/139‑12 (MQ=255)
                                   ccACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaag                                                                                                <  2:1083362‑M1/139‑14 (MQ=255)
                                                   ggTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcg                                                                                >  2:332873‑M1/1‑110 (MQ=255)
                                                   ggTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcg                                                                                >  2:1521189‑M1/1‑110 (MQ=255)
                                                   ggTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaag                                                                                                >  2:1513651‑M1/1‑110 (MQ=255)
                                                   ggTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaag                                                                                                <  1:1513651‑M1/123‑14 (MQ=255)
                                                     tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacca                                                                                                   >  2:606009‑M1/1‑108 (MQ=255)
                                                     tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacca                                                                                                   <  1:606009‑M1/118‑11 (MQ=255)
                                                        aGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacatt                                                                           >  1:831369‑M1/1‑105 (MQ=255)
                                                          cacaACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattga                                                                         >  2:66549‑M1/1‑103 (MQ=255)
                                                             aaCCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggc                                                                      >  2:178548‑M1/1‑100 (MQ=255)
                                                               cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccg                                                                    >  1:220175‑M1/1‑98 (MQ=255)
                                                                 cAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTTCGGTAAAACTTAGTCAATCAGACGTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgcc                                                                  <  1:188656‑M1/139‑44 (MQ=255)
                                                                    tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                               <  1:66549‑M1/139‑47 (MQ=255)
                                                                    tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                               >  1:232369‑M1/1‑93 (MQ=255)
                                                                         aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGATCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                                                          <  2:1076574‑M1/139‑52 (MQ=255)
                                                                         aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                                                          >  1:499371‑M1/1‑88 (MQ=255)
                                                                         aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                                                          >  1:1641878‑M1/1‑88 (MQ=255)
                                                                           ttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgt                                                        <  1:1521189‑M1/139‑54 (MQ=255)
                                                                              ggTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgta                                                     <  1:889282‑M1/139‑57 (MQ=255)
                                                                                         tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                                                     <  1:1385370/65‑1 (MQ=255)
                                                                                         tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCg                                                                                                                     >  2:1385370/1‑65 (MQ=255)
                                                                                              gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                               <  2:680649‑M1/113‑47 (MQ=255)
                                                                                              gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                               >  1:680649‑M1/1‑67 (MQ=255)
                                                                                               cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagc                                          <  2:773832‑M1/133‑68 (MQ=255)
                                                                                               cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagc                                          >  1:773832‑M1/1‑66 (MQ=255)
                                                                                                 tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCa                                                                                                                             >  1:263720/1‑49 (MQ=255)
                                                                                                 tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCa                                                                                                                             <  2:263720/49‑1 (MQ=255)
                                                                                                 tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagca                                         >  2:128160‑M1/1‑64 (MQ=255)
                                                                                                 tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagca                                         <  1:128160‑M1/132‑69 (MQ=255)
                                                                                                 tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagaga                                                                                             >  1:1439750‑M1/1‑64 (MQ=255)
                                                                                                 tGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagaga                                                                                             <  2:1439750‑M1/80‑17 (MQ=255)
                                                                                                  ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa                                 <  2:1502891‑M1/139‑77 (MQ=255)
                                                                                                  ggTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaa                                 <  2:499371‑M1/139‑77 (MQ=255)
                                                                                                   gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAAGTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggacgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaa                                >  2:329750‑M1/1‑62 (MQ=255)
                                                                                                   gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCACCATGTTcaccacaccaaagagatacattacggtcgacattgaggccg                                                                    <  2:1347177‑M1/103‑42 (MQ=255)
                                                                                                   gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCACCATGTTcaccacaccaaagagatacattacggtcgacattgaggccg                                                                    >  1:1347177‑M1/1‑62 (MQ=255)
                                                                                                       gctgacatcGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaa                            <  1:706663‑M1/135‑82 (MQ=255)
                                                                                                            cATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacagga                       <  2:1641878‑M1/139‑87 (MQ=255)
                                                                                                            cATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                                                            >  1:127585‑M1/1‑53 (MQ=255)
                                                                                                            cATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                                                            <  2:127585‑M1/102‑50 (MQ=255)
                                                                                                              tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaag                     >  2:342714‑M1/1‑51 (MQ=255)
                                                                                                              tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacag                                           <  1:1385089‑M1/117‑67 (MQ=255)
                                                                                                              tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacag                                           >  2:1385089‑M1/1‑51 (MQ=255)
                                                                                                                 ggATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggta                  <  1:332873‑M1/139‑92 (MQ=255)
                                                                                                                  gATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGttcacc                                                                                                         >  2:602885/1‑47 (MQ=255)
                                                                                                                  gATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGttcacc                                                                                                         <  1:602885/51‑5 (MQ=255)
                                                                                                                    tCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTT‑‑‑CACCACACCAAAgaga                                                                                          >  2:210877/1‑61 (MQ=255)
                                                                                                                    tCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTT‑‑‑CACCACACCAAAgaga                                                                                          <  1:210877/61‑1 (MQ=255)
                                                                                                                         ttttGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccag          <  2:373730‑M1/139‑100 (MQ=255)
                                                                                                                         ttttGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccag          <  2:691879‑M1/139‑100 (MQ=255)
                                                                                                                         ttttGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattcc            <  1:62568‑M1/137‑98 (MQ=255)
                                                                                                                         ttttGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattcc            >  2:62568‑M1/1‑40 (MQ=255)
                                                                                                                                gTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcg   >  1:168414‑M1/1‑33 (MQ=255)
                                                                                                                                 tAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcgg  <  2:831369‑M1/139‑108 (MQ=255)
                                                                                                                                 tAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcgg  <  1:1599212‑M1/139‑108 (MQ=255)
                                                                                                                                       |                                                                                                                                     
GCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTT  >  NC_000913/3862362‑3862630

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