Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,198,736 Δ1 bp 100% coding (133/729 nt) cobA2 ← uroporphyrinogen‑III C‑methyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,198,7360C→.100.0% 617.8 142coding (133/729 nt)cobA2uroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  new base (74/66):  ref base (0/1):  total (74/68)

Only 100 of 161 total aligned reads displayed.

GCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCCGGGCGGGCATGGCCCCAGGATCGGCCCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGCCC  >  NC_014034/2198636‑2198833
                                                                                                    |                                                                                                 
gcATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGACGACCGCGATGAGTTCGG.CCCCCGGACGGGCATGGcc                                                                                                  >  2:291106/1‑101 (MQ=255)
  aTGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGccca                                                                                                <  1:2169057/101‑2 (MQ=255)
   tGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccag                                                                                               >  1:2919576/1‑99 (MQ=255)
   tGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccag                                                                                               <  1:4035990/101‑3 (MQ=255)
     aTCGACGAGCAAGCGGCGCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccagga                                                                                             <  2:2691263/101‑5 (MQ=255)
     aTCGAAGAGCAAGAGGCTAACATGATCATGCTTCGGCGAGGGCCGACCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccagaa                                                                                             >  1:2814155/1‑97 (MQ=255)
       cgacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATc                                                                                          <  1:3341876/101‑1 (MQ=255)
         acgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCgg                                                                                        >  2:2948452/1‑101 (MQ=255)
          cgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGc                                                                                       <  1:256603/101‑1 (MQ=255)
          cgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGc                                                                                       >  2:3021832/1‑101 (MQ=255)
           gaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcc                                                                                      >  2:2375603/1‑101 (MQ=255)
           gaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcc                                                                                      <  2:2143714/101‑1 (MQ=255)
            aGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGccc                                                                                     >  1:2441083/1‑101 (MQ=255)
             gCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGGGGGCATGG.CCCAGGATCGGcccg                                                                                    >  1:521049/1‑100 (MQ=255)
              cAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccgg                                                                                   <  1:1600025/101‑3 (MQ=255)
               aaGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccggc                                                                                  <  2:316898/101‑4 (MQ=255)
               aaGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccggc                                                                                  >  2:3511103/1‑98 (MQ=255)
                    gCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGaaa                                                                            <  1:388464/101‑1 (MQ=255)
                     cgcacaTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAg                                                                           >  2:3958185/3‑101 (MQ=255)
                     cTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAg                                                                           >  1:2880783/1‑101 (MQ=255)
                         caTGATCCTGCTTCGGCGAGGGCCGCCCGGGCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcg                                                                       <  2:2479457/101‑1 (MQ=255)
                         caTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcg                                                                       >  1:1376385/1‑101 (MQ=255)
                           tGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGCAAGAtcgtc                                                                     >  1:837190/1‑101 (MQ=255)
                             aTCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGa                                                                   >  1:1762594/1‑101 (MQ=255)
                             aTCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGa                                                                   >  1:3054650/1‑101 (MQ=255)
                              tCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGaa                                                                  <  1:782935/101‑1 (MQ=255)
                               ccTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAAc                                                                 >  1:1313568/1‑101 (MQ=255)
                                cTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACa                                                                <  1:2013097/101‑1 (MQ=255)
                                cTGCTTCGGCGAGGGCCGCCCGGCCCGCGTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACa                                                                <  1:3719278/101‑1 (MQ=255)
                                 tGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAg                                                               <  2:2615854/101‑1 (MQ=255)
                                     tCGGCGAGGGCCGCGCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcacc                                                           >  2:3505992/1‑101 (MQ=255)
                                      cGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcacca                                                          >  1:29652/1‑101 (MQ=255)
                                       ggCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCTGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGGCGAACAGaaccac                                                         >  2:1399253/1‑101 (MQ=255)
                                       ggCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcaccac                                                         >  1:2327761/1‑101 (MQ=255)
                                       ggCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcaccac                                                         <  1:1376982/101‑1 (MQ=255)
                                         cGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACAt                                                       <  2:3243051/101‑1 (MQ=255)
                                         cGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACAt                                                       <  2:957219/101‑1 (MQ=255)
                                          gAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATc                                                      <  2:1253949/101‑1 (MQ=255)
                                          gAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGAACCACATc                                                      >  2:1032393/1‑101 (MQ=255)
                                           aGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCg                                                     <  1:1446626/101‑1 (MQ=255)
                                            gggCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCgg                                                    <  1:904374/101‑1 (MQ=255)
                                              gccgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGcc                                                  >  1:3164514/1‑101 (MQ=255)
                                              gccgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGcc                                                  >  1:3239910/1‑101 (MQ=255)
                                               ccgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCt                                                 >  2:2520772/1‑101 (MQ=255)
                                               ccgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCt                                                 <  2:2813936/101‑1 (MQ=255)
                                                cgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCtt                                                <  1:496549/101‑1 (MQ=255)
                                                cgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCtt                                                <  2:2184359/101‑1 (MQ=255)
                                                 gcccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTg                                               <  2:3679863/101‑1 (MQ=255)
                                                 gcccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTg                                               <  2:3054425/101‑1 (MQ=255)
                                                   ccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGc                                             >  1:4022344/1‑101 (MQ=255)
                                                    cggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCATGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCg                                            <  2:2309139/101‑1 (MQ=255)
                                                     ggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGa                                           <  2:3668947/101‑1 (MQ=255)
                                                      gcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAg                                          <  1:4176533/101‑1 (MQ=255)
                                                        ccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGcc                                        <  1:1182956/101‑1 (MQ=255)
                                                        ccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGcc                                        <  2:3432318/101‑1 (MQ=255)
                                                        ccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGcc                                        <  1:3858492/101‑1 (MQ=255)
                                                         cgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCg                                       <  1:2619683/101‑1 (MQ=255)
                                                           cTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGt                                     >  2:710846/1‑101 (MQ=255)
                                                            ttGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTc                                    <  1:546535/101‑1 (MQ=255)
                                                             tGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCg                                   >  1:3581875/1‑101 (MQ=255)
                                                             tGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCg                                   >  1:2154957/1‑101 (MQ=255)
                                                              gCCGACCGCGATGAGTTCGG.CGCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCgg                                  >  1:2402614/1‑101 (MQ=255)
                                                              gCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCgg                                  >  1:3686391/1‑101 (MQ=255)
                                                                cGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCa                                >  2:3971779/1‑101 (MQ=255)
                                                                 gACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAg                               >  2:2654698/1‑101 (MQ=255)
                                                                 gACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAg                               <  2:3896714/101‑1 (MQ=255)
                                                                  aCCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGc                              >  2:2535311/1‑101 (MQ=255)
                                                                    cgcgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCtt                            <  1:1252940/101‑1 (MQ=255)
                                                                    cgcgATGAGTT.GGGCCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCtt                            <  2:1112827/101‑1 (MQ=255)
                                                                     gcgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCttt                           >  1:3802576/1‑101 (MQ=255)
                                                                      ggATGAGTTCGG.CCCGCGGGGGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTc                          <  1:1042761/100‑1 (MQ=255)
                                                                       gATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCa                         >  2:2758607/1‑101 (MQ=255)
                                                                       gATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCa                         >  1:980539/1‑101 (MQ=255)
                                                                       gATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCa                         >  2:2971430/1‑101 (MQ=255)
                                                                         tGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGc                       <  1:2901157/101‑1 (MQ=255)
                                                                         tGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGc                       >  2:883036/1‑101 (MQ=255)
                                                                            gTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTc                    <  1:1021899/101‑1 (MQ=255)
                                                                            gTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTc                    >  1:4201811/1‑101 (MQ=255)
                                                                             ttCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCCCCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCa                   <  1:1944061/101‑1 (MQ=255)
                                                                             ttCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCa                   >  2:3920099/1‑101 (MQ=255)
                                                                               tcggcccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAg                  >  2:3729346/5‑101 (MQ=255)
                                                                                cggcccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGa                 >  1:2117725/4‑101 (MQ=255)
                                                                                cggcccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGa                 >  2:1825230/4‑101 (MQ=255)
                                                                                cggcccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGa                 <  2:1762475/98‑1 (MQ=255)
                                                                                 ggcccccGGGGGGGCATGG.CCAAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCCCCCCATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGaa                <  2:3239669/99‑1 (MQ=255)
                                                                                  gcccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTTGGCAGCTTTCAGCGTCAGAAg               >  1:1046105/2‑101 (MQ=255)
                                                                                   cccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGc              <  2:2867768/101‑1 (MQ=255)
                                                                                      ccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTcc           <  2:1420549/101‑1 (MQ=255)
                                                                                       cgggcgggCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCg          <  1:4089512/101‑1 (MQ=255)
                                                                                        gggcgggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCgg         >  1:3065375/1‑101 (MQ=255)
                                                                                         ggcgggcATGG.CCCAGGATCGG.CCCGTCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCTAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGa        <  2:1376285/101‑1 (MQ=255)
                                                                                          gcgggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGAt       <  1:760515/101‑1 (MQ=255)
                                                                                           cgggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCCCCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATc      <  2:1743345/101‑1 (MQ=255)
                                                                                            gggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCg     >  1:1060614/1‑101 (MQ=255)
                                                                                            gggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCg     >  2:1616156/1‑101 (MQ=255)
                                                                                            gggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCg     >  2:4069322/1‑101 (MQ=255)
                                                                                            gggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCg     >  2:1394719/1‑101 (MQ=255)
                                                                                            gggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCg     >  2:259183/1‑101 (MQ=255)
                                                                                              gcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGcc   <  1:3022057/101‑1 (MQ=255)
                                                                                               gcatgCCCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGcc   <  1:3753129/97‑1 (MQ=255)
                                                                                                    |                                                                                                 
GCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCCGGGCGGGCATGGCCCCAGGATCGGCCCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGCCC  >  NC_014034/2198636‑2198833

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG