BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,732,696	+C	100%	intergenic (‑120/‑283)	RCAP_rcc01601 ← / → tgt	FeS assembly SUF system protein/queuine tRNA‑ribosyltransferase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,732,696	1	.→C	100.0%	26.3	84	intergenic (‑120/‑283)	RCAP_rcc01601/tgt	FeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand): new base (5/69): ref base (6/4): total (11/73)
Fisher's exact test strand distribution p-value = 1.90e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.62e-03

Only 100 of 102 total aligned reads displayed.


TAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..CC..GGGC..GGGCGGGGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGC  >  NC_014034/1732597‑1732779


tAGAGTTTATATAGGGTATGGGGGGGAGAGGTAAAGGCGGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..c                                                                                        <  1:1751716/101‑1 (MQ=255)


tAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..c                                                                                        >  1:1662586/1‑101 (MQ=255)


tAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..c                                                                                        <  1:2367714/101‑1 (MQ=255)


 agagTTTATATAGGGGAGGGGGGGGAGAGGTAAAGGCCGTGGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..cc                                                                                       <  1:164070/101‑1 (MQ=255)


  gagTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCGGCGGG.G.A..ccc                                                                                      >  2:2384193/1‑100 (MQ=255)


  gagTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGG.G.A..ccc                                                                                      >  1:3796336/1‑100 (MQ=255)


  gagTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..ccc                                                                                      <  2:4096140/101‑2 (MQ=255)


   agTTTATATAGGGGATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..cccg                                                                                     <  2:1241113/101‑3 (MQ=255)


     tttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..Cc..cggg                                                                                 >  2:277566/1‑100 (MQ=255)


      ttATATAGGGTATCGGGGGGAGAGGTAAAGGCCTTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..Cc..cgggc                                                                                <  1:4936542/101‑3 (MQ=255)


      ttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..Cc..cgggc                                                                                <  2:2679004/101‑3 (MQ=255)


      ttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..Cc..cgggc                                                                                <  2:4024364/101‑3 (MQ=255)


       tatataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..g                                                                              >  2:5276873/1‑101 (MQ=255)


        atataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..GCGCGGCGGGGGCGCGGCGG.C..Cc..gggc..gg                                                                             >  1:193102/1‑101 (MQ=255)


        atataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..gg                                                                             <  2:1112151/101‑1 (MQ=255)


        atataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..gg                                                                             <  1:3299226/101‑1 (MQ=255)


        atataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..gg                                                                             <  2:4294711/101‑1 (MQ=255)


         tataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..ggg                                                                            >  1:1770843/1‑101 (MQ=255)


         aatggggTAT.GGGGGGGGAGGGTAAGGGCGGCGGTTGCGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..gggg                                                                           <  2:2988330/97‑2 (MQ=255)


          ataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..gggg                                                                           <  1:2148889/101‑2 (MQ=255)


          ataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..gggg                                                                           >  1:3448412/1‑100 (MQ=255)


          ataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..Cc..cggggcggg                                                                            >  1:3450859/1‑95 (MQ=255)


           taGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..ggggg                                                                          <  2:375780/101‑3 (MQ=255)


           taGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..ggggg                                                                          <  2:5351600/101‑3 (MQ=255)


           taGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..Cc..cggg..gggggg                                                                         >  1:2256537/1‑101 (MQ=255)


            aGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..gggggc                                                                         <  2:4536558/101‑4 (MQ=255)


             gggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..gggggct                                                                        <  2:3454115/101‑5 (MQ=255)


              ggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..gggggctc                                                                       <  1:3055720/101‑6 (MQ=255)


              ggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..CC..CGGG..GGGGggggg                                                                      >  1:3589708/1‑100 (MQ=255)


               gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..gggggctcc                                                                      <  2:2042548/101‑7 (MQ=255)


               gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..gggggctcc                                                                      <  2:3448207/101‑7 (MQ=255)


               gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGC..CGG.GGGGcc                                                                     >  1:2756947/1‑101 (MQ=255)


               gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..CC..CGGG..GGGGgggggg                                                                     >  1:2411964/1‑99 (MQ=255)


                tATCGGGGGGAGAGGTAAGGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..gggggctccg                                                                     <  2:193082/101‑8 (MQ=38)


                 aTCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.Cc..gggc..gggggctccgc                                                                    <  2:2624089/101‑9 (MQ=38)


                   cGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..CC..CGG...GGCCggggggcggcc                                                                >  2:5137377/1‑94 (MQ=255)


                    ggggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.ACCCC..GGGC..GGGGgggggcgcc                                                                  >  2:2160145/1‑96 (MQ=255)


                    ggggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGC..GGGGGGGGcccccc                                                                 >  2:1932832/1‑100 (MQ=255)


                        ggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCCCGGGGGGGGGGCGGGGG.G.C..CC..GGGG..GGGCgggggggggcgcc                                                              >  1:4914475/1‑92 (MQ=255)


                           gagGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..CC..CGGG..GGGGGGGGCCCCCCCCcccccc                                                         >  2:3472166/1‑100 (MQ=255)


                            agGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGGCGCGGC.G.A..CCCGGGGG..GGGGGGGGCCCCCCCCccccc                                                          >  1:5336502/1‑101 (MQ=255)


                             ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCAccccc                                                          <  2:5540195/101‑1 (MQ=25)


                             ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCAccccc                                                          <  2:5346912/101‑1 (MQ=25)


                             ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCAccccc                                                          <  1:1622517/101‑1 (MQ=25)


                             ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCAccccc                                                          <  2:3796094/101‑1 (MQ=25)


                             ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCAccccc                                                          <  1:5029265/101‑1 (MQ=25)


                              gTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCg                                                         <  1:848036/101‑1 (MQ=25)


                              gTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCg                                                         <  2:4850876/101‑1 (MQ=25)


                              gTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCg                                                         <  2:5143971/101‑1 (MQ=25)


                               tAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAGGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCgg                                                        <  1:4942276/101‑1 (MQ=25)


                               tAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCgg                                                        <  1:4552891/101‑1 (MQ=25)


                                aaaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGc                                                       <  2:3406267/101‑1 (MQ=25)


                                  aGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..CC..CGGG..GGGGGGGCCCCGCCCCccccccccccccg                                                  >  2:3502122/1‑93 (MQ=255)


                                     ccGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtt                                                  <  2:5485969/101‑1 (MQ=25)


                                     ccGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtt                                                  <  2:2610916/101‑1 (MQ=25)


                                      cGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAttt                                                 <  1:5188937/101‑1 (MQ=25)


                                      cGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAttt                                                 <  1:1884046/101‑1 (MQ=25)


                                       gTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTg                                                <  1:824374/101‑1 (MQ=21)


                                          ggTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTgcgc                                             <  2:399398/101‑1 (MQ=21)


                                           gTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGcc                                            <  2:1744332/101‑1 (MQ=21)


                                            tGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCa                                           <  2:5349549/101‑1 (MQ=21)


                                              cGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCaag                                         <  2:4421544/101‑1 (MQ=21)


                                                 ggCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCaagaag                                      <  2:5262807/101‑1 (MQ=18)


                                                  gCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAagaaga                                     <  1:1155130/101‑1 (MQ=18)


                                                    ggAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGaaa                                   <  1:5302905/101‑1 (MQ=18)


                                                    ggAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGaaa                                   <  1:84661/101‑1 (MQ=18)


                                                    ggAAAGCCGCGTGCCAAGCAGGCTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGaaa                                   <  1:4141648/101‑1 (MQ=18)


                                                     gAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAc                                  <  2:5347890/101‑1 (MQ=18)


                                                      aaaGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACa                                 <  2:4423281/101‑1 (MQ=18)


                                                       aaGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACaa                                <  1:2213836/101‑1 (MQ=18)


                                                        aGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAg                               <  2:2199034/101‑1 (MQ=18)


                                                        aGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAg                               <  2:3039614/101‑1 (MQ=18)


                                                         gCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGc                              <  2:3410919/101‑1 (MQ=18)


                                                           cgcgTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCgg                            <  1:978252/101‑1 (MQ=18)


                                                           cgcgTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCgg                            <  1:5137663/101‑1 (MQ=18)


                                                            gcgTGCGAAGCAGGGTGTCG..CCGCGGCGGGGCCCCGGCGG.C..CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCggg                           <  1:3917827/101‑1 (MQ=17)


                                                            gcgTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCggg                           <  1:1541876/101‑1 (MQ=18)


                                                             cgTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGa                          <  2:3589491/101‑1 (MQ=18)


                                                             cgTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGa                          <  2:391310/101‑1 (MQ=18)


                                                             cgTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGa                          <  2:4109848/101‑1 (MQ=18)


                                                              gTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGAc                         <  2:2089767/101‑1 (MQ=18)


                                                               tGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGAcc                        <  1:2018076/101‑1 (MQ=18)


                                                               tGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGAcc                        <  1:2478254/101‑1 (MQ=18)


                                                                gCCAAGCAGTGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCa                       <  1:1828841/101‑1 (MQ=18)


                                                                 cttccgatcttttCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAg                      <  2:4332799/91‑1 (MQ=17)


                                                                 ccAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAg                      <  2:2756786/101‑1 (MQ=18)


                                                                 ccAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAg                      <  2:3570389/101‑1 (MQ=18)


                                                                   aaGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGc                    <  1:277599/101‑1 (MQ=18)


                                                                     gCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGccc                  <  2:4184079/101‑1 (MQ=18)


                                                                     gCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGccc                  <  2:1867666/101‑1 (MQ=18)


                                                                         ggTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCt              <  2:3223585/101‑1 (MQ=18)


                                                                         ggTTTCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCt              <  1:4235101/101‑1 (MQ=18)


                                                                          gTTTCG..CCGCGGGGGGGCCGCGGC.GCC..CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTg             <  2:1615486/101‑1 (MQ=18)


                                                                            ttCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCa           <  1:3923805/101‑1 (MQ=18)


                                                                            ttCG..CCGCGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCa           <  1:1662443/101‑1 (MQ=18)


                                                                               g..ccgcGGCGGGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGcc        <  1:24039/101‑1 (MQ=18)


                                                                                     gggcggGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCAt    <  2:3251576/100‑1 (MQ=18)


                                                                                     cggcggGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCAt    <  1:41899/101‑1 (MQ=18)


                                                                                      gggggggCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATg   <  1:3937703/98‑1 (MQ=18)


                                                                                       gcggGGCCGCGGC.G.AC.CC..GGGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGc  <  1:2160267/101‑1 (MQ=18)


TAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCG..CCGCGGCGGGGCCGCGGC.G.A..CC..GGGC..GGGCGGGGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGC  >  NC_014034/1732597‑1732779

Base quality scores:


ATCG

< 3 ≤


ATCG

< 20 ≤


ATCG

< 26 ≤


ATCG

< 33 ≤


ATCG

< 39 ≤


ATCG