Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,820,393 +C 2.1% intergenic (+84/‑102) RCAP_rcc01683 → / → RCAP_rcc01684 terminase‑like family protein/HK97 family phage portal protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,820,3932.→C 2.1% 2.2 ‑3.7 ‑0.2144intergenic (+84/‑102)RCAP_rcc01683/RCAP_rcc01684terminase‑like family protein/HK97 family phage portal protein
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (102/39):  total (102/42)
Fisher's exact test strand distribution p-value = 2.36e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.94e-01

Only 100 of 175 total aligned reads displayed.

AGGTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGAACGG  >  NC_014034/1820293‑1820493
                                                                                                        |                                                                                                       
aGGTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGaaa                                                                                                          >  1:299096/1‑101 (MQ=255)
aGGTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGaaa                                                                                                          >  1:2348253/1‑101 (MQ=255)
aGGTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGaaa                                                                                                          >  1:1901149/1‑101 (MQ=255)
 ggTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..c                                                                                                       <  1:2892928/101‑1 (MQ=255)
 ggTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..c                                                                                                       <  2:2925281/101‑1 (MQ=255)
    gcgcGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..cccc                                                                                                    <  2:687525/101‑1 (MQ=255)
      gcgGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..g                                                                                                <  1:1220373/101‑1 (MQ=255)
        gggCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCt                                                                                              >  1:1383355/1‑101 (MQ=255)
            tttGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTg                                                                                          >  1:3042969/1‑101 (MQ=255)
            tttGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTg                                                                                          >  2:1478767/1‑101 (MQ=255)
              tGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGcc                                                                                        >  1:302584/1‑101 (MQ=255)
              tGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGcc                                                                                        >  1:3630248/1‑101 (MQ=255)
              tGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGcc                                                                                        >  2:1324856/1‑101 (MQ=255)
               gAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGCGCGGGGGTGTGAAA..CCCCC..GCTTCTGCcg                                                                                       >  1:636293/1‑101 (MQ=255)
                  cgcagaaGGGCGTATTTGGACCAGGAAAAA.GTGGCTTGTTTTTTCTT.GGTA.AAAATCCGCCCGGGGGGGGCGGGGGTGTAAAC..CCCCC..GCTTCTGCCGCgg                                                                                    <  2:1844980/96‑1 (MQ=18)
                     gcAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGaaa                                                                                 >  1:1518258/1‑101 (MQ=255)
                     gcAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGaaa                                                                                 >  1:127867/1‑101 (MQ=255)
                     gcAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGaaa                                                                                 >  1:2858774/1‑101 (MQ=255)
                      aaaGGGCGTATTTGGACCAGAAAAAA.GCTGGTTTTTTTTTCTG.GTTA.AAAATACGCCGGGGGGGTGCGGGGGTGTGAAC..CCCCC..GCTTCTGCCGCGGAAAt                                                                                <  2:302539/100‑1 (MQ=21)
                        ggggCTTATTTGGACCAGAAAAAG.GGTGTGTTTTTTTTTGT.GGTG.AAAATGCGCCGGGGGGGTGGGGGGGTGTGAAC..CCCCC..GCTTCTGCCGCGGAAATCt                                                                              <  1:2742653/100‑1 (MQ=255)
                        aGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCt                                                                              >  2:520506/1‑101 (MQ=255)
                        aGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCt                                                                              >  2:2778053/1‑101 (MQ=255)
                            tgTATTTGGACCAGGAAATA.GTTGCGTGTTTTTTCTT.GGTG.AAAATACCCCACGGGGGTGCGGGGATGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTc                                                                          <  2:1929211/100‑1 (MQ=37)
                            cGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTc                                                                          >  2:4021619/1‑101 (MQ=255)
                            cGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTc                                                                          >  1:1951570/1‑101 (MQ=255)
                             gTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCg                                                                         >  1:1697249/1‑101 (MQ=255)
                             gTATTTGGACAAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTA.AAAATACGCCAGGGGGGTGGGGGGGTGTAAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCg                                                                         <  2:333417/101‑1 (MQ=37)
                             gTATTTGGACAAGAAAAAG.GCTGGTTGTTTTTTCTT.GGTG.AAAATCCGCCAGGGGGGGGGGGGGGTGTAAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCg                                                                         <  2:1126590/101‑1 (MQ=25)
                              tATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGa                                                                        >  1:990464/1‑101 (MQ=255)
                               ttttGGGCCAAGAAAAAGGTTGTTTTTTTTTTTTG.GGAA.AAAAGCCGCCGGGGGGGGGGGGGGGTGTAAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGa                                                                        <  1:3968621/100‑1 (MQ=255)
                               ttttGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGGGGGGGTGTGAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGaa                                                                       <  1:2040081/100‑1 (MQ=255)
                               aTTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGaa                                                                       >  2:2457600/1‑101 (MQ=255)
                               aTTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGaa                                                                       >  1:3336254/1‑101 (MQ=255)
                                tttGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGaaa                                                                      >  2:1770409/1‑101 (MQ=255)
                                  tGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATATGGTCGAAAgg                                                                    >  2:1302158/1‑101 (MQ=255)
                                   ggagagagAAAAA.GGTG.GTTTTTTTTTTT.GGGGAAAAAGGCGCGGGGGGGGGGGGGGGGTGGAAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAggg                                                                   <  1:1757427/94‑1 (MQ=255)
                                   agACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGG.GGGAAA..ACCCC..GGCTCTTCCCCCGGAATCTGGGCGGAAgggg                                                                  >  2:846070/2‑100 (MQ=21)
                                    gACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGt                                                                  >  1:3910866/1‑101 (MQ=255)
                                        agaaaaag.GTTGTGTGTTTTTTTTG.GGAA.AAAACCCGCCGGGGGGGGGGGGGGGTGAAAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCggg                                                              <  2:4119883/98‑1 (MQ=255)
                                          gAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGt                                                            >  1:2226376/1‑101 (MQ=255)
                                          gAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGt                                                            >  1:301819/1‑101 (MQ=255)
                                            aaag.GCTGTGTTTTTTTTTTT.GGTA.AAAATCCGCCAGGGGGGTGGGGGGGTGTAAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGttt                                                          <  1:2636724/101‑1 (MQ=255)
                                            aaag.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGttt                                                          <  2:1347390/101‑1 (MQ=255)
                                            aaaa.GCTGGGTGTTTTTCCTT.GGTA.AAAATACGCCAGGGGGGTGCGGGGGTGTGAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGttt                                                          <  2:1243714/101‑1 (MQ=37)
                                            aaaa.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGttt                                                          >  1:3344645/1‑101 (MQ=255)
                                            aaaa.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGGGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGACGAAAGGGTCGGGGttt                                                          >  1:3123613/1‑101 (MQ=255)
                                              aa.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTc                                                        >  1:563274/1‑101 (MQ=255)
                                              aa.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTc                                                        <  2:1037926/101‑1 (MQ=255)
                                               a.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTcc                                                       >  1:236359/1‑101 (MQ=255)
                                                 gCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTccc                                                      <  2:1383255/101‑1 (MQ=255)
                                                 gCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGGTTTccc                                                      >  1:2524494/1‑101 (MQ=255)
                                                 actggggtttttttCTTGGGGG.AAAATCCGCCGGGGGGGGGGGGGGGTGGGAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTcc                                                       <  2:3996397/93‑1 (MQ=255)
                                                  cggggtgtTTTTTCTG.CGCG.AAAATGCGCCAGGGGGGTGGGGGGGTGTGAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCg                                                     <  1:700666/97‑1 (MQ=255)
                                                  cTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGttttagat                                                     >  2:2614883/1‑97 (MQ=255)
                                                    gctttttttttTGT.GGAA.AAAAGCCGCCGGGGGGGTGGGGGGGTGTAAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGc                                                   <  2:3171133/96‑1 (MQ=255)
                                                     cGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGcc                                                  >  2:3198691/1‑101 (MQ=255)
                                                      gtgtTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGCTTTCCCCGccc                                                 >  1:684284/1‑101 (MQ=255)
                                                       ggtTTTTTCTG.GGGA.AAAATCGGCCCGGGGGGGGGGGGGGAGACACC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCt                                                <  1:1134190/100‑1 (MQ=255)
                                                        gtTTTTTCTT.GGTG.AAAATGCGCCACGGGGGCGCGGGGGTGTGAAG..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCtt                                               <  2:636216/101‑1 (MQ=255)
                                                         ttttttCTT.GGTA.AAAATACGCCAGGGGGGTGGGGGGGTGTGAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCttt                                              <  2:397212/101‑1 (MQ=255)
                                                          tttttCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCtttt                                             >  1:3720249/1‑101 (MQ=255)
                                                          tttttCTG.GGTA.AAAATCCGCCAGGGGGGTGGGGGGTGGTGAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCtttt                                             <  2:4234329/101‑1 (MQ=255)
                                                              tgtt.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCtt                                         >  1:2810268/3‑101 (MQ=255)
                                                               cTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCttt                                        >  2:3388510/1‑101 (MQ=255)
                                                               cTT.GGTG.AAAAGCCGCCAGGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCttt                                        <  2:2847342/101‑1 (MQ=255)
                                                                tt.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTg                                       >  2:4071376/1‑101 (MQ=255)
                                                                tt.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GATTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTg                                       >  2:3089792/1‑101 (MQ=255)
                                                                tt.GGTG.AAAATACGCCACGGGGGGGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGACAGGGCCGGGGTTTTCCCGGCCCTTATTCTTTg                                       >  1:729990/1‑101 (MQ=38)
                                                                   gtgg.AAAATACGCCAGGGGGGTGGGGGGGTGTGAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGc                                     <  2:2693061/98‑1 (MQ=255)
                                                                   ggTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGc                                     >  1:1010455/1‑101 (MQ=255)
                                                                   ggTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGc                                     >  2:2070771/1‑101 (MQ=255)
                                                                   ggTA.AAAATCCGCCGGGGGGGGGGGGGGGTGTGAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGc                                     <  2:120015/101‑1 (MQ=255)
                                                                    gTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGcc                                    >  1:1671271/1‑101 (MQ=255)
                                                                     tG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAT..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCg                                   <  1:520563/101‑1 (MQ=255)
                                                                      g.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTgagatc                                  >  1:980922/1‑96 (MQ=255)
                                                                        aaaaTACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCggg                                 >  2:2957874/1‑101 (MQ=255)
                                                                          aaaacggcgcgggggggggggggtgTGGAA..CCCCCCGGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCggg                                 <  2:1621995/92‑1 (MQ=255)
                                                                           aTACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGcc                              >  1:1847047/1‑101 (MQ=255)
                                                                                 cACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggc                        >  2:1812051/1‑101 (MQ=255)
                                                                                  gggggggggggggggggTGAACCCCCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTggg                         <  1:2731375/99‑1 (MQ=255)
                                                                                  aCGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggcc                       >  1:3880567/1‑101 (MQ=255)
                                                                                   cgggggtgcgggggtgTGAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggccg                      <  1:1324956/101‑1 (MQ=255)
                                                                                    gggggtgggggggtgGGAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggccgg                     <  1:2795041/101‑1 (MQ=255)
                                                                                    gggggtgcgggggtgTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggccgg                     >  1:532269/1‑101 (MQ=255)
                                                                                          gggggggtgTGAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCcgcg               <  1:846151/99‑1 (MQ=255)
                                                                                          gcgggggtgTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCcgcg               >  1:4215751/1‑101 (MQ=255)
                                                                                           cgggggtgTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCcgcgc              >  1:3271842/1‑100 (MQ=255)
                                                                                            gggggagTGAAA..CCCCC..GCTTCTGCCGCGGCAATCTGGTCGAAAGGGTCGGGGTTTTCCCCGCCCTTTTTCTTTGGCCGGGGCCTCGGGCTGCGCCGcaac             >  2:601902/1‑101 (MQ=255)
                                                                                              gggtgTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGcgacga           >  2:319314/1‑101 (MQ=255)
                                                                                               gttgtaAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACgag          <  1:1077685/99‑1 (MQ=255)
                                                                                               gttgtaAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACgag          <  2:2348085/99‑1 (MQ=255)
                                                                                               ggtgTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACgag          >  2:3850202/1‑101 (MQ=255)
                                                                                               ggtgTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACgag          >  2:3509376/1‑101 (MQ=255)
                                                                                               ggggggAAC..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACgag          <  1:317716/96‑1 (MQ=255)
                                                                                                gtgtGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACgagg         <  2:1276012/101‑1 (MQ=255)
                                                                                                 tgtgAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACgagga        >  1:3294553/1‑101 (MQ=255)
                                                                                                  gtgAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACgaggag       >  1:3223058/1‑101 (MQ=255)
                                                                                                     aaa..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGAAc    >  2:2873876/1‑101 (MQ=255)
                                                                                                      aa..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGAACg   >  2:1925930/1‑101 (MQ=255)
                                                                                                       a..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGAACgg  >  1:3216066/1‑101 (MQ=255)
                                                                                                        |                                                                                                       
AGGTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAA.GCTGCGTGTTTTTTCTT.GGTG.AAAATACGCCACGGGGGTGCGGGGGTGTGAAA..CCCCC..GCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGAACGG  >  NC_014034/1820293‑1820493

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG