Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,640,204 A→G 13.0% intergenic (‑180/+96) thrS ← / ← RCAP_rcc02462 threonyl‑tRNA synthetase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,640,2040A→G13.0% 133.6 ‑3.6 ‑0.1418intergenic (‑180/+96)thrS/RCAP_rcc02462threonyl‑tRNA synthetase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (33/21):  ref base (160/199):  total (194/224)
Fisher's exact test strand distribution p-value = 2.79e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.82e-01

Only 100 of 478 total aligned reads displayed.

CGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GG.C.TT.CGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTCCGGCTCAGA  >  NC_014034/2640105‑2640304
                                                                                                         |                                                                                                        
cGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAA.GGGGGGGGCTGCCCCCCCC.CT.GGCC.CT.CGGCCCATTCCCCCCCCCAGGGATTTTGAGCCAGAAGaaaa                                                                                                          >  1:5175079/1‑101 (MQ=255)
      caCCGTGGCGCTTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATg                                                                                                   <  2:4783237/100‑1 (MQ=21)
       aCCGTGGCGCTTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATg..g                                                                                                <  2:8519836/101‑1 (MQ=25)
          gTGGCGCTTTTGCGCAGGCAAAAGGGGGGGGGCGCCCCCCCCC.TT.GC.C.TTCCGGCCAA.TCACCCCCCCCGGCAGTTTGAGCGAGATCAAAA.ACCATG..GG..tt                                                                                           >  1:7501151/1‑101 (MQ=255)
          gTGGCGCTTGTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..tt                                                                                           <  1:439316/101‑1 (MQ=21)
           tGGCGCTTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGG.C.TT.CGGCCAATTCACCCCCCGAGGGATTTGGAGCAAGATG.AAAGACCCTG..GG..GTc                                                                                          >  1:2771956/1‑101 (MQ=21)
            ggCGCTTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCa                                                                                         <  1:8619762/101‑1 (MQ=25)
            gccgcTTTGCCGCCGGCAAGG.AGGGGCCGCGCCGCCCCTC.TT.GGCC.TT.CGCCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCa                                                                                         <  1:4243960/99‑1 (MQ=16)
             gcgcTTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCAt                                                                                        <  2:789557/101‑1 (MQ=25)
             gcgcTTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCAt                                                                                        <  1:6736064/101‑1 (MQ=25)
                cTTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCtt                                                                                     <  2:7009971/101‑1 (MQ=25)
                cTTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCtt                                                                                     <  1:4678157/101‑1 (MQ=25)
                 ttttGCGCAGGCAAAC.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCACG..GG..TTCATCTTg                                                                                    <  2:6693678/101‑1 (MQ=21)
                 ttttGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCCCTTT.GC.C.TT.TGGCCAATTCACCCCCCCACGCATTTTGGGCAAAATG.AAAGAACATG..GG..GTCATCttc                                                                                    >  2:3776823/1‑100 (MQ=12)
                    tGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCt                                                                                 <  1:6159898/101‑1 (MQ=25)
                    tGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGGATTTGGAGCCAGATG.AAAGACCAAG..GG..TTCATCTTGCCt                                                                                 >  2:2577157/1‑101 (MQ=21)
                      cgcAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCa                                                                               <  2:6146053/101‑1 (MQ=21)
                        cAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCaaa                                                                             >  2:1252188/1‑101 (MQ=21)
                          ggCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCGCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAata                                                                           <  2:2811729/101‑1 (MQ=21)
                          cccgaaa.AGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GG..TTTATCTTGCCCCAAatat                                                                          >  2:2694194/5‑101 (MQ=17)
                           gCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAatat                                                                          >  1:7100347/1‑101 (MQ=21)
                           ccgaaa.AGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GG..TTTATCTTGCCCCAAATATc                                                                         >  1:3090450/4‑101 (MQ=17)
                            cgaaa.AGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GG..TTTATCTTGCCCCAAATATcc                                                                        >  2:4542344/3‑101 (MQ=17)
                             gaaa.AGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGAC..GG..GG..TTTATCTTGCCCCAAATATccc                                                                       >  1:3313207/2‑101 (MQ=17)
                              aaa.AGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GG..TTTATCTTGCCCCAAATATCCCg                                                                      >  1:7898947/1‑101 (MQ=17)
                              aaa.AGGGGGCGCTGCCCCCCTC.TT.GG.CTTT.CGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GG..TTTATCTTGCCCCAAATATCCCg                                                                      <  2:524788/101‑1 (MQ=18)
                               aa.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCg                                                                      <  2:7566000/101‑1 (MQ=21)
                                a.GGGGGGCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCgg                                                                     <  2:5403377/101‑1 (MQ=25)
                                  ggggggCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCggg                                                                    >  1:8418753/1‑101 (MQ=21)
                                    ggggCGCTG.CCCCCTC.TT.GGCC.TT.TGGCCAATTC.ACCCCTGGGATATTTCGGGCAAGATGAA...ATCATGCCCC..TTCATCTTGCCGCAAATATCCCGGGGGGa                                                                <  2:3525400/101‑1 (MQ=2)
                                    ggggCGCTG.CCCCCTC.TT.GGCC.TT.TGGCCAATTC.ACCCCTGGGATATTTCGGGCAAGATGAA...ATCATG..CCCCTTCATCTTGCCGCAAATATCCCGGGGGGa                                                                <  2:343243/101‑1 (MQ=2)
                                     gggCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GG..TTTATCTTGCCCCAAATATCCCgggggtg                                                                <  1:969375/101‑3 (MQ=18)
                                     gggCGCTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GG..TTTATCTTGCCCCAAATATCCCgggggtg                                                                <  2:8851178/101‑3 (MQ=18)
                                       gcgcTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GG..TTTATCTTGCCCCAAATATCCCgggggtgag                                                              <  1:3905923/101‑5 (MQ=21)
                                       gcgcTGCCCCCCTC.GT.GGCC.TT.CGGCCAATTCGCCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GG..TTTATCTTGCCCCAAATATCCCgggggtgag                                                              <  1:6932700/101‑5 (MQ=18)
                                         gcTGCCCCCCTC.TT.GGCC.TT.CGGGCAATTCACCCCCCGGGGTATTTTGAGCAAGAAG.AAAGAACATG..GG..GTCATCTTGCCTCCAATATCCCCGGGGGGGTc                                                             >  2:1721476/1‑101 (MQ=12)
                                         gcTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGATGA....GG..TTTATCTTGCCCCAAATATCCCgggggtgaggc                                                            <  2:3269523/101‑7 (MQ=21)
                                         gcTGCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTc                                                             >  2:4100169/1‑101 (MQ=25)
                                            gCCCCCCTC.TT.GGCC.TT.CTGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTg                                                          <  1:6339846/101‑1 (MQ=25)
                                            gCCCCCCTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTg                                                          <  1:3640360/101‑1 (MQ=25)
                                             ccccccTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGa                                                         <  2:7368722/101‑1 (MQ=37)
                                              cccccAC.CT.GGCC.TT.CGGCCAATTCACCCCACGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGaa                                                        <  1:2927739/101‑1 (MQ=25)
                                               ccccTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGaaa                                                       >  1:1422510/1‑101 (MQ=37)
                                                cccTC.TT.GGCC.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                                      >  2:7141025/1‑101 (MQ=37)
                                                       ttggcc.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgg                                                 >  2:4649357/6‑101 (MQ=37)
                                                        tggcc.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                                                >  1:8754458/5‑101 (MQ=37)
                                                         ggcc.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggg                                               >  2:5649680/4‑101 (MQ=37)
                                                         ggcc.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggg                                               >  2:4981206/4‑101 (MQ=37)
                                                          gcc.TT.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggggg                                              >  1:2032158/3‑101 (MQ=37)
                                                               t.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCag                                          <  1:5809196/101‑1 (MQ=255)
                                                               t.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCag                                          >  1:7725284/1‑101 (MQ=37)
                                                               t.CGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCag                                          <  1:1764971/101‑1 (MQ=255)
                                                                  ggCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCagag                                        >  1:5033319/1‑101 (MQ=37)
                                                                  ggCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCagag                                        >  1:4831850/1‑101 (MQ=37)
                                                                      aaTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccc                                    <  2:1111937/101‑1 (MQ=255)
                                                                      aaTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGAACATG..GG..TTTATCTTGCCTCAAATATCTCGGGGGGAAACAGGAAAGGAGGGGGGGCACGGccca                                    >  2:2931737/1‑100 (MQ=16)
                                                                       aTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccccc                                   >  1:1093358/1‑101 (MQ=37)
                                                                       aTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccccc                                   >  1:5471368/1‑101 (MQ=37)
                                                                        ttCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTCCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccccc                                  <  1:5409568/101‑1 (MQ=37)
                                                                        ttCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccccc                                  >  2:9118210/1‑101 (MQ=25)
                                                                        ttCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccccc                                  <  2:3188428/101‑1 (MQ=255)
                                                                         tcaccccccGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCt                                 >  1:1196910/4‑101 (MQ=25)
                                                                         tcaccccccGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCt                                 <  1:5725110/98‑1 (MQ=255)
                                                                         tcaccccccGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCt                                 >  2:2398873/4‑101 (MQ=25)
                                                                          caccccccGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc                                >  1:7289234/3‑101 (MQ=255)
                                                                          caccccccGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc                                >  1:2242472/3‑101 (MQ=255)
                                                                           cccccccGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTcc                               <  1:2820490/100‑1 (MQ=255)
                                                                            ccccccGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTccc                              <  1:3037537/101‑1 (MQ=255)
                                                                            ccccccGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTccc                              <  2:8331401/101‑1 (MQ=255)
                                                                                ccGAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGt                          <  2:7947076/101‑1 (MQ=255)
                                                                                  gAGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCt                        <  2:1423796/101‑1 (MQ=255)
                                                                                   aGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCtt                       <  2:4343866/101‑1 (MQ=255)
                                                                                   aGGTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCtt                       >  2:5114242/1‑101 (MQ=255)
                                                                                    ggTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCttt                      >  2:8158060/1‑101 (MQ=255)
                                                                                    ggTATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCttt                      <  1:8258630/101‑1 (MQ=255)
                                                                                      tATTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtct                    >  2:3881474/1‑101 (MQ=255)
                                                                                       aTTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtctc                   <  2:8754053/101‑1 (MQ=255)
                                                                                       aTTTGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtctc                   <  2:2338201/101‑1 (MQ=255)
                                                                                         ttGGAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCaa                 >  1:8001882/1‑101 (MQ=255)
                                                                                         tctgagcAAGATG.AAAGCCCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCaa                 <  2:6513311/95‑1 (MQ=25)
                                                                                           ggAGCAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGc               <  2:5464834/101‑1 (MQ=255)
                                                                                            gagcAAGATG.AACGACCATG..GG..TTCATCATGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCaaaga              >  1:6513571/4‑98 (MQ=25)
                                                                                            gagcAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCg              >  1:3727543/4‑101 (MQ=255)
                                                                                              gcAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCACGCGtt            <  2:6235407/100‑1 (MQ=255)
                                                                                              gcAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGtt            <  2:2080606/100‑1 (MQ=255)
                                                                                               cAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTc           >  2:6319075/1‑101 (MQ=255)
                                                                                               cAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTc           >  2:6815769/1‑101 (MQ=255)
                                                                                               cAAGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTc           >  1:6572159/1‑101 (MQ=255)
                                                                                                aaGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTcc          <  1:9050118/101‑1 (MQ=255)
                                                                                                aaGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTcc          >  1:1341909/1‑101 (MQ=255)
                                                                                                aaGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTcc          <  1:4584204/101‑1 (MQ=255)
                                                                                                aaGATG.AAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTcc          >  1:2570351/1‑101 (MQ=255)
                                                                                                  agagtaaaGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTCCg         <  2:4242161/96‑1 (MQ=255)
                                                                                                   gatgaaaGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTCCgg        <  2:3105355/97‑1 (MQ=255)
                                                                                                   gatgaaaGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTCCgg        >  2:5848926/5‑101 (MQ=255)
                                                                                                    atgaaaGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTCCGGc       >  1:3433116/4‑101 (MQ=255)
                                                                                                       aaaGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTCCGGCTCa    >  1:523221/1‑101 (MQ=255)
                                                                                                       aaaGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTCCGGCTCa    <  1:4132570/101‑1 (MQ=255)
                                                                                                        aaGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTCCGGCTCAg   >  2:7749470/1‑101 (MQ=255)
                                                                                                         aGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTCCGGCTCAGa  >  2:190916/1‑101 (MQ=255)
                                                                                                         |                                                                                                        
CGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TT.GG.C.TT.CGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATG..GG..TTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTCCGGCTCAGA  >  NC_014034/2640105‑2640304

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG