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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,198,7180C→.97.3% 8.5 ‑0.4 ‑0.1185coding (151/729 nt)cobA2uroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  new base (94/86):  ref base (2/3):  total (96/89)
Fisher's exact test strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.41e-01

Only 100 of 218 total aligned reads displayed.

ACCCGTTGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCCGGGCGGGCATGGCCCCAGGATCGGCCCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAA  >  NC_014034/2198618‑2198819
                                                                                                    |                                                                                                     
accCGTTGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGc                                                                                                       >  2:1652153/1‑101 (MQ=255)
 cccGTTGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGcc                                                                                                      <  2:1359218/101‑1 (MQ=255)
  ccGTTGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGccc                                                                                                     <  2:3509519/101‑1 (MQ=255)
    gTTGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGccccc                                                                                                   >  1:2969561/1‑101 (MQ=255)
    gTTGCCCGGTCAGCGCATGAGCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGccccc                                                                                                   <  2:2106472/101‑1 (MQ=255)
     ttGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGcccccg                                                                                                  <  2:1559500/101‑2 (MQ=255)
     ttGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGcccccg                                                                                                  >  2:4104368/1‑100 (MQ=255)
     ttGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGcccccg                                                                                                  <  2:427313/101‑2 (MQ=255)
      tGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGcccccgg                                                                                                 >  1:2789812/1‑99 (MQ=255)
       gCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCggg                                                                                                <  1:1265231/101‑1 (MQ=255)
        cccGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCgggc                                                                                               <  1:2346281/101‑2 (MQ=255)
           ggTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCcgggcggg                                                                                           >  2:4210416/1‑101 (MQ=255)
            gTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCgggccggc                                                                                           >  2:1443590/1‑100 (MQ=255)
              cAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCAt                                                                                        >  2:4551510/1‑101 (MQ=255)
              cAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCAt                                                                                        <  2:3497108/101‑1 (MQ=255)
                   cATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGccc                                                                                   >  1:4668238/1‑101 (MQ=255)
                   cATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGccc                                                                                   >  1:670066/1‑101 (MQ=255)
                     tGATCAACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccag                                                                                 <  2:2115826/101‑3 (MQ=255)
                        tCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGAt                                                                             <  2:585060/101‑1 (MQ=255)
                        tCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGAt                                                                             <  1:1089396/101‑1 (MQ=255)
                         cgacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATc                                                                            <  2:317849/101‑1 (MQ=255)
                         cgacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATc                                                                            >  1:75317/1‑101 (MQ=255)
                          gacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCg                                                                           >  1:605598/1‑101 (MQ=255)
                          gacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCg                                                                           <  2:4357333/101‑1 (MQ=255)
                          gacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CACCCGGGCGGGCATGG.CCCAGGATCg                                                                           >  2:1562917/1‑101 (MQ=255)
                           acgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCgg                                                                          >  1:4099991/1‑101 (MQ=255)
                             gaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcc                                                                        >  2:1249873/1‑101 (MQ=255)
                             gaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcc                                                                        <  2:2821927/101‑1 (MQ=255)
                                  aGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccggcg                                                                   >  2:2727088/1‑97 (MQ=255)
                                   gCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCcggcgg                                                                 <  2:2049195/101‑1 (MQ=255)
                                    cGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGa                                                                <  2:2323147/101‑1 (MQ=255)
                                    cGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGa                                                                <  2:839615/101‑1 (MQ=255)
                                     ggCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGaa                                                               >  1:2397059/1‑101 (MQ=255)
                                      gCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGaaa                                                              >  2:1928660/1‑101 (MQ=255)
                                       cTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCTGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAg                                                             >  1:1602864/1‑101 (MQ=255)
                                        tCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGa                                                            <  2:3316442/101‑1 (MQ=255)
                                          acaTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCACGATCGG.CCCGGCGGAAAGAtc                                                          <  1:2639406/101‑1 (MQ=255)
                                           cgtGATCCTGCTTCGGGGAGGGCCGCCTGGCCCGCTTGCCGACCGCGATGGGTT.GGGCCGCCGGGGGGGCAT.GCCCCAGGATCGG.CCCGGCGGAAAGAtcg                                                         <  2:2265635/99‑1 (MQ=255)
                                           caTGATCCTTCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcg                                                         <  2:670004/101‑1 (MQ=255)
                                            aTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcgt                                                        >  2:1608016/1‑101 (MQ=255)
                                            aTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcgt                                                        >  2:147151/1‑101 (MQ=255)
                                            aTGATCCTGCTTCGGCGAGGCCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGCAAGAtcgc                                                        >  2:2371706/1‑100 (MQ=255)
                                              gATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcgtcg                                                      <  2:47122/101‑1 (MQ=255)
                                              gATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcgtcg                                                      <  1:1080467/101‑1 (MQ=255)
                                               aTCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGa                                                     >  2:2643294/1‑101 (MQ=255)
                                                tCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGaa                                                    <  1:4551809/101‑1 (MQ=255)
                                                tCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGGCGaa                                                    >  2:3547816/1‑101 (MQ=255)
                                                 ncTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCNNCNnan                                                   >  1:1835330/2‑100 (MQ=255)
                                                 ccTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAAc                                                   <  1:952949/101‑1 (MQ=255)
                                                  cTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACa                                                  <  1:1439632/101‑1 (MQ=255)
                                                    gCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGc                                                >  1:1769833/1‑101 (MQ=255)
                                                      ttCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcac                                              >  1:4398760/1‑101 (MQ=255)
                                                        cGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcaccc                                            >  1:3311414/1‑100 (MQ=255)
                                                        cGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcacca                                            >  2:1966788/1‑101 (MQ=255)
                                                        cGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcacca                                            >  2:245516/1‑101 (MQ=255)
                                                        cGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcacca                                            <  2:2588154/101‑1 (MQ=255)
                                                           cGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACAt                                         <  2:1769714/101‑1 (MQ=255)
                                                           cGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACAt                                         >  1:3520037/1‑101 (MQ=255)
                                                            gAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATc                                        <  1:3912517/101‑1 (MQ=255)
                                                            gAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATc                                        >  2:924114/1‑101 (MQ=255)
                                                            gAGGG.CGCCCGGCCCGCTTGCCGACCGCGGTGAGGTTGG.GCCCCGGGCGGGCATTG.GCCAGGGTCGG.CCCGGCGGGAAGGTTGGCCGACAGgaccaccacg                                       >  2:3720966/1‑97 (MQ=255)
                                                             aGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCg                                       <  1:4118082/101‑1 (MQ=255)
                                                              gggCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCgg                                      <  1:1554642/101‑1 (MQ=255)
                                                              gggCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCgg                                      <  1:4104654/101‑1 (MQ=255)
                                                                gccgccCGGCCCGCTTGCCGACCGCGCTGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGcc                                    <  1:2803903/101‑1 (MQ=255)
                                                                gccgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGcc                                    >  2:184964/1‑101 (MQ=255)
                                                                  cgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCtt                                  >  2:909317/1‑101 (MQ=255)
                                                                    cccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTgg                                <  2:3783627/101‑1 (MQ=255)
                                                                    attggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATAGGCCCTgg                                >  2:1629900/4‑101 (MQ=255)
                                                                      cggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCg                              >  2:2035509/1‑101 (MQ=255)
                                                                      cggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCg                              <  2:3363729/101‑1 (MQ=255)
                                                                       ggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGa                             >  1:1188280/1‑101 (MQ=255)
                                                                       ggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGa                             >  1:1360987/1‑101 (MQ=255)
                                                                        gcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAg                            >  1:2498806/1‑101 (MQ=255)
                                                                        gcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAg                            <  1:3397300/101‑1 (MQ=255)
                                                                          ccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGcc                          <  1:2839623/101‑1 (MQ=255)
                                                                           cgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCg                         <  1:4515810/101‑1 (MQ=255)
                                                                            gCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCgg                        >  2:2483906/1‑101 (MQ=255)
                                                                            gCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCgg                        >  1:887856/1‑101 (MQ=255)
                                                                             cTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGt                       >  1:402986/1‑101 (MQ=255)
                                                                              ttGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTc                      >  2:3439158/1‑101 (MQ=255)
                                                                               tGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCg                     >  2:797623/1‑101 (MQ=255)
                                                                                  cGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACACCCCCACATCGGCCTTGGCGAGCCGGTCGGCa                  <  1:4157881/101‑1 (MQ=255)
                                                                                   gACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCTTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAg                 <  2:4205302/101‑1 (MQ=255)
                                                                                   gACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGGGAGCCGACCGGCAg                 >  2:1915887/1‑101 (MQ=255)
                                                                                      cgcgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCtt              <  1:2986768/101‑1 (MQ=255)
                                                                                       gcgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCttt             <  2:2418142/101‑1 (MQ=255)
                                                                                        cgATGAGTTCGG.CCCCCGGTGGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTc            <  2:516185/101‑1 (MQ=255)
                                                                                        cgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGGAGCTTTc            >  1:824172/1‑101 (MQ=255)
                                                                                        cgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTc            >  2:111674/1‑101 (MQ=255)
                                                                                        cgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTc            >  1:3795154/1‑101 (MQ=255)
                                                                                         gATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCa           <  1:1563022/101‑1 (MQ=255)
                                                                                         gATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCa           >  2:59085/1‑101 (MQ=255)
                                                                                          aTGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAg          >  1:3648458/1‑101 (MQ=255)
                                                                                          aTGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAg          >  1:4612917/1‑101 (MQ=255)
                                                                                             aGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGt       <  1:1966923/101‑1 (MQ=255)
                                                                                             aGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGt       <  2:3024458/101‑1 (MQ=255)
                                                                                               ttCGG.CCCCCGGGGGGGCA.GGGCCCGGGGTCGG.CCCGGCGGAAAGATCGTCGAACAGCCCCCCATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCa     <  1:3143118/101‑1 (MQ=255)
                                                                                               ttCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCa     <  1:2540014/101‑1 (MQ=255)
                                                                                                   gCCCCCCGGGCGGGCCTCGCCCC.GGAAC.GCCCCCGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGaa  <  2:824096/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
ACCCGTTGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCCGGGCGGGCATGGCCCCAGGATCGGCCCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAA  >  NC_014034/2198618‑2198819

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG