Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,198,487 A→C 100% S128A (TCC→GCC)  cobA2 ← uroporphyrinogen‑III C‑methyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,198,4870A→C100.0% 537.8 149S128A (TCC→GCC) cobA2uroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  new base (68/81):  ref base (0/0):  total (68/81)

Only 100 of 177 total aligned reads displayed.

CGCCCAGTTCAGCCCCTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGACGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGTTCTTCCTCAAGCCGGCCG  >  NC_014034/2198389‑2198584
                                                                                                  |                                                                                                 
cgccCAGTTCAGCCCCTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCcg                                                                                                 <  1:1879690/101‑1 (MQ=255)
 gccCAGTTCAGCCCCTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCcgc                                                                                                >  1:4014605/1‑101 (MQ=255)
   ccAGTTCAGCCCCTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCgcgg                                                                                              <  1:1440843/101‑1 (MQ=255)
     ggTTCAGCCCCTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCgcggcg                                                                                            <  1:1209159/100‑1 (MQ=255)
      gTTCAGCCCCTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCgcggcgg                                                                                           >  1:3757721/1‑101 (MQ=255)
       ttCAGCCCCTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCgcggcggc                                                                                          >  2:1021774/1‑101 (MQ=255)
       ttCAGCCCCTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCgcggcggc                                                                                          <  2:1304725/101‑1 (MQ=255)
        tCAGCCCCTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCgcggcgggg                                                                                         >  2:1008172/1‑99 (MQ=255)
          aGCCCCTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGGAGCCCGGGCccggcggcggc                                                                                       >  2:2621932/1‑101 (MQ=255)
          aGCCCCTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCgcggcggcggc                                                                                       >  1:3861058/1‑101 (MQ=255)
             cccTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGGCAGCGGCAGCCCGGCCgcggcggcggcggc                                                                                    >  1:293846/1‑101 (MQ=255)
               cTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCCCCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGgcgc                                                                                  <  1:4618790/101‑1 (MQ=255)
                ccGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCa                                                                                 <  2:3210191/100‑1 (MQ=255)
                 cGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAg                                                                                >  1:1867130/1‑101 (MQ=255)
                    gCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGCCGGCGGCGCAGGcc                                                                             <  2:2934520/101‑1 (MQ=255)
                    gCAGACCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGcc                                                                             <  1:2630836/101‑1 (MQ=255)
                      aGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCGGGCCGa                                                                           >  1:4277027/1‑101 (MQ=255)
                      aGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCgg                                                                           >  2:3665409/1‑100 (MQ=255)
                      aGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGa                                                                           >  2:205744/1‑101 (MQ=255)
                        aTCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAgg                                                                         <  1:4471309/101‑1 (MQ=255)
                        aTCGCCCGTCACATCGGCCCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAgg                                                                         <  2:734979/101‑1 (MQ=255)
                          cGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGc                                                                       >  1:3190134/1‑101 (MQ=255)
                           gCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCa                                                                      <  2:3885931/101‑1 (MQ=255)
                            cccGTCACATCGGCGCCGGTGACGAATTGCCCCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCAc                                                                     <  1:177199/101‑1 (MQ=255)
                            cccGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCAc                                                                     <  1:2264002/101‑1 (MQ=255)
                            cccGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCAc                                                                     <  1:2076937/101‑1 (MQ=255)
                            cccGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCAc                                                                     <  1:1236880/101‑1 (MQ=255)
                             ccGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGGGGCGCAGGCCCAgggggcg                                                                    >  1:1770273/1‑97 (MQ=255)
                             ccGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCAcc                                                                    >  2:900844/1‑100 (MQ=255)
                             ccGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACg                                                                    <  1:227591/101‑1 (MQ=255)
                             ccGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACg                                                                    <  2:4042343/101‑1 (MQ=255)
                              cGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCCCGc                                                                   >  1:169632/1‑101 (MQ=255)
                               gTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGGACGcc                                                                  >  2:3292295/1‑101 (MQ=255)
                               gTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGACAGCGGCAGCCCGGCCGCCGCGGCGGCCGGGCCAGCCGAGGGCACggc                                                                  >  1:3573658/1‑99 (MQ=255)
                                tCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCg                                                                 <  2:248368/101‑1 (MQ=255)
                                    aTCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGGGCAGGCCGAGGGCACGCCggggg                                                             >  1:3593376/1‑100 (MQ=255)
                                    aTCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGGACGCCGGGga                                                             >  2:4398648/1‑101 (MQ=255)
                                    aTCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCggggg                                                             >  2:2114436/1‑100 (MQ=255)
                                      cGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGGGAGGGCCGAGGGCACGCCGGGgagg                                                           >  2:1841761/1‑99 (MQ=255)
                                       ggCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGgatga                                                          >  1:1583076/1‑101 (MQ=255)
                                             ggTGACGAATTGCCCCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGt                                                    <  1:98948/101‑1 (MQ=255)
                                              gTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTa                                                   >  2:2450440/1‑101 (MQ=255)
                                              gTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCCTa                                                   >  1:595643/1‑101 (MQ=255)
                                               tGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAg                                                  >  2:642539/1‑101 (MQ=255)
                                               tGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGAGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAg                                                  >  2:1727434/1‑101 (MQ=255)
                                                gACGAATTGCCCCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGc                                                 <  2:4612056/101‑1 (MQ=255)
                                                gACGAATTGCACCCGCCCCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGc                                                 <  2:4260000/101‑1 (MQ=255)
                                                  cGAATTGCCCCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCg                                               <  2:3561005/101‑1 (MQ=255)
                                                  cGAATGGCCCCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCg                                               <  2:275697/101‑1 (MQ=255)
                                                   gAATTGCCCCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGa                                              <  2:549172/101‑1 (MQ=255)
                                                   gAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGa                                              >  2:2047753/1‑101 (MQ=255)
                                                    aaTTGCGCCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGAt                                             <  2:708142/101‑1 (MQ=255)
                                                    aaTTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGAt                                             >  1:3725554/1‑101 (MQ=255)
                                                     aTTGCCCCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATg                                            <  2:3857010/101‑1 (MQ=255)
                                                       tGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGcc                                          >  2:1334655/1‑101 (MQ=255)
                                                       tGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGcc                                          <  2:4698079/101‑1 (MQ=255)
                                                        gcccccGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGccc                                         <  2:2981800/98‑1 (MQ=255)
                                                        gCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCTGCGGCGCATGTCTAGGGCACGCCGGGGATGATCCCGTATCCGATGcac                                         >  2:3130596/1‑99 (MQ=255)
                                                         cACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCcg                                        >  1:2584995/1‑101 (MQ=255)
                                                         cACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCcg                                        >  1:2141127/1‑101 (MQ=255)
                                                          aCCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCcgc                                       >  1:1674420/1‑101 (MQ=255)
                                                            ccgccgCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCgcgc                                     <  2:2258278/101‑1 (MQ=255)
                                                             cgccgcGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGcc                                    <  1:1987454/101‑1 (MQ=255)
                                                              gccgcGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCg                                   <  1:595422/101‑1 (MQ=255)
                                                              gccgcGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACCCCGGGGATGATCTCGTAGGCGATGCCCGCGCCg                                   >  1:4094028/1‑101 (MQ=255)
                                                               ccgcGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCgg                                  <  2:3255260/101‑1 (MQ=255)
                                                               ccgcGCGGGCAGACGGGGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCgg                                  <  1:1003286/101‑1 (MQ=255)
                                                                cgcgcgGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGt                                 <  2:1616861/101‑1 (MQ=255)
                                                                cgcgcgGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGt                                 <  1:4178726/101‑1 (MQ=255)
                                                                 gcgcgGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTc                                <  1:107703/101‑1 (MQ=255)
                                                                 gcgcgGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGGAGCCGCAGCCCGCGCCgggc                                >  1:1264968/1‑99 (MQ=255)
                                                                  cgcgGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCa                               <  1:158153/101‑1 (MQ=255)
                                                                  cgcgGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGCCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCa                               <  2:3915446/101‑1 (MQ=255)
                                                                    cgGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGGTGCCCGCGCCGGTCAgc                             >  2:1147003/1‑101 (MQ=255)
                                                                    cgGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAgc                             <  2:3455766/101‑1 (MQ=255)
                                                                    cgGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGGCAgc                             >  1:4091858/1‑101 (MQ=255)
                                                                    cgGTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCAGAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAgc                             >  2:4105373/1‑101 (MQ=255)
                                                                     ggTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAgcg                            <  1:1008251/101‑1 (MQ=255)
                                                                     ggTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAgcg                            >  1:3882118/1‑101 (MQ=255)
                                                                      gTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAgcgc                           <  1:4742266/101‑1 (MQ=255)
                                                                      gTCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCAGATGCCCGCGCCGGTCAgcgc                           >  2:137525/1‑101 (MQ=255)
                                                                       tCAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGTCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGcgcg                          <  1:3930534/101‑1 (MQ=255)
                                                                        cAGACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGt                         >  1:74925/1‑101 (MQ=255)
                                                                          gACGGCGGGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTcc                       <  1:3409229/101‑1 (MQ=255)
                                                                            cggcggGTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAg                     <  1:3490700/101‑1 (MQ=255)
                                                                                gggTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGttct                 <  2:3075741/101‑1 (MQ=255)
                                                                                 ggTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGttctt                <  1:3620310/101‑1 (MQ=255)
                                                                                 ggTCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGttctt                <  1:1742625/101‑1 (MQ=255)
                                                                                 ggTCAGCGGCAGCCCGGCCGCGGCGCCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGttctt                <  2:3593157/101‑1 (MQ=255)
                                                                                   tCAGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGTTCTTcc              >  1:3591734/1‑101 (MQ=255)
                                                                                     aGCGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGACCGCCTGGCAGGATCTCGTAGCCGGTGTCCGCGCCGGTCAGCCTGTCCAGTTCTTCCTc            >  2:4414724/1‑101 (MQ=255)
                                                                                       cGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGTTCTTCCTCaa          <  2:1757948/101‑1 (MQ=255)
                                                                                       cGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGTTCTTCCTCaa          <  2:3798845/101‑1 (MQ=255)
                                                                                       cGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGTTCTTCCTCaa          <  1:1908107/101‑1 (MQ=255)
                                                                                       cGGCAGCCCGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGTTCTTCCTCaa          <  2:106212/101‑1 (MQ=255)
                                                                                             cccGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGTTCTTCCTCAAgccggc    >  1:4464818/1‑101 (MQ=255)
                                                                                             cccGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGTTCTCCCTCAAgccggc    >  2:285135/1‑101 (MQ=255)
                                                                                              ccGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGTTCTTCCTCAAgccggcc   >  2:2617311/1‑101 (MQ=255)
                                                                                               cGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGTTCTTCCTCAAgccggccg  >  1:1112624/1‑101 (MQ=255)
                                                                                               cGGCCGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGTTCTTCCTCAAgccggccg  >  1:836391/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
CGCCCAGTTCAGCCCCTCGGGCAGATCGCCCGTCACATCGGCGCCGGTGACGAATTGCACCCGCCGCGCGGTCAGACGGCGGGTCAGCGGCAGCCCGGACGCGGCGGCGGCGGCGCAGGCCGAGGGCACGCCGGGGATGATCTCGTAGCCGATGCCCGCGCCGGTCAGCGCGTCCAGTTCTTCCTCAAGCCGGCCG  >  NC_014034/2198389‑2198584

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 40 ≤  ATCG