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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,198,7330T→. 1.4% 2.9 ‑3.2 0.0287coding (136/729 nt)cobA2uroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (157/124):  total (159/128)
Fisher's exact test strand distribution p-value = 3.96e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

CAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCCGGGCGGGCATGGCCCCAGGATCGGCCCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGCCC  >  NC_014034/2198632‑2198833
                                                                                                     |                                                                                                    
cAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCAt                                                                                                      >  1:1738200/1‑101 (MQ=255)
cAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCAt                                                                                                      >  1:5535169/1‑101 (MQ=255)
cAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCAt                                                                                                      <  2:2129226/101‑1 (MQ=255)
 aGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATg                                                                                                     >  2:6292162/1‑101 (MQ=255)
    ccATGATGGGCAAGGAAGCGGCTGACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcc                                                                                                  <  2:5072686/100‑1 (MQ=255)
     cATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGccc                                                                                                 <  1:1419282/101‑1 (MQ=255)
       tGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccag                                                                                               <  1:3290863/101‑3 (MQ=255)
       tGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccag                                                                                               <  2:5568909/101‑3 (MQ=255)
        gATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccagg                                                                                              <  1:5687541/101‑4 (MQ=255)
        gATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccagg                                                                                              >  2:235949/1‑98 (MQ=255)
         aTCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGGCGCCCGGCCCGCTTGCCGACCCCGATGAGTTCGG.CCCCCGGGCGGGCGTGGcccagga                                                                                             >  1:6288451/1‑97 (MQ=255)
         aTCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGGcccagga                                                                                             >  2:1813686/1‑97 (MQ=255)
           cgacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATc                                                                                          >  1:1159731/1‑101 (MQ=255)
           cgacgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCAGGCGGGCATGG.CCCAGGATc                                                                                          >  1:4361053/1‑101 (MQ=255)
             acgaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCgg                                                                                        >  2:5282173/1‑101 (MQ=255)
               gaGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcc                                                                                      >  1:5400659/1‑101 (MQ=255)
                aGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGccc                                                                                     <  1:869633/101‑1 (MQ=255)
                aGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGccc                                                                                     >  2:2797990/1‑101 (MQ=255)
                  cAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccgg                                                                                   <  2:6979755/101‑3 (MQ=255)
                  cAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccgg                                                                                   >  2:817765/1‑99 (MQ=255)
                   aaGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccggc                                                                                  <  1:2714669/101‑4 (MQ=255)
                   aaGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccggc                                                                                  <  1:2143251/101‑4 (MQ=255)
                    aGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccggcg                                                                                 >  1:263325/1‑97 (MQ=255)
                    aGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGGcccggcg                                                                                 <  1:1811908/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                                                                            >  2:1206026/1‑101 (MQ=255)
                        gCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGaaa                                                                            <  1:6884436/101‑1 (MQ=255)
                         cTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAg                                                                           <  2:1689107/101‑1 (MQ=255)
                           cacaTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAt                                                                         >  2:3039638/1‑101 (MQ=255)
                            acaTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtc                                                                        <  1:5858619/101‑1 (MQ=255)
                            acaTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtc                                                                        <  2:4344774/101‑1 (MQ=255)
                              aTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcgt                                                                      >  2:4958906/1‑101 (MQ=255)
                              aTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcgt                                                                      >  2:881449/1‑101 (MQ=255)
                              aTGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcgt                                                                      >  1:198271/1‑101 (MQ=255)
                                gATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGAtcgtcg                                                                    >  2:1080865/1‑101 (MQ=255)
                                  tCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGaa                                                                  >  2:200216/1‑101 (MQ=255)
                                   ccTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAAc                                                                 >  1:2801282/1‑101 (MQ=255)
                                     tGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAg                                                               <  2:4002543/101‑1 (MQ=255)
                                     tGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAg                                                               <  1:3539754/101‑1 (MQ=255)
                                      gCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGc                                                              >  1:78119/1‑101 (MQ=255)
                                         tCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcacc                                                           <  2:56739/101‑1 (MQ=255)
                                         tCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGcacc                                                           <  1:1498892/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                                                      >  2:2484934/1‑101 (MQ=255)
                                              gAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATc                                                      <  1:5116433/101‑1 (MQ=255)
                                               aGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCg                                                     <  2:4467247/101‑1 (MQ=255)
                                                gggCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATctg                                                    >  2:360094/1‑99 (MQ=255)
                                                gggCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCgg                                                    <  1:664224/101‑1 (MQ=255)
                                                   ccgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.TCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCt                                                 <  1:3376646/101‑1 (MQ=255)
                                                   ccgccCGGCCCGCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCt                                                 >  1:4755098/1‑101 (MQ=255)
                                                     gcccggcccgTTCGCCGATCGCTATGTGTCCGG.CCCCCGTGCGGGCTTTG.ACCAGTATCGG.ACCGTCGGAATGATCGTCGAACAGCACCACATCGGCCTTg                                               <  2:1271587/101‑1 (MQ=255)
                                                       ccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGc                                             <  1:3810330/101‑1 (MQ=255)
                                                       ccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGc                                             <  2:1585807/101‑1 (MQ=255)
                                                       ccggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGc                                             >  2:2882949/1‑101 (MQ=255)
                                                         ggccctCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGa                                           <  1:1001551/101‑1 (MQ=255)
                                                         ggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGa                                           >  1:401136/1‑101 (MQ=255)
                                                         ggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGa                                           >  1:2332542/1‑101 (MQ=255)
                                                         ggcccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGa                                           >  1:4650629/1‑101 (MQ=255)
                                                           cccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGc                                         >  1:4417144/1‑101 (MQ=255)
                                                           cccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGc                                         <  1:3197609/101‑1 (MQ=255)
                                                            ccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGcc                                        >  2:131669/1‑101 (MQ=255)
                                                            ccgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGcc                                        <  2:7052837/101‑1 (MQ=255)
                                                             cgCTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCg                                       <  1:5187400/101‑1 (MQ=255)
                                                             cccTTGCCGCCCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCg                                       <  2:348601/99‑1 (MQ=255)
                                                               cTTGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGt                                     >  2:4845671/1‑101 (MQ=255)
                                                                ttGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTc                                    <  1:4959158/101‑1 (MQ=255)
                                                                 tGCCGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCg                                   >  2:2972725/1‑101 (MQ=255)
                                                                   ccGACCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGc                                 <  1:425814/101‑1 (MQ=255)
                                                                      aCCGCGATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGc                              >  2:6127567/1‑101 (MQ=255)
                                                                         gcgATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCttt                           >  1:3964621/1‑101 (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:997522/101‑1 (MQ=255)
                                                                           gATGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCa                         <  1:650824/101‑1 (MQ=255)
                                                                             tGAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGc                       >  1:6197993/1‑101 (MQ=255)
                                                                              gAGTTCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCg                      >  2:4226613/1‑101 (MQ=255)
                                                                                 ttCGG.CCCCCGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCa                   <  2:6039549/101‑1 (MQ=255)
                                                                                    cggcccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGa                 >  1:2556886/4‑101 (MQ=255)
                                                                                      gcccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAg               >  2:4076173/2‑101 (MQ=255)
                                                                                       cccccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGc              <  2:1652961/101‑1 (MQ=255)
                                                                                         cccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCCCCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTc            <  1:3100160/101‑1 (MQ=255)
                                                                                         cccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTc            >  2:6621763/1‑101 (MQ=255)
                                                                                          ccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTcc           <  2:5120161/101‑1 (MQ=255)
                                                                                          ccGGGCGGGCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTcc           <  2:5846994/101‑1 (MQ=255)
                                                                                          ccGGGCGGGCA.GGGCCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACCGCACCCCATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTcc           <  2:3859722/101‑1 (MQ=255)
                                                                                           cgggcgggCATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCg          >  1:7092886/1‑101 (MQ=255)
                                                                                            gggcgggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCgg         <  1:4503929/101‑1 (MQ=255)
                                                                                            gggcgggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCgg         >  2:5362944/1‑101 (MQ=255)
                                                                                             ggcgggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGa        >  2:6489173/1‑101 (MQ=255)
                                                                                               cgggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATc      >  2:3792138/1‑101 (MQ=255)
                                                                                                gggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATcc     >  2:4203446/1‑100 (MQ=255)
                                                                                                gggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCg     >  1:656358/1‑101 (MQ=255)
                                                                                                gggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCg     >  2:6706566/1‑101 (MQ=255)
                                                                                                 ggcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGc    >  2:6886467/1‑101 (MQ=255)
                                                                                                  gcATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGcc   <  2:2687333/101‑1 (MQ=255)
                                                                                                   cATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGccc  <  1:1687015/101‑1 (MQ=255)
                                                                                                   cATGG.CCCAGGATCGG.CCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGccc  >  1:2180609/1‑101 (MQ=255)
                                                                                                     |                                                                                                    
CAGCGCATGATCGACGAGCAAGCGGCTCACATGATCCTGCTTCGGCGAGGGCCGCCCGGCCCGCTTGCCGACCGCGATGAGTTCGGCCCCCCGGGCGGGCATGGCCCCAGGATCGGCCCCGGCGGAAAGATCGTCGAACAGCACCACATCGGCCTTGGCGAGCCGGTCGGCAGCTTTCAGCGTCAGAAGCTCCGGATCGCCC  >  NC_014034/2198632‑2198833

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