Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014035	34,742	Δ1 bp	4.6%	coding (400/1041 nt)	RCAP_rcp00049 ←	iron siderophore/cobalamin ABC transporter permease

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014035	34,742	0	T→.	4.6%	4.9 ‑3.0 0.0	87	coding (400/1041 nt)	RCAP_rcp00049	iron siderophore/cobalamin ABC transporter permease
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (45/38):  total (45/42)
Fisher's exact test strand distribution p-value = 5.03e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCGACCCCCGCGAGGATCGTGTGGTTCGGCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCATGGGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCCCGCCGAGACGCCCAGCACGAAGGGTTC  >  NC_014035/34642‑34835                                                                                                     |                                                                                              gcgACCCCCGCGAGGATCGTGTGGTTCGGCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAt                                                                                               <  1:1691281/101‑1 (MQ=32)  cgACCCCCGCGAGGATCGTGTGGTTCGGCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCATg                                                                                              >  2:4881508/1‑101 (MQ=32)   gACCCCCGCGAGGATCGTGTGGTTCGGCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCATgg                                                                                             <  1:1615939/101‑1 (MQ=32)   gACCCCCGCGAGGATCGTGTGGTTCGGCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCATgg                                                                                             <  2:1247535/101‑1 (MQ=32)             gggATCGTGTGGTTCGGCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCA.GGgacagcgacag                                                                                  <  1:6297592/100‑1 (MQ=2)                        gTTCGGCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGcc                                                                       <  2:663548/101‑1 (MQ=11)                         ttaggCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGACCACAT.GGACAGCGACAGACCGCCCGGCCg                                                                      >  2:121583/4‑101 (MQ=11)                         ttCGGCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCg                                                                      <  1:3585551/101‑1 (MQ=14)                         ttCGGCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCg                                                                      <  1:242120/101‑1 (MQ=14)                          tCGGCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGa                                                                     <  1:228975/101‑1 (MQ=17)                          tCGGCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGa                                                                     >  2:3637905/1‑101 (MQ=17)                            ggCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATg                                                                   >  2:1789850/1‑101 (MQ=17)                            ggCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATg                                                                   <  1:4733144/101‑1 (MQ=17)                            ggCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCACAT.GGACAGCGACAGCCCGCCCGGCCGATg                                                                   >  1:5334350/1‑101 (MQ=12)                             gCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGc                                                                  >  2:7079050/1‑101 (MQ=17)                              cccGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGcc                                                                 >  1:261364/1‑101 (MQ=17)                               ccGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGccc                                                                <  1:920777/101‑1 (MQ=17)                               ccGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGccc                                                                >  2:4915313/1‑101 (MQ=17)                               ccGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGccc                                                                >  1:1184498/1‑101 (MQ=17)                               ccGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGccc                                                                <  1:3083980/101‑1 (MQ=17)                                cGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCa                                                               <  2:2310022/101‑1 (MQ=17)                                cGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCa                                                               <  2:1090400/101‑1 (MQ=17)                                 gCTGTTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAg                                                              >  1:3095140/1‑101 (MQ=12)                                 gCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAg                                                              >  1:4159982/1‑101 (MQ=17)                                  cTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGc                                                             <  2:4427655/101‑1 (MQ=17)                                  cTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGc                                                             <  2:280974/101‑1 (MQ=17)                                  cTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGc                                                             >  2:2094153/1‑101 (MQ=17)                                   tGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCa                                                            >  1:4972045/1‑101 (MQ=17)                                   tGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCa                                                            >  2:4335673/1‑101 (MQ=17)                                    ggTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCAc                                                           >  1:4132530/1‑101 (MQ=17)                                    ggTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCAc                                                           >  2:2504960/1‑101 (MQ=17)                                    ggTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCAc                                                           >  2:1194941/1‑101 (MQ=17)                                      gcgcgCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAGGGCGCCCGCGAAAGCCCCCA.GGGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGa                                                         <  1:4402269/100‑1 (MQ=255)                                        gcgcCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACACCGACAGCCCGCCCGGCCGATGCCCAGCACGATg                                                       <  1:6224901/101‑1 (MQ=12)                                          gcCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGAc                                                     <  1:79321/101‑1 (MQ=17)                                          gcCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGAc                                                     >  2:2816010/1‑101 (MQ=18)                                          gcCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGACGAc                                                     >  1:3691246/1‑101 (MQ=12)                                           ccGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGAcc                                                    <  1:3907889/101‑1 (MQ=17)                                             gTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGc                                                  >  1:5180379/1‑101 (MQ=17)                                             gTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGc                                                  >  1:1390115/1‑101 (MQ=17)                                             gTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGc                                                  >  1:3361665/1‑101 (MQ=18)                                             gTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGc                                                  >  1:697875/1‑101 (MQ=17)                                              ttCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGCCACCGACAGCCCGCCCGGCCGATGCCCAGCACGATTACCGcc                                                 >  1:1224827/1‑101 (MQ=11)                                              ttCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGcc                                                 <  2:2566739/101‑1 (MQ=17)                                                cGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCAc                                               >  2:2446535/1‑101 (MQ=18)                                                cGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCAc                                               <  2:6567515/101‑1 (MQ=17)                                                  aCAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACcg                                             >  1:924965/1‑101 (MQ=18)                                                     aaaGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCgcgc                                          >  2:5495490/1‑101 (MQ=18)                                                      aaGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGcc                                         <  2:6185409/101‑1 (MQ=17)                                                      aaGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGcc                                         <  1:2181386/101‑1 (MQ=17)                                                       aGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCg                                        <  1:4881752/101‑1 (MQ=17)                                                       aGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGAAAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCg                                        >  1:6432449/1‑101 (MQ=255)                                                           ccaccaGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCg                                    >  2:5227702/1‑101 (MQ=18)                                                            caccaGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGTCCAGCACGATGACCGCCACCGCGCCGGTCGa                                   >  2:6826595/1‑101 (MQ=255)                                                            caccaGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGa                                   >  1:1419332/1‑101 (MQ=18)                                                            caccaGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGa                                   >  2:754404/1‑101 (MQ=18)                                                              ccaGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAgg                                 >  2:5628142/1‑101 (MQ=18)                                                              ccaGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAgg                                 >  2:6083689/1‑101 (MQ=18)                                                               caGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGc                                >  1:6365985/1‑101 (MQ=17)                                                               caGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGc                                >  1:4915636/1‑101 (MQ=17)                                                               caGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGACACCGCGCCGGTCGAGGc                                >  2:358882/1‑101 (MQ=12)                                                                aGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGGCAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGcc                               <  1:1922313/101‑1 (MQ=12)                                                                aGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGcc                               <  1:2720328/101‑1 (MQ=17)                                                                 gcgcGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGccc                              <  2:1112794/101‑1 (MQ=17)                                                                  cgcgAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGcccc                             <  1:3638054/101‑1 (MQ=17)                                                                  cgcgAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGcccc                             <  1:2810459/101‑1 (MQ=17)                                                                  cgcgAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGcccc                             <  1:5363002/101‑1 (MQ=17)                                                                  cgcgAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGcccc                             <  1:2500997/101‑1 (MQ=17)                                                                   gcgAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGccccc                            <  2:261360/101‑1 (MQ=17)                                                                   gcgAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGccccc                            <  2:4567920/101‑1 (MQ=17)                                                                    cgAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGcccccc                           >  2:7015196/1‑100 (MQ=12)                                                                    cgAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCccg                           <  1:7079424/101‑1 (MQ=17)                                                                    cgAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCccg                           <  2:6039016/101‑1 (MQ=17)                                                                    cgAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGACGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCccg                           <  1:1240573/101‑1 (MQ=12)                                                                     gAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.TGACAGCGACAGCCCGCCCGGCCGATTCCCTGCACGATGACTTCCACCGCGCCGGTCCAGGACCccgc                          >  2:4348906/1‑101 (MQ=11)                                                                     gAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCccgc                          >  1:2945352/1‑101 (MQ=17)                                                                     gAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCccgc                          <  2:6540673/101‑1 (MQ=17)                                                                     gAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCccgc                          >  1:6507666/1‑101 (MQ=17)                                                                     gAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCccgc                          >  2:3512301/1‑101 (MQ=17)                                                                      aaGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCccgcc                         >  2:4163139/1‑101 (MQ=17)                                                                      aaGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCccgcc                         <  2:934033/101‑1 (MQ=17)                                                                       aGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCccgccg                        <  1:3191301/101‑1 (MQ=17)                                                                       aGGCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCccgccg                        <  1:5883223/101‑1 (MQ=17)                                                                         gCAGCGAAGGCGCCGGCG.AAGCCCCCCAGGGACGGCGACGGCCCGCCCGGCCGATGCCCAGCACGATGCCCGCCACCGCGCCGGTCGAGGCCCCCGCCgag                      <  2:6082470/101‑1 (MQ=11)                                                                         gCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCCCGCCgag                      >  2:5164778/1‑101 (MQ=17)                                                                         gCAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCCCGCCgag                      >  1:1331367/1‑101 (MQ=17)                                                                          cAGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCTACCGCGCCGGTCGAGGCCCCCGCCgaga                     <  1:6282496/101‑1 (MQ=9)                                                                           aGCGAAGGCGCCCGCGAAAGCCCCCAT.GGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCCCGCCGAGAc                    >  2:6441929/1‑101 (MQ=11)                                                                                            aaGCCCCGA.GGGACAGCGACCGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCCCGCCGAGACGCCCAGCACGAAGGGtt   <  2:5603031/101‑1 (MQ=255)                                                                                            aaGCCCCCA.GGGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCCCGCCGAGACGCCCAGCACGAAGGGtt   <  2:7070992/101‑1 (MQ=17)                                                                                             aGCCCCCA.GGGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCCCGCCGAGACGCCCAGCACGAAGGGTTc  <  2:2945243/101‑1 (MQ=17)                                                                                                     |                                                                                              GCGACCCCCGCGAGGATCGTGTGGTTCGGCCCGCTGGTCGCGCCGTTCGACAAAAGCGCCACCAGCGCGAAGGCAGCGAAGGCGCCCGCGAAAGCCCCCATGGGACAGCGACAGCCCGCCCGGCCGATGCCCAGCACGATGACCGCCACCGCGCCGGTCGAGGCCCCCGCCGAGACGCCCAGCACGAAGGGTTC  >  NC_014035/34642‑34835 Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG