Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,198,718 Δ1 bp 96.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→.96.3% 20.4 0.0 ‑0.2242coding (151/729 nt)cobA2uroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  new base (110/123):  ref base (4/5):  total (114/128)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.91e-01

Only 100 of 262 total aligned reads displayed.

ACCCGTTGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCCGGGCGGGCATGGCCCCAGGATCGGCCCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCA  >  NC_014034/2198618‑2198816
                                                                                                    |                                                                                                  
  ccGTTGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCAGACCGCGATGAGTTCGGccc                                                                                                  <  1:5598613/101‑1 (MQ=255)
   cGTTGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGcccc                                                                                                 <  1:1064162/101‑1 (MQ=255)
   cGTTGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGcccc                                                                                                 >  2:5144955/1‑101 (MQ=255)
   cGTTGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCACCCGGCCCGCTTGCCGACCGCGATGAGTTCGGcccc                                                                                                 <  2:301531/101‑1 (MQ=255)
   cGTTGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTACTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGCGTTCGGcccc                                                                                                 >  1:4899919/1‑101 (MQ=255)
       gCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCggg                                                                                             <  2:2696910/101‑1 (MQ=255)
        cccGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCgggc                                                                                            <  2:5490212/101‑2 (MQ=255)
        cccGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCgggc                                                                                            <  1:378833/101‑2 (MQ=255)
        cccGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCgggc                                                                                            >  2:4671922/1‑100 (MQ=255)
         ccGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCcgggcg                                                                                          >  1:5286605/1‑101 (MQ=255)
          cGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCcgggcgg                                                                                         >  2:165378/1‑101 (MQ=255)
          cGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCcgggcgg                                                                                         >  2:2761214/1‑101 (MQ=255)
             tCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCa                                                                                      >  1:3813717/1‑101 (MQ=255)
             tCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCa                                                                                      >  2:525385/1‑101 (MQ=255)
              cAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCAt                                                                                     <  1:4598988/101‑1 (MQ=255)
               aGCGCATGATCGGCGAGAAAGCGGCTCCTATGATGCTGCTTCGGCAAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGGTCGG.CCCCGGGGCGGGCATg                                                                                    >  2:1475372/1‑101 (MQ=255)
               aGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATg                                                                                    <  1:3206542/101‑1 (MQ=255)
                  gcATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATggac                                                                                 >  2:3898340/1‑99 (MQ=255)
                  gcATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcc                                                                                 >  2:4336044/1‑101 (MQ=255)
                   cATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGccc                                                                                <  2:1231491/101‑1 (MQ=255)
                    aTGATCGACTAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGccca                                                                               >  1:3562286/1‑100 (MQ=255)
                    aTGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGccca                                                                               <  2:680920/101‑2 (MQ=255)
                      gATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccagg                                                                             <  1:2047306/101‑4 (MQ=255)
                        tCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGAt                                                                          <  1:4513299/101‑1 (MQ=255)
                         cgacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATc                                                                         <  2:437273/101‑1 (MQ=255)
                         cgacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGCTCGG.CCCCCGGGCGGGCATGG.CCCAGGATc                                                                         <  1:5145241/101‑1 (MQ=255)
                          gacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCC.GGCGGGCATGG.CCCAGGATCg                                                                        <  1:3286252/101‑1 (MQ=255)
                          gacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCg                                                                        <  1:4137913/101‑1 (MQ=255)
                          gacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCCCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCg                                                                        >  2:456225/1‑101 (MQ=255)
                            cgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGc                                                                      >  1:2436755/1‑101 (MQ=255)
                             gaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.GCCAGGATCGGcc                                                                     >  2:1882259/1‑101 (MQ=255)
                              aGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGGCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGAGCGGccc                                                                    >  2:2639112/1‑101 (MQ=255)
                              aGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGccc                                                                    <  2:1778217/101‑1 (MQ=255)
                               gCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccg                                                                   >  2:2165202/1‑100 (MQ=255)
                                 aaGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccggc                                                                 >  2:1373839/1‑98 (MQ=255)
                                 aaGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccggc                                                                 <  1:1460813/101‑4 (MQ=255)
                                 aaGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccggc                                                                 >  2:5035035/1‑98 (MQ=255)
                                    cGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGa                                                             >  1:5723726/1‑101 (MQ=255)
                                      gCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGaaa                                                           >  1:2618449/1‑101 (MQ=255)
                                       cTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAg                                                          <  1:4059157/101‑1 (MQ=255)
                                       cTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAg                                                          >  2:3554835/1‑101 (MQ=255)
                                       cTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAg                                                          <  1:1647870/101‑1 (MQ=255)
                                        tCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGCCCGG.CCCCCGGGCGGCCAT.GCACCAGGTTCGG.CCAGGCGGAGAGa                                                         >  1:839191/1‑101 (MQ=255)
                                         cacaTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAt                                                        >  2:616339/1‑101 (MQ=255)
                                          acaTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtc                                                       <  2:1953235/101‑1 (MQ=255)
                                          acaTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtc                                                       <  2:5546231/101‑1 (MQ=255)
                                           caTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcg                                                      <  1:1725697/101‑1 (MQ=255)
                                           caTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcg                                                      >  1:4817843/1‑101 (MQ=255)
                                               aTCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGa                                                  <  1:2767100/101‑1 (MQ=255)
                                                 ccTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAAc                                                <  1:4672198/101‑1 (MQ=255)
                                                    gCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGc                                             >  2:5767995/1‑101 (MQ=255)
                                                     cTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGca                                            >  1:2406385/1‑101 (MQ=255)
                                                      cgcGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcac                                           <  1:3871502/99‑1 (MQ=255)
                                                        cGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcacca                                         >  1:3347671/1‑101 (MQ=255)
                                                        cGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcacca                                         >  1:1082387/1‑101 (MQ=255)
                                                         ggCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcaccac                                        >  2:2442249/1‑101 (MQ=255)
                                                         ggCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcaccac                                        >  2:3800166/1‑101 (MQ=255)
                                                         ggCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcaccac                                        <  2:5399356/101‑1 (MQ=255)
                                                         ggCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcaccac                                        <  2:4817562/101‑1 (MQ=255)
                                                          gCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACcaca                                       <  1:2737595/101‑1 (MQ=255)
                                                           cGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACAt                                      <  2:4964803/101‑1 (MQ=255)
                                                            gAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATc                                     >  2:5667665/1‑101 (MQ=255)
                                                            gAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATc                                     >  2:5374080/1‑101 (MQ=255)
                                                            gAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATc                                     <  1:2664709/101‑1 (MQ=255)
                                                            gAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCAACACATc                                     >  2:2863294/1‑101 (MQ=255)
                                                              gggCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGCACAGCCCCACATCgg                                   >  2:529630/1‑101 (MQ=255)
                                                              gggCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCgg                                   <  1:311772/101‑1 (MQ=255)
                                                               ggCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGc                                  <  1:3497917/101‑1 (MQ=255)
                                                                gccgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCTGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGcc                                 <  2:69116/101‑1 (MQ=255)
                                                                gccgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGcc                                 <  2:2278852/101‑1 (MQ=255)
                                                                gccgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGcc                                 <  1:3055558/101‑1 (MQ=255)
                                                                 ccgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCt                                <  1:658093/101‑1 (MQ=255)
                                                                  cgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCtt                               <  1:4626293/101‑1 (MQ=255)
                                                                   gcccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTg                              <  2:3408002/101‑1 (MQ=255)
                                                                   gcccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTg                              <  2:5350836/101‑1 (MQ=255)
                                                                   gcccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTg                              >  2:4024391/1‑101 (MQ=255)
                                                                    cccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTgg                             <  2:3562069/101‑1 (MQ=255)
                                                                    cccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTgg                             <  1:1522456/101‑1 (MQ=255)
                                                                     ccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGc                            >  1:3701122/1‑101 (MQ=255)
                                                                      cggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCg                           <  1:997706/101‑1 (MQ=255)
                                                                         cccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGc                        <  1:230428/101‑1 (MQ=255)
                                                                          ccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGcc                       <  1:4872427/101‑1 (MQ=255)
                                                                              ttGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTc                   >  2:2539550/1‑101 (MQ=255)
                                                                              ttGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTc                   >  1:1589157/1‑101 (MQ=255)
                                                                               tGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCg                  >  1:4172015/1‑101 (MQ=255)
                                                                                gCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCCCCACATCGGCCTTGGCGAGCCGGTCgg                 <  2:884199/101‑1 (MQ=255)
                                                                                    aCCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGc             >  2:3427375/1‑101 (MQ=255)
                                                                                    aCCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGc             >  2:1914435/1‑101 (MQ=255)
                                                                                     ccGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCt            <  1:4336310/101‑1 (MQ=255)
                                                                                      cgcgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCtt           <  1:5213680/101‑1 (MQ=255)
                                                                                      cgcgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCtt           >  2:3272039/1‑101 (MQ=255)
                                                                                      cgcgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCtt           >  2:3076208/1‑101 (MQ=255)
                                                                                       gcgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCttt          >  1:1256433/1‑101 (MQ=255)
                                                                                        cgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTc         >  1:3992624/1‑101 (MQ=255)
                                                                                        cgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTc         >  2:1563333/1‑101 (MQ=255)
                                                                                          aTGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAg       <  2:2070820/101‑1 (MQ=255)
                                                                                            gAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCg     >  1:2150561/1‑101 (MQ=255)
                                                                                            gAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCg     <  2:4255577/101‑1 (MQ=255)
                                                                                             aGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGt    <  1:1357169/101‑1 (MQ=255)
                                                                                                ttcgCCCCCCGGCCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGACCAGCACCACTTCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCa  <  1:5498903/98‑1 (MQ=255)
                                                                                                    |                                                                                                  
ACCCGTTGCCCGGTCAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCCGGGCGGGCATGGCCCCAGGATCGGCCCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCA  >  NC_014034/2198618‑2198816

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