Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,398,015 A→C 1.5% intergenic (+319/+1579) RCAP_rcc01307 → / ← RCAP_rcc01309 hypothetical protein/IS66 family transposase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,398,0150A→C 1.5% 2.1 0.0 ‑0.7265intergenic (+319/+1579)RCAP_rcc01307/RCAP_rcc01309hypothetical protein/IS66 family transposase
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (137/124):  total (139/126)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.00e-01

Only 100 of 300 total aligned reads displayed.

GCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAA  >  NC_014034/1397915‑1398115
                                                                                                    |                                                                                                       
 cgGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAg                                                                                                        <  1:5134851/101‑1 (MQ=255)
  gggggCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAACAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGc                                                                                                       >  2:2088068/1‑101 (MQ=255)
   ggggCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGcc                                                                                                      >  2:2902800/1‑101 (MQ=255)
   ggggCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGcc                                                                                                      >  2:1129772/1‑101 (MQ=255)
    gggCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTTCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGccc                                                                                                     <  2:477408/101‑1 (MQ=255)
    gggCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGccc                                                                                                     <  1:3970796/101‑1 (MQ=255)
       cattgATCGAAACTCTATGGGAGGCTAGGGAGCCTACGCGCCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACCGCCGAGCAAAAAGCCCGcc                                                                                                  <  2:5711432/97‑1 (MQ=37)
       cGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGcc                                                                                                  >  1:4334448/1‑101 (MQ=255)
        gTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCa                                                                                                 <  2:5107129/101‑1 (MQ=255)
        gTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCa                                                                                                 <  1:3072198/101‑1 (MQ=255)
          ggATCGAAACTCTATGGGAGGCTAGGGGGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCCCCaag                                                                                               >  1:926616/1‑101 (MQ=255)
          ggATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGACCGCCaag                                                                                               >  1:3429815/1‑101 (MQ=255)
             tCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCaagaag                                                                                            <  1:4256939/101‑1 (MQ=255)
              cGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGc                                                                                           <  2:1062787/101‑1 (MQ=255)
               gAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCg                                                                                          >  2:5597232/1‑101 (MQ=255)
                gatctctATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCgg                                                                                         <  1:1969606/98‑1 (MQ=255)
                aaaCTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCgg                                                                                         <  1:1800518/101‑1 (MQ=255)
                 aaCTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..a                                                                                      <  1:4466876/101‑1 (MQ=255)
                 aaCTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..a                                                                                      >  2:4015382/1‑101 (MQ=255)
                 aaCTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..a                                                                                      >  1:1776610/1‑101 (MQ=255)
                    tctATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCa                                                                                   <  1:1188626/101‑1 (MQ=255)
                    tctATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCa                                                                                   >  2:4156431/1‑101 (MQ=255)
                     ctATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCa.a                                                                                 <  2:3263659/101‑1 (MQ=255)
                      tATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.Ag                                                                                <  2:1197715/101‑1 (MQ=255)
                      tATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.Ag                                                                                >  1:3864434/1‑101 (MQ=255)
                       aTGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGt                                                                               >  1:854999/1‑101 (MQ=255)
                       aTGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGt                                                                               <  2:5669236/101‑1 (MQ=255)
                         gggaggctagggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCg                                                                             >  1:3647061/1‑101 (MQ=255)
                         gggaggctagggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCg                                                                             <  2:5559986/101‑1 (MQ=255)
                             ggctagggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTcc                                                                         >  2:2970068/1‑101 (MQ=255)
                               ctagggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTcccc                                                                       >  2:44793/1‑101 (MQ=255)
                               ctagggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTcccc                                                                       <  1:2280731/101‑1 (MQ=255)
                                tagggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTccccc                                                                      <  2:2011962/101‑1 (MQ=255)
                                tagggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTccccc                                                                      >  1:5552466/1‑101 (MQ=255)
                                tagggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTccccc                                                                      >  1:2114681/1‑101 (MQ=255)
                                 agggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTcccccc                                                                     >  1:766962/1‑101 (MQ=255)
                                 agggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTcccccc                                                                     <  1:735557/101‑1 (MQ=255)
                                 agggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTcccccc                                                                     <  2:584119/101‑1 (MQ=255)
                                  gggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCcg                                                                    >  1:160378/1‑101 (MQ=255)
                                  gggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCcg                                                                    >  2:1942128/1‑101 (MQ=255)
                                   ggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCcgc                                                                   >  2:4655791/1‑101 (MQ=255)
                                     gggctaCGCCCCTTCAAAAGTCGCTCCCCCCGCGGGGCCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCgg                                                                 <  2:1367250/100‑1 (MQ=255)
                                      ggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCggg                                                                <  1:1511433/101‑1 (MQ=255)
                                         taCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGc                                                             >  1:5700832/1‑101 (MQ=255)
                                         taCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGc                                                             >  2:5519230/1‑101 (MQ=255)
                                            gCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTg                                                          <  1:4736493/101‑1 (MQ=255)
                                             cACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTgg                                                         <  1:4677303/101‑1 (MQ=255)
                                              cccTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGa                                                        <  2:5361984/100‑1 (MQ=255)
                                               ccTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGAt                                                       <  1:4828444/101‑1 (MQ=255)
                                               ccTTCAAAAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGAt                                                       <  1:4159940/101‑1 (MQ=255)
                                                cTTCAAATGGCGCTCCCCCCGCGGGGGCGGGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTTGATc                                                      <  2:3734870/101‑1 (MQ=38)
                                                cTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAGCCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATc                                                      <  2:3429610/101‑1 (MQ=255)
                                                cTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATc                                                      >  1:4065089/1‑101 (MQ=255)
                                                cTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATc                                                      <  2:5650296/101‑1 (MQ=255)
                                                 ttCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                     >  2:2746258/1‑101 (MQ=255)
                                                  tgcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAA.CGCTGCTGACACCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATc                                                      <  2:3545145/95‑1 (MQ=2)
                                                  tgcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAA.CGCTGCTGACACCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATc                                                      <  1:2512478/95‑1 (MQ=2)
                                                   cAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGGGGATCGTa                                                   >  1:2599025/1‑101 (MQ=255)
                                                    aaaaGCCGCTCCCCCCGCGGGGGCGTGGATCAAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTaa                                                  <  1:1062868/101‑1 (MQ=255)
                                                          cGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTg                                            >  2:796321/1‑101 (MQ=255)
                                                            cTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCa                                          >  1:3250909/1‑101 (MQ=255)
                                                              ccccccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATc                                        <  2:3504729/101‑1 (MQ=255)
                                                              ccccccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATc                                        >  2:2456415/1‑101 (MQ=255)
                                                                 cccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGc                                     >  2:92177/1‑101 (MQ=255)
                                                                  ccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGcc                                    >  2:2147088/1‑101 (MQ=255)
                                                                  ccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGcc                                    <  2:4359007/101‑1 (MQ=255)
                                                                    gcgGGGGCGTGGATCGAACCCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCg                                  <  1:1024253/101‑1 (MQ=255)
                                                                      gggggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGAt                                <  1:5173747/101‑1 (MQ=255)
                                                                      gggggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGAt                                <  1:5403303/101‑1 (MQ=255)
                                                                       ggggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATa                               >  2:1772302/1‑101 (MQ=255)
                                                                       ctggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATa                               <  1:1240763/99‑1 (MQ=255)
                                                                        gggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGGGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGAGGATCGTAAGCGGTGCCTCGGCCCGAGaa                              >  2:3059849/1‑101 (MQ=39)
                                                                        gggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATaa                              >  2:482868/1‑101 (MQ=255)
                                                                          gCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAAc                            <  1:5297224/101‑1 (MQ=255)
                                                                          gCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAAc                            <  1:5647062/101‑1 (MQ=255)
                                                                           cGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAAcc                           <  2:5256289/101‑1 (MQ=255)
                                                                            gTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAAccc                          <  2:1859343/101‑1 (MQ=255)
                                                                               gATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGa                       >  1:3171510/1‑101 (MQ=255)
                                                                               gATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGa                       >  1:2389469/1‑101 (MQ=255)
                                                                                aTCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCCAt                      >  2:4416878/1‑101 (MQ=255)
                                                                                  cGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGAttt                    <  1:3919952/101‑1 (MQ=255)
                                                                                   gAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGAtttt                   >  1:2447907/1‑101 (MQ=255)
                                                                                   gAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGAtttt                   >  1:2164931/1‑101 (MQ=255)
                                                                                   gAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGAtttt                   <  2:221605/101‑1 (MQ=255)
                                                                                   gAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGGAAGCGGTGCATCGGCCCGATAAAACCCGAtttt                   >  1:4931815/1‑101 (MQ=255)
                                                                                     acccgccgCCCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCg                 <  1:2902967/99‑1 (MQ=255)
                                                                                     aaCCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCg                 <  2:4691261/101‑1 (MQ=255)
                                                                                            gcgaAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTg          <  2:5456310/101‑1 (MQ=255)
                                                                                             cgcAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGa         <  2:4688055/101‑1 (MQ=255)
                                                                                               cAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAgg       <  2:504374/101‑1 (MQ=255)
                                                                                               cAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAgg       >  2:5186830/1‑101 (MQ=255)
                                                                                               cAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAgg       <  1:5505232/101‑1 (MQ=255)
                                                                                               cAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAgg       <  1:1117954/101‑1 (MQ=255)
                                                                                                caaaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAggg      <  2:2560171/100‑1 (MQ=255)
                                                                                                aaaaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCTGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAggg      <  2:3962217/101‑1 (MQ=255)
                                                                                                aaaaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAAGCCCGGTTTTCGGCCCCTGAggg      >  2:3769469/1‑101 (MQ=255)
                                                                                                 aaaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAgggg     >  2:3346866/1‑101 (MQ=255)
                                                                                                 aaaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGCGGATCGTAGGCGGTGCATTTGCCCGTTAAACCCCGCTTTTCGCCCCCTGAgggg     >  2:2658518/1‑101 (MQ=37)
                                                                                                  aaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGa    <  1:2785032/101‑1 (MQ=255)
                                                                                                    aGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGaaa  <  2:1523487/101‑1 (MQ=255)
                                                                                                    |                                                                                                       
GCGGGGGCGTGGATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAA  >  NC_014034/1397915‑1398115

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG