Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,198,736 +G 1.5% 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,7360.→G 1.5% 4.7 ‑3.1 0.0259coding (133/729 nt)cobA2uroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (140/115):  total (140/119)
Fisher's exact test strand distribution p-value = 4.33e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

CGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCCGGGCGGGCATGGCCCCAGGATCGGCCCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGCCC  >  NC_014034/2198635‑2198833
                                                                                                                                                                                                       
  cATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGccc                                                                                                 >  1:4002718/1‑101 (MQ=255)
    tGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccag                                                                                               <  2:5568909/101‑3 (MQ=255)
     gATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccagg                                                                                              >  2:235949/1‑98 (MQ=255)
     gATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccagg                                                                                              <  2:2115073/101‑4 (MQ=255)
     gATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccagg                                                                                              >  2:1687145/1‑98 (MQ=255)
     gATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccagg                                                                                              >  1:3571354/1‑98 (MQ=255)
     gATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccagg                                                                                              <  1:5687541/101‑4 (MQ=255)
      aTCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGGCGCCCGGCCCGCTTGCCGACCCCGATGAGTTCGG.CCCCCGGGCGGGCGTGGcccagga                                                                                             >  1:6288451/1‑97 (MQ=255)
      aTCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGGCCGGCCGGCTTGCCGACACGGATGAGTTGGGCGCCCC.GGCGGGCATGGCCCCGgg                                                                                             >  2:1579271/1‑101 (MQ=255)
      aTCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCC.GGCGGGCATGGcccagga                                                                                             >  2:4760958/1‑97 (MQ=255)
       tCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGAt                                                                                           <  2:4226291/101‑1 (MQ=255)
        cgacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCAGGCGGGCATGG.CCCAGGATc                                                                                          >  1:4361053/1‑101 (MQ=255)
         gacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCg                                                                                         <  2:1306398/101‑1 (MQ=255)
          acgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCgg                                                                                        >  1:5160051/1‑101 (MQ=255)
            gaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcc                                                                                      <  1:961415/101‑1 (MQ=255)
             aGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGccc                                                                                     >  2:2797990/1‑101 (MQ=255)
              gCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccg                                                                                    <  1:3227615/101‑2 (MQ=255)
               cAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccgg                                                                                   <  1:79599/101‑3 (MQ=255)
                 aGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccggcg                                                                                 >  1:263325/1‑97 (MQ=255)
                 aGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccggcg                                                                                 >  1:2036736/1‑97 (MQ=255)
                 aGCGGCTCACAGGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccggcg                                                                                 <  1:1411669/101‑5 (MQ=255)
                    ggCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGaa                                                                             <  2:5429550/101‑1 (MQ=255)
                    ggCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGaa                                                                             >  2:1121626/1‑101 (MQ=255)
                     gCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGaaa                                                                            <  1:6884436/101‑1 (MQ=255)
                     gCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGaaa                                                                            >  2:1206026/1‑101 (MQ=255)
                       tCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGa                                                                          >  2:838358/1‑101 (MQ=255)
                        cacaTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAt                                                                         >  2:5116173/1‑101 (MQ=255)
                        cacaTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAt                                                                         >  2:3039638/1‑101 (MQ=255)
                         acaTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtc                                                                        <  1:5858619/101‑1 (MQ=255)
                           aTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcgt                                                                      >  2:1632656/1‑101 (MQ=255)
                            tGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcgtc                                                                     >  2:3822281/1‑101 (MQ=255)
                             gATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcgtcg                                                                    >  2:4905304/1‑101 (MQ=255)
                             gATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcgtcg                                                                    >  2:1080865/1‑101 (MQ=255)
                                  ggCTCCGGCGAGGGCCGCCCGGCGCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAg                                                               <  2:4986540/100‑1 (MQ=255)
                                   gCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGc                                                              <  2:5952452/101‑1 (MQ=255)
                                     ttCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcac                                                            <  1:1813750/101‑1 (MQ=255)
                                      tCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcacc                                                           <  1:1498892/101‑1 (MQ=255)
                                      tCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcacc                                                           <  2:56739/101‑1 (MQ=255)
                                      tCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcacc                                                           <  2:3698372/101‑1 (MQ=255)
                                        ggCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcaccac                                                         <  1:1579332/101‑1 (MQ=255)
                                         tcctgggCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACcaca                                                        >  2:5162813/5‑101 (MQ=255)
                                         gCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACcaca                                                        >  1:4100270/1‑101 (MQ=255)
                                           gAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATc                                                      <  1:317994/101‑1 (MQ=255)
                                           gAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATc                                                      <  1:72964/101‑1 (MQ=255)
                                           gAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATc                                                      >  2:2484934/1‑101 (MQ=255)
                                           gAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATc                                                      >  2:4311338/1‑101 (MQ=255)
                                            aGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.GCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCCCATCg                                                     >  2:4178784/1‑101 (MQ=255)
                                            aGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCg                                                     >  1:4667948/1‑101 (MQ=255)
                                             gggCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCgg                                                    <  1:664224/101‑1 (MQ=255)
                                              ggCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGc                                                   >  2:6715480/1‑101 (MQ=255)
                                              ggCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGAT.GGGCCCGGCGGAAAGATCGTCGAACAGCACCACATCGGc                                                   >  1:5093427/1‑101 (MQ=255)
                                                 cgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCtt                                                <  1:5512431/101‑1 (MQ=255)
                                                  gcccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTg                                               >  2:6340054/1‑101 (MQ=255)
                                                   cccggcccgCTTGCGGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTgg                                              <  2:1263722/101‑1 (MQ=255)
                                                   cccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTgg                                              >  2:5317827/1‑101 (MQ=255)
                                                   cccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTgg                                              <  2:3558851/101‑1 (MQ=255)
                                                    ccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGc                                             <  1:4761198/101‑1 (MQ=255)
                                                    ccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGc                                             <  1:5210790/101‑1 (MQ=255)
                                                    ccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGc                                             >  2:2882949/1‑101 (MQ=255)
                                                     cggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGcc                                            >  2:3864439/1‑100 (MQ=255)
                                                      ggccctCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGa                                           <  1:1001551/101‑1 (MQ=255)
                                                      ggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGa                                           >  1:401136/1‑101 (MQ=255)
                                                        cccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGc                                         >  1:4417144/1‑101 (MQ=255)
                                                         ccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGcc                                        <  2:7091118/101‑1 (MQ=255)
                                                          cgCTTGCCGGCCGCGATGAGTGCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCg                                       <  2:1738136/101‑1 (MQ=255)
                                                           gCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCgg                                      >  2:886724/1‑101 (MQ=255)
                                                              tGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCg                                   >  2:2972725/1‑101 (MQ=255)
                                                              tGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCg                                   >  1:3428518/1‑101 (MQ=255)
                                                                ccGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGc                                 <  1:425814/101‑1 (MQ=255)
                                                                      gcgACGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCttt                           <  2:4523471/101‑1 (MQ=255)
                                                                       ctaTGAGTTCGG.CCCCCGGGCGGGCATGG.CCCGGGATC.GCCCCGGCGGAAATATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTc                          <  2:3279221/99‑1 (MQ=255)
                                                                       cgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTc                          <  1:292199/101‑1 (MQ=255)
                                                                       cgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTc                          <  2:6668120/101‑1 (MQ=255)
                                                                        gATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCa                         >  1:2353368/1‑101 (MQ=255)
                                                                        gATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCa                         >  2:4376268/1‑101 (MQ=255)
                                                                        gATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCa                         <  1:650824/101‑1 (MQ=255)
                                                                         aTGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAg                        <  1:2798092/101‑1 (MQ=255)
                                                                          tGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGc                       >  1:6197993/1‑101 (MQ=255)
                                                                             gTTCGG.CCCCGGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTc                    <  1:2319574/101‑1 (MQ=255)
                                                                             gTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTc                    <  2:6120862/101‑1 (MQ=255)
                                                                              ttCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCa                   <  1:2148209/101‑1 (MQ=255)
                                                                              gtgcggcccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTc                    <  1:3325798/95‑1 (MQ=255)
                                                                                 cggcccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGa                 >  1:2556886/4‑101 (MQ=255)
                                                                                 cggcccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGa                 >  2:6911054/4‑101 (MQ=255)
                                                                                 cggcccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGa                 <  1:1752246/98‑1 (MQ=255)
                                                                                  ggcccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGaa                >  2:6988443/3‑101 (MQ=255)
                                                                                  ggcccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGAGTCAGac                >  1:5335763/3‑100 (MQ=255)
                                                                                    cccccGGGCGGGCA.GGGCCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGc              <  1:1115742/101‑1 (MQ=255)
                                                                                       ccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTcc           <  2:5846994/101‑1 (MQ=255)
                                                                                       ccGGGCGGGCA.GGGCCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTcc           <  2:4065594/101‑1 (MQ=255)
                                                                                        cgggcgggCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCg          >  1:7092886/1‑101 (MQ=255)
                                                                                        cgggcgggCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCg          >  2:3357017/1‑101 (MQ=255)
                                                                                         gggcgggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCgg         <  1:4503929/101‑1 (MQ=255)
                                                                                         gggcgggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCgg         >  2:5362944/1‑101 (MQ=255)
                                                                                           gcgggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGCCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGAt       <  2:6534737/101‑1 (MQ=255)
                                                                                             gggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCg     >  1:656358/1‑101 (MQ=255)
                                                                                             gggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCg     <  1:3810731/101‑1 (MQ=255)
                                                                                             gggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCg     <  2:181442/101‑1 (MQ=255)
                                                                                               gcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGcc   <  2:2687333/101‑1 (MQ=255)
                                                                                                cATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGccc  <  1:1687015/101‑1 (MQ=255)
                                                                                                                                                                                                       
CGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCCGGGCGGGCATGGCCCCAGGATCGGCCCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGCCC  >  NC_014034/2198635‑2198833

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 28 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG