Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,398,029 C→T 1.5% intergenic (+333/+1565) 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,0290C→T 1.5% 2.1 0.0 ‑0.9272intergenic (+333/+1565)RCAP_rcc01307/RCAP_rcc01309hypothetical protein/IS66 family transposase
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (132/135):  total (135/137)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.33e-01

Only 100 of 293 total aligned reads displayed.

CGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCTTGGCACAA  >  NC_014034/1397929‑1398129
                                                                                                    |                                                                                                       
  aaaCTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCgg                                                                                                       <  1:4031122/101‑1 (MQ=255)
  aaaCTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCgg                                                                                                       <  1:1800518/101‑1 (MQ=255)
      tctATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCa                                                                                                 >  1:4329449/1‑101 (MQ=255)
      tctATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCa                                                                                                 >  2:4736212/1‑101 (MQ=255)
        tATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.Ag                                                                                              >  1:3864434/1‑101 (MQ=255)
        tATGGGAGCCTAGGGAAGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.Ag                                                                                              <  2:4437200/101‑1 (MQ=255)
          tGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTc                                                                                            >  2:1829977/1‑101 (MQ=255)
          tGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTc                                                                                            <  1:4722982/101‑1 (MQ=255)
          tGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTc                                                                                            >  2:5587591/1‑101 (MQ=255)
                  tagggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTccccc                                                                                    >  1:5552466/1‑101 (MQ=255)
                  tagggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTccccc                                                                                    <  2:2011962/101‑1 (MQ=255)
                   agggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTcccccc                                                                                   <  2:402247/101‑1 (MQ=255)
                   agggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTcccccc                                                                                   >  1:2310153/1‑101 (MQ=255)
                   agggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTcccccc                                                                                   <  2:584119/101‑1 (MQ=255)
                    gggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCcg                                                                                  >  2:1942128/1‑101 (MQ=255)
                     ggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCcgc                                                                                 >  2:5687864/1‑101 (MQ=255)
                     ggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCcgc                                                                                 >  2:4655791/1‑101 (MQ=255)
                      gaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCCAAAAGCCCGCCAAGAAGCGG..AGCA.AGGCGCTCCCCCcgcg                                                                                >  2:1891482/1‑101 (MQ=255)
                       aggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCgg                                                                               <  1:2172995/101‑1 (MQ=255)
                        ggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCggg                                                                              <  1:1511433/101‑1 (MQ=255)
                          ctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCggggg                                                                            >  2:439287/1‑101 (MQ=255)
                          ctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCggggg                                                                            >  1:3221633/1‑101 (MQ=255)
                           taCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGc                                                                           >  1:5700832/1‑101 (MQ=255)
                           taCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGc                                                                           >  1:4782577/1‑101 (MQ=255)
                           taCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGc                                                                           >  1:1100677/1‑101 (MQ=255)
                           taCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGc                                                                           >  1:3962461/1‑101 (MQ=255)
                           taCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGc                                                                           <  1:1583655/101‑1 (MQ=255)
                                aCCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGa                                                                      >  2:5214270/1‑101 (MQ=255)
                                 ccTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                     <  2:5667172/101‑1 (MQ=255)
                                  cTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAGCCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                    <  2:3429610/101‑1 (MQ=255)
                                  cTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                    >  2:3743021/1‑101 (MQ=255)
                                   ttCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                   >  2:2746258/1‑101 (MQ=255)
                                    tgcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAA.CGCTGCTGACACCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                    >  2:4972853/7‑101 (MQ=2)
                                    tgcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAA.CGCTGCTGACACCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                    >  1:1775134/7‑101 (MQ=2)
                                    tgcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAA.CGCTGCTGACACCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                    <  1:2512478/95‑1 (MQ=2)
                                     cAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTa                                                                 >  2:5140408/1‑101 (MQ=255)
                                     cAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGGGGATCGTa                                                                 >  1:2599025/1‑101 (MQ=255)
                                         aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGTGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCg                                                             >  1:5597775/1‑101 (MQ=255)
                                           tCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCggg                                                           >  2:1026179/1‑100 (MQ=255)
                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGc                                                         <  1:3352492/101‑1 (MQ=255)
                                              cTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCa                                                        >  2:5646770/1‑101 (MQ=255)
                                               cccccccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCAt                                                       <  2:2114564/100‑1 (MQ=255)
                                                ccccccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATc                                                      <  1:1949584/101‑1 (MQ=255)
                                                    ccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGcc                                                  >  2:2147088/1‑101 (MQ=255)
                                                     cgcgGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGccc                                                 <  1:5136228/101‑1 (MQ=255)
                                                      gcgGGGGCGTGGATCGAACCCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCg                                                <  1:1024253/101‑1 (MQ=255)
                                                        gggggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGAt                                              >  2:3511591/1‑101 (MQ=255)
                                                        gggggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGAt                                              >  1:76050/1‑101 (MQ=255)
                                                        gggggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGAt                                              <  1:5403303/101‑1 (MQ=255)
                                                        gggggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGAt                                              >  1:2422232/1‑101 (MQ=255)
                                                        gggggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGAt                                              <  1:5173747/101‑1 (MQ=255)
                                                         gtgtcGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATa                                             <  2:5531245/97‑1 (MQ=255)
                                                          gggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATaa                                            <  2:3025197/101‑1 (MQ=255)
                                                          gggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATaa                                            <  1:4423123/101‑1 (MQ=255)
                                                           ggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATaag                                           >  2:1240680/1‑100 (MQ=255)
                                                            gCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAAc                                          <  1:5297224/101‑1 (MQ=255)
                                                              gTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGCCCCGATAAAccc                                        <  1:3142571/101‑1 (MQ=255)
                                                               tGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAAcccc                                       <  2:5287733/101‑1 (MQ=255)
                                                                ggATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCg                                      <  2:3450872/101‑1 (MQ=255)
                                                                   tCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGAtt                                   <  1:2970235/101‑1 (MQ=255)
                                                                    cGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGAttt                                  <  1:3919952/101‑1 (MQ=255)
                                                                     gAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGAtttt                                 >  1:881922/1‑101 (MQ=255)
                                                                     gAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGAtttt                                 <  2:221605/101‑1 (MQ=255)
                                                                     gAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGAtttt                                 <  1:2088185/101‑1 (MQ=255)
                                                                     gAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACACCGAtttt                                 <  1:4218316/101‑1 (MQ=255)
                                                                     gAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGGAAGCGGTGCATCGGCCCGATAAAACCCGAtttt                                 >  1:4931815/1‑101 (MQ=255)
                                                                      aaaCCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGCTTTTc                                >  2:3320016/1‑101 (MQ=255)
                                                                      aaaCCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTc                                <  2:4782296/101‑1 (MQ=255)
                                                                       acccgccgCCCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCg                               <  1:2902967/99‑1 (MQ=255)
                                                                       aaCCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.ATTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGAttttag                               <  2:2962563/101‑3 (MQ=255)
                                                                       aaCCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCg                               <  2:4691261/101‑1 (MQ=255)
                                                                       aaCCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCg                               >  1:698187/1‑101 (MQ=255)
                                                                         ccgccgCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGc                             >  1:2115626/1‑101 (MQ=255)
                                                                          cgccgcGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGcc                            >  2:1427812/1‑101 (MQ=255)
                                                                             cgcgCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCt                         >  2:357529/1‑101 (MQ=255)
                                                                              gcgcAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTg                        <  2:3967236/101‑1 (MQ=255)
                                                                              gcgcAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTg                        <  2:854930/101‑1 (MQ=255)
                                                                              gcgaAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTg                        <  2:5456310/101‑1 (MQ=255)
                                                                                 cAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAgg                     <  1:5505232/101‑1 (MQ=255)
                                                                                 cAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAgg                     <  2:679336/101‑1 (MQ=255)
                                                                                 cAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAgg                     <  2:4156830/101‑1 (MQ=255)
                                                                                  caaaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAggg                    <  2:2560171/100‑1 (MQ=255)
                                                                                  aaaaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGCggg                    >  2:2585579/1‑101 (MQ=255)
                                                                                  aaaaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAggg                    <  1:1046921/101‑1 (MQ=255)
                                                                                    aaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGa                  <  1:2785032/101‑1 (MQ=255)
                                                                                    aaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGa                  >  1:1199097/1‑101 (MQ=255)
                                                                                     aaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGaa                 <  1:2453068/101‑1 (MQ=255)
                                                                                      aGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGaaa                >  2:1571752/1‑101 (MQ=255)
                                                                                      aGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGaaa                <  2:1122667/101‑1 (MQ=255)
                                                                                       gCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAg               <  2:1351735/101‑1 (MQ=255)
                                                                                        cccGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGa              <  2:2463380/101‑1 (MQ=255)
                                                                                          cGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGAtt            >  1:3801362/1‑101 (MQ=255)
                                                                                          cGCCAAGAAGCGG..AGCA.AGTCGCCCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGAtt            <  1:796389/101‑1 (MQ=255)
                                                                                             cAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCtt         <  2:3250709/101‑1 (MQ=255)
                                                                                             cAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCtt         <  1:1129853/101‑1 (MQ=255)
                                                                                               agaagCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCTTgg       <  1:1274917/101‑1 (MQ=255)
                                                                                                gaagCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCTTGGc      <  1:4673582/101‑1 (MQ=255)
                                                                                                 aagCGG..AGCC.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCTTGGca     <  2:4286140/101‑1 (MQ=255)
                                                                                                  agCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCTTGGcac    >  2:3940353/1‑101 (MQ=255)
                                                                                                    cGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCTTGGCACaa  <  1:3599907/101‑1 (MQ=255)
                                                                                                    |                                                                                                       
CGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCTTGGCACAA  >  NC_014034/1397929‑1398129

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG