Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,732,696	+C	89.1%	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	89.1%	13.5 ‑4.3 ‑1.0	55	intergenic (‑120/‑283)	RCAP_rcc01601/tgt	FeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (3/46):  ref base (3/3):  total (6/49)
Fisher's exact test strand distribution p-value = 1.33e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.05e-01

CTAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.CC..GGGC..GGGCGG..G..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATG  >  NC_014034/1732596‑1732778                                                                                                         |                                                                                           cTAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGCTTCGCCGCGGC.GGGGCCGGGGG.G.a                                                                                             >  1:1265537/1‑101 (MQ=255)  tAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGG.G.A.c                                                                                           >  2:3834247/1‑101 (MQ=255)   agagTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.cc                                                                                          <  1:3684548/101‑1 (MQ=255)     agTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.cccg                                                                                        >  2:3278322/1‑99 (MQ=255)     agTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.cccg                                                                                        <  2:3055747/101‑3 (MQ=255)      gTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.Cc..cgg                                                                                     <  2:3756049/101‑1 (MQ=255)       tttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.Cc..cggg                                                                                    >  1:3437276/1‑100 (MQ=255)       tttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.Cc..cggg                                                                                    <  1:3517643/101‑2 (MQ=255)        ttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.Cc..cgggg                                                                                   >  2:2313585/1‑99 (MQ=255)        ttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.Cc..cgggc                                                                                   <  1:3919695/101‑3 (MQ=255)        ttATATAGGGTATCGGGGGGAGAAGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.cccttgc                                                                                     <  2:205545/101‑6 (MQ=255)         tatataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGAGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGATTCGCCGCGGC.GGGGCCGCGGCGG.C.Cc..gggc..g                                                                                 >  2:3195183/1‑101 (MQ=255)            ctaGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCc..gggc..gggg                                                                              <  2:1168233/100‑2 (MQ=255)            ataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCc..gggc..gggg                                                                              <  2:2131535/101‑2 (MQ=255)            ataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCc..gggc..gggg                                                                              <  2:1494699/101‑2 (MQ=255)             taGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.Cc..cggggcgggg                                                                              >  2:3103507/1‑94 (MQ=255)              aGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGGTTCGCCGCGGCGGGGGCCCGGGG.G.C.Cc..gggc..gggggg                                                                            >  1:1168327/1‑101 (MQ=255)               gggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.CC..CGGG..GGGGgg..g..g                                                                      >  2:2690075/1‑101 (MQ=255)                ggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCc..gggc..gggggctc                                                                          <  1:3070559/101‑6 (MQ=255)                ggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCc..gggc..gggggctc                                                                          <  1:265525/101‑6 (MQ=255)                ggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCc..gggc..gggggctc                                                                          <  2:4193958/101‑6 (MQ=255)                 gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCc..gggc..gggggctcc                                                                         <  2:1002518/101‑7 (MQ=255)                 gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCc..gggc..gggggctcc                                                                         <  1:1477786/101‑7 (MQ=255)                  tATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCc..gggc..gggggctccg                                                                        <  1:3251709/101‑8 (MQ=39)                  tATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCc..gggc..gggggctccg                                                                        <  1:3862549/101‑8 (MQ=39)                   aTCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCc..gggc..gggggctccgc                                                                       <  2:2286278/101‑9 (MQ=38)                   aTCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..CGGC..GGGGgg..g..gggg                                                                   >  1:2256525/1‑98 (MQ=255)                    tCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.CCCGGGGC..GGGGgg..g..ggcg                                                                   >  2:1792202/1‑98 (MQ=255)                     cGGGGGGAGAGGTAAAGGCCGTCGGTGGGGGGCGGAAAGCTGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCc..gggc..gggggctccgccc                                                                     <  2:2116995/101‑11 (MQ=37)                     cGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCc..gggc..gggggctccgccc                                                                     <  1:1632546/101‑11 (MQ=38)                       gggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGGCGCGGC.C.A.CC..CGGG..GGGGgg..g..ggggggccc                                                              >  1:2766531/1‑93 (MQ=255)                        ggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGGCGCGGC.G.A.CC..CGGG..GGGGGG..G..Tccccccccc                                                             >  1:4151882/1‑97 (MQ=255)                           gagaGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.CC..CCGG..GGGGGG..G..GGCCCCCCCcccc                                                          >  1:970843/1‑101 (MQ=255)                              agGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCAcccc                                                          <  2:4151620/101‑1 (MQ=25)                               ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCAccccc                                                         <  1:3287269/101‑1 (MQ=25)                               ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCAccccc                                                         <  2:3159941/101‑1 (MQ=25)                               ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCAccccc                                                         <  2:1609783/101‑1 (MQ=25)                               ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCAccccc                                                         <  2:1678316/101‑1 (MQ=25)                                 tAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCgg                                                       <  1:3278563/101‑1 (MQ=25)                                 tAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCgg                                                       <  2:2756478/101‑1 (MQ=25)                                 tAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGC..GGGGGG..G..CCGCCCCCCccccccccc                                                     >  1:414817/1‑97 (MQ=255)                                  aaaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGC..GGGGGG..G..GCGCCCCCCcccccccccg                                                    >  2:2528008/1‑96 (MQ=255)                                        cGTCGGTGCGGGGCGGACGGCCGGGTGCCAAGCAGGGTTTCGCCGCGGG.GGGGCCGCGGC.GCC.CC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTAttt                                                <  2:1530854/101‑1 (MQ=16)                                         gTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTg                                               <  1:1894797/101‑1 (MQ=21)                                          tCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTgc                                              <  1:1054370/101‑1 (MQ=21)                                          tCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGC..GGGGGG..GCCGCCCCCGCCCCCCCCGCggcctttg                                              >  2:1030967/1‑94 (MQ=255)                                              tGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCa                                          <  2:3485951/101‑1 (MQ=21)                                                 ggggCGGAAAGCCGCGTGCCAGGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGTGCCaaga                                       <  2:1265441/101‑1 (MQ=18)                                                 ggggCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTCTTTGCGCCaaga                                       <  1:2298605/101‑1 (MQ=21)                                                 ggggCAGAAGGCCGCGTGCCCAGAATGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGT..GGGGTGTCC..GCCCCCCCACCTCCGGCTTATTTGCGTCaaga                                       <  2:1080930/101‑1 (MQ=255)                                                  gggCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCaagaa                                      <  2:941662/101‑1 (MQ=21)                                                   gcgagAAAGCCGCGTGCCAGGCGGGGGTGCGCGGCGGG.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCaagaag                                     <  2:1469560/97‑1 (MQ=9)                                                     cGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAgaagaa                                   <  1:3872664/101‑1 (MQ=18)                                                      ggAGAGCCGCGTGCCAAGCGGGGTTTCGCCGGGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGaaa                                  <  1:2921000/101‑1 (MQ=17)                                                      ggAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGaaa                                  <  2:759145/101‑1 (MQ=18)                                                      ggAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGaaa                                  <  2:3776171/101‑1 (MQ=18)                                                       gAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAc                                 <  1:2126933/101‑1 (MQ=18)                                                       gAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAc                                 <  2:163316/101‑1 (MQ=18)                                                            ccGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCg                            <  2:1785524/101‑1 (MQ=18)                                                               cgTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGa                         <  2:3437024/101‑1 (MQ=18)                                                                gTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGAc                        <  2:2892352/101‑1 (MQ=18)                                                                 tGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGAcc                       <  1:1954841/101‑1 (MQ=18)                                                                   ccAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAg                     <  1:2690281/101‑1 (MQ=18)                                                                    cAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAgg                    <  2:2133135/101‑1 (MQ=18)                                                                     aaGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGc                   <  2:2256366/101‑1 (MQ=18)                                                                      aGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGcc                  <  1:3210292/101‑1 (MQ=18)                                                                          gggTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCTGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGcccccc              <  2:3061714/101‑1 (MQ=18)                                                                            gTTTCGCCGCGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTg            <  1:3254114/101‑1 (MQ=18)                                                                                cgccgcGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGc        <  2:2680733/101‑1 (MQ=18)                                                                                  ccgcGGC.GGGGCCGCGGC.GTG.CC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCg      <  1:2523589/101‑1 (MQ=21)                                                                                  ccgcGGC.GGGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCg      <  1:428965/101‑1 (MQ=18)                                                                                     cggc.ggGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCAt   <  2:2857051/101‑1 (MQ=18)                                                                                      ggc.ggGGCCGCGGC.G.ACCC..GGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATg  <  1:2528192/101‑1 (MQ=18)                                                                                                         |                                                                                           CTAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGC.GGGGCCGCGGC.G.A.CC..GGGC..GGGCGG..G..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATG  >  NC_014034/1732596‑1732778 Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG