Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 735,233 +C 2.4% intergenic (‑139/‑346) bchF ← / → ppaA 2‑vinyl bacteriochlorophyllide hydratase/regulatory protein PpaA

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034735,2331.→C 2.4% 7.1 ‑1.7 ‑0.0209intergenic (‑139/‑346)bchF/ppaA2‑vinyl bacteriochlorophyllide hydratase/regulatory protein PpaA
Reads supporting (aligned to +/- strand):  new base (1/4):  ref base (115/89):  total (116/93)
Fisher's exact test strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.36e-01

Only 100 of 249 total aligned reads displayed.

TAAAACTGGTCTTGTCACGAGAAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTGGCCCGTATCTTGCGTCGTGGC  >  NC_014034/735133‑735332
                                                                                                     |                                                                                                   
tAAAACTGGTCTTGTCACGAGAAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCgg                                                                                                      >  2:1036023/1‑101 (MQ=255)
tAAAACTGGTCTTGTCACGAGAAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCgg                                                                                                      >  1:2052818/1‑101 (MQ=255)
 aaaaCTGGTCTTGTCACGAGAAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.c                                                                                                    >  1:4082414/1‑101 (MQ=255)
    aCTGGTCTTGTCACGAGAAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.cccc                                                                                                 <  1:2578620/101‑1 (MQ=255)
      tGGTCTTGTCACGAGAAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.cccccc                                                                                               >  1:3004827/1‑101 (MQ=255)
        gTCTTGTCACGAGAAAGTGTCAATGAAAACTTATACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.cccccccc                                                                                             <  1:3803469/101‑1 (MQ=255)
        gTCTTGTCACGAGAAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.cccccccc                                                                                             <  1:3156134/101‑1 (MQ=255)
          cTTGTCACGAGAAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCg                                                                                           >  2:366308/1‑101 (MQ=255)
            tGTCACGAGAAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGc                                                                                         >  2:2411475/1‑101 (MQ=255)
            tGTCACGAGAAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGc                                                                                         >  1:1607295/1‑101 (MQ=255)
             gTCACGAGAAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCt                                                                                        <  1:3023849/101‑1 (MQ=255)
             gTCACGAGAAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCt                                                                                        <  1:1570552/101‑1 (MQ=255)
                  gagaAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAAcc                                                                                   >  2:236661/1‑101 (MQ=255)
                    gaAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACcgc                                                                                 >  2:2788140/1‑101 (MQ=255)
                    gaAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACcgc                                                                                 <  2:876026/101‑1 (MQ=255)
                    gaAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACcgc                                                                                 >  1:1045499/1‑101 (MQ=255)
                     aaaGTGTCAATGAAAACTTACCCTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTCATGCAGCGG.CCCCCCCCCGGCTAAACcgcg                                                                                <  1:3970741/101‑1 (MQ=255)
                       aGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGt                                                                              <  2:59139/101‑1 (MQ=255)
                        gtgtCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTc                                                                             <  1:1741519/101‑1 (MQ=255)
                         tgtCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCa                                                                            >  1:3820309/1‑101 (MQ=255)
                         tgtCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCa                                                                            <  2:3741392/101‑1 (MQ=255)
                         tgtCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCa                                                                            >  1:1070630/1‑101 (MQ=255)
                         ggtCAAGAAAAACTTACCCTCAGGGGCGGGGCGAAACGCCCCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGGCCCCCCCCCCGGCTAAACCGCGGTc                                                                             <  1:347412/100‑1 (MQ=255)
                           tCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACg                                                                          <  2:1409706/101‑1 (MQ=255)
                             aaTGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGAc                                                                        >  1:540056/1‑101 (MQ=255)
                               tGAAAACTTTCACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTc                                                                      <  2:2911570/101‑1 (MQ=255)
                               tGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTCAACCGCGGGCACGACTc                                                                      >  2:110257/1‑101 (MQ=255)
                               ngANANNTTACACGCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTc                                                                      <  1:1822380/100‑1 (MQ=255)
                                gAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCg                                                                     <  2:2684276/101‑1 (MQ=255)
                                gAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCg                                                                     >  2:94930/1‑101 (MQ=255)
                                gAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCg                                                                     >  2:1814614/1‑101 (MQ=255)
                                 aaaaCTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGt                                                                    <  1:2612163/101‑1 (MQ=255)
                                   aaCTTACCCTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTcc                                                                  <  1:1722925/101‑1 (MQ=255)
                                    aCTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCg                                                                 <  2:3049085/101‑1 (MQ=255)
                                    aCTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCg                                                                 <  2:22541/101‑1 (MQ=255)
                                    aCTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCg                                                                 <  1:1786977/101‑1 (MQ=255)
                                     cTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCgg                                                                >  1:3355300/1‑101 (MQ=255)
                                     cTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCgg                                                                >  2:451184/1‑101 (MQ=255)
                                      ttACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCggg                                                               >  1:2012637/1‑101 (MQ=255)
                                         cacTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGCCTCGTCCGGGCAg                                                            >  1:2186349/1‑101 (MQ=255)
                                         cacTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCCGGCAg                                                            >  1:3513379/1‑101 (MQ=255)
                                         cacTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG..CCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGa                                                           >  1:1441850/1‑101 (MQ=255)
                                          acTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGCCACGACTCCTCCGGGCGGa                                                           >  1:2321238/1‑101 (MQ=255)
                                           cTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGaa                                                          >  1:73868/1‑101 (MQ=255)
                                           cTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGaa                                                          >  2:113531/1‑101 (MQ=255)
                                            tCAGGGGCGGGACGAACCGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCGGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGaaa                                                         <  2:4400009/101‑1 (MQ=255)
                                            tCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGaaa                                                         <  2:2564555/101‑1 (MQ=255)
                                            tCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGaaa                                                         <  1:3278345/101‑1 (MQ=255)
                                            tCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACCCGTCCGGGCAGaaa                                                         >  1:1910204/1‑101 (MQ=255)
                                            tCAGGGGCGGGACCAAACGCCGCTGCCCTTGCGCGCAAGGGCTCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGaaa                                                         <  2:416669/101‑1 (MQ=255)
                                              aGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAAcc                                                       >  2:795372/1‑101 (MQ=255)
                                              aGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAAcc                                                       >  2:3862979/1‑101 (MQ=255)
                                               ggggCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACcg                                                      <  1:3901194/101‑1 (MQ=255)
                                                   cGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACcgggcg                                                  >  1:1209660/1‑101 (MQ=255)
                                                   cGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGACCGGGCCGAAACcgggcg                                                  >  1:895755/1‑101 (MQ=255)
                                                    gggACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACcgggcgg                                                 >  1:2368469/1‑101 (MQ=255)
                                                     ggACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACcgggcggg                                                <  2:2320082/101‑1 (MQ=255)
                                                     ggACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACcgggcggg                                                <  2:1283279/101‑1 (MQ=255)
                                                      gACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCgggcaggg                                               <  2:979151/101‑1 (MQ=255)
                                                      gACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGGCACGACTCGGCCGGGCAGAAACCgggcgggg                                               >  1:4414979/1‑101 (MQ=255)
                                                         gAAACTCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGAAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACCACTCTCCCGGCCAGAAACCGggcggggcgg                                            >  1:1990349/1‑101 (MQ=255)
                                                         gAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGggcggggcgg                                            >  1:4049611/1‑101 (MQ=255)
                                                          aaaCGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGgcggggcggg                                           >  2:3293798/1‑101 (MQ=255)
                                                          aaaCGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAAAAACCGGgcgggggggg                                           >  1:3693265/1‑101 (MQ=255)
                                                            aCGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGGCCCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGa                                          <  2:197373/101‑1 (MQ=255)
                                                            aCGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG..CCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATg                                        <  1:3714517/101‑1 (MQ=255)
                                                             cgccgcTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATg                                        >  1:607288/1‑101 (MQ=255)
                                                               ccgcTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAAAAACCGGGCGGGGCGGGATgcg                                      >  2:2585571/1‑101 (MQ=255)
                                                                 gcTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGGGGCGcc                                    >  1:4392827/1‑101 (MQ=255)
                                                                 gcTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGGAGCGcc                                    >  1:4326659/1‑101 (MQ=255)
                                                                  cTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGGGGGGCGGGATGCGCCt                                   >  1:2895422/1‑101 (MQ=255)
                                                                    gCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGCTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCttt                                 <  2:1340326/101‑1 (MQ=255)
                                                                     cccTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCtttt                                >  2:2090218/1‑101 (MQ=255)
                                                                     cccTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCtttt                                <  1:3370078/101‑1 (MQ=255)
                                                                      ccTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTtc                               <  2:3434975/101‑1 (MQ=255)
                                                                       cgtGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTtct                              <  2:4362569/99‑1 (MQ=255)
                                                                       cTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCCGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGGTGCGACtttttc                              >  2:3270342/1‑99 (MQ=255)
                                                                         tGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCa                            <  2:1534128/101‑1 (MQ=255)
                                                                         tGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCa                            >  1:980706/1‑101 (MQ=255)
                                                                          gcgcgcAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGGGCCTTTTCCCaa                           >  1:4293344/1‑101 (MQ=255)
                                                                            gcgcAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAAcc                         <  2:1233492/101‑1 (MQ=255)
                                                                              gcAGGGTCGCCTTGGTGCAGCGG.CCCTCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGAGGCGCCTTTTCTCAACCTg                       <  2:736208/101‑1 (MQ=255)
                                                                              gcAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTg                       <  2:3359419/101‑1 (MQ=255)
                                                                               cAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTgg                      >  2:4203655/1‑101 (MQ=255)
                                                                                  gggCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGGGCCTTTTCTCAACCTGGccc                   >  1:2236185/1‑101 (MQ=255)
                                                                                  gggCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTGGccc                   <  1:1175592/101‑1 (MQ=255)
                                                                                   ggCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTGGCCCg                  <  2:1733221/101‑1 (MQ=255)
                                                                                    gcgcCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGGGCCTTTTCTCAACCCGGCCCGt                 >  1:1346845/1‑101 (MQ=255)
                                                                                    gcgcCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTGGCCCGt                 >  2:3558374/1‑101 (MQ=255)
                                                                                     cgcCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTGGCCCGTa                <  1:1278356/101‑1 (MQ=255)
                                                                                      gcCTTGCTGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTGGCCCGTAt               <  2:1607189/101‑1 (MQ=255)
                                                                                        cTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGGTGCGCCTTTTCTCAACCTGGCCCCGATCt             >  1:3183124/1‑101 (MQ=255)
                                                                                        cTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTGGCCCGTATCt             <  2:1312270/101‑1 (MQ=255)
                                                                                         ttGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGGTGCGCCTTTTCTCAACCTGGCCCGTATCtt            >  2:3138640/1‑101 (MQ=255)
                                                                                          tGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTGGCCCGTATCTTg           >  2:2879007/1‑101 (MQ=255)
                                                                                          tGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGAGGCGCCTTTTCTCAACCTGGCCCGTATCTTg           >  1:826863/1‑101 (MQ=255)
                                                                                           gATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTGGCCCGTATCTTGc          <  1:3294030/101‑1 (MQ=255)
                                                                                           gATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTGGCCCGTATCTTGc          <  1:2851983/101‑1 (MQ=255)
                                                                                            aTGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCCTTTCTCAACCTGGCCCCTCTCTTGcg         >  2:4437381/1‑101 (MQ=255)
                                                                                               cAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTGGCCCGTATCTTGcgtcg      <  2:2012508/101‑1 (MQ=255)
                                                                                                     |                                                                                                   
TAAAACTGGTCTTGTCACGAGAAAGTGTCAATGAAAACTTACACTCAGGGGCGGGACGAAACGCCGCTGCCCTTGCGCGCAAGGGCGCCTTGATGCAGCGG.CCCCCCCCCGGCTAAACCGCGGTCACGACTCGTCCGGGCAGAAACCGGGCGGGGCGGGATGCGCCTTTTCTCAACCTGGCCCGTATCTTGCGTCGTGGC  >  NC_014034/735133‑735332

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG