BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	939,928:1	+C	7.8%	coding (5697/18933 nt)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	939,928	1	.	C	7.8%	782.9 / 82.4	255	coding (5697/18933 nt)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand): ref base . (100/135); new base C (10/10); total (110/145)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 257 total aligned reads displayed.

GAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCCAG  >  NC_002947/939782‑940077

gAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑Gtg                                                                                                                                                     <  2:5613803/150‑1 (MQ=255)

gAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑Gtg                                                                                                                                                     <  2:4998077/150‑1 (MQ=255)

         cctgccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGGCAG‑CAC‑CGTACGCGTGCTGGCCgcc                                                                                                                                            <  2:1083521/150‑3 (MQ=14)

          ctgccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGtt                                                                                                                                           <  1:1394479/150‑1 (MQ=255)

           tgccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTc                                                                                                                                          >  2:2669827/1‑150 (MQ=255)

            gccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCg                                                                                                                                         <  1:1588268/150‑1 (MQ=255)

            gccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCg                                                                                                                                         <  2:4695400/150‑1 (MQ=255)

                cATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCCGTTCGgcg‑c                                                                                                                                    >  1:3905415/1‑150 (MQ=255)

                 aTCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGcg‑cg                                                                                                                                   <  2:1276422/150‑1 (MQ=255)

                      tGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAC‑CACCCCGAAC‑GTGCTGCCCCGTTCGGCG‑CGGcccg                                                                                                                              >  1:5428299/1‑148 (MQ=14)

                       ggCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTCCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGGACG‑GTGCTGCCCGGTTCGGCC‑CGGCCGGt                                                                                                                             >  2:5524095/1‑150 (MQ=255)

                         ccGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGa                                                                                                                           <  2:3139088/150‑1 (MQ=255)

                         ccGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAC‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGGGa                                                                                                                           >  1:1044409/1‑150 (MQ=255)

                          cGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCG‑CACGCCTACG‑GTGCTGCCCCGTTCGGCG‑CGGCCGGTgtc                                                                                                                          >  1:5509237/1‑148 (MQ=255)

                          cGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGAc                                                                                                                          >  2:1743733/1‑150 (MQ=255)

                            tGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACCCGTACG‑GTGCTGCCCCGTTCGGCG‑CCGCCCGGGAccc                                                                                                                        >  1:4842782/1‑149 (MQ=14)

                                ttGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCt                                                                                                                    <  2:1958602/150‑1 (MQ=255)

                                   aaGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCAGcc                                                                                                                 >  2:2048459/1‑150 (MQ=255)

                                   aaGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGcc                                                                                                                 <  1:2587535/150‑1 (MQ=40)

                                   aaGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGcc                                                                                                                 >  2:6533701/1‑150 (MQ=40)

                                   aaGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCGACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGcc                                                                                                                 <  2:5735648/150‑1 (MQ=255)

                                        ccTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCgg                                                                                                            <  2:3870411/150‑1 (MQ=38)

                                         cTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGc                                                                                                           <  2:2941019/150‑1 (MQ=38)

                                         cTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAC‑CACCCGGACG‑GTGCTGCCCCGGTCGGCG‑CGGCCGGTGACCACGCGGCCggcggc                                                                                                           >  2:5957032/1‑150 (MQ=11)

                                           cGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCAt                                                                                                         <  2:113627/150‑1 (MQ=35)

                                           cGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCAt                                                                                                         <  1:2478209/150‑1 (MQ=35)

                                           cGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCAt                                                                                                         <  1:3153515/150‑1 (MQ=35)

                                            gCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCAtt                                                                                                        <  1:4375605/150‑1 (MQ=35)

                                             ccGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTg                                                                                                       <  2:2175492/150‑1 (MQ=34)

                                               gTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGAt                                                                                                     <  2:1198983/150‑1 (MQ=34)

                                                  tATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCAc                                                                                                  <  1:4917538/150‑1 (MQ=34)

                                                  tATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCAc                                                                                                  <  1:1859654/150‑1 (MQ=34)

                                                   aTCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGGGGGCAG‑CACGCCGACG‑GGGCTGCCCGGTGGGGCC‑CGGCCCGTGGCCACCCGGCCGTCGGCAGTGGTAAcc                                                                                                 >  1:1497989/1‑150 (MQ=11)

                                                    tCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCa                                                                                                <  2:3473782/150‑1 (MQ=34)

                                                           cGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTtcggt                                                                                         <  1:3035276/150‑1 (MQ=34)

                                                             ggTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCGCCACCCG‑ACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCAGTGATCAACAGgtcggtcg                                                                                       >  2:2364528/1‑150 (MQ=11)

                                                             ggTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTtcggtcg                                                                                       <  1:5950013/150‑1 (MQ=34)

                                                             ggTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTtcggtcg                                                                                       <  1:5859239/150‑1 (MQ=34)

                                                                 gAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGc                                                                                   <  2:3733002/150‑1 (MQ=34)

                                                                 gAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGc                                                                                   <  2:4516836/150‑1 (MQ=34)

                                                                    ggTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCG‑CACGCGTACG‑GTGCTGCCCCGTTTGGCG‑CGGCCGGGGACCACCCAGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGc                                                                                >  2:717378/1‑150 (MQ=11)

                                                                    ggTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGGACG‑GTGCTGCCCGGGTCGGCG‑CGGCCGGGGGCCACCCAGCCGGCGGCATTTATGACCACTTCGGTCGGTGGGGc                                                                                >  1:1401929/1‑150 (MQ=11)

                                                                      tGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGGTCGGCG‑CGGCCGGTGTCCACGCTGCCGGCGGCATTGGTAAACACTGCGGTCGGggcggccg                                                                              >  2:2470295/1‑148 (MQ=11)

                                                                         tGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCCCTTCGGTCGGTGCGGCAgggg                                                                           >  1:6819064/1‑150 (MQ=21)

                                                                            tGAAGCTGCCGGTAGCGCCATCTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCgg                                                                        <  1:3391142/150‑1 (MQ=21)

                                                                            tGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCgg                                                                        <  1:5686603/150‑1 (MQ=34)

                                                                            tGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCgg                                                                        <  2:6022956/150‑1 (MQ=34)

                                                                            tGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCAATGATCACCAGTTCGGTCGGTGCGGCAGGGGCgg                                                                        >  2:5269345/1‑150 (MQ=21)

                                                                            tGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGGTCGGCG‑CGGGCGGTGACCCCGCGGCCGGCGGCAGTGTTCACCAGTTCGGTCGGGGCGGCGGGGGCgg                                                                        >  1:1024111/1‑150 (MQ=11)

                                                                            tGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGGACG‑GGGCTGCCCCGGTCGGCG‑CGGCCGGGGGCCACGCTGCCGGCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCgg                                                                        >  1:4460213/1‑150 (MQ=11)

                                                                             gAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGt                                                                       <  1:2547845/150‑1 (MQ=34)

                                                                             gAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGt                                                                       <  1:4194361/150‑1 (MQ=34)

                                                                                  tGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCG‑CACGCGTACG‑GTGCTGGCCGGGTCGGCG‑CGGCCGGTGACCACGCGGCCGTCGGCATTTATAACCAGTTCGGGCGGTGCGGCAGGGGCGGTGatat                                                                  >  1:2783915/1‑150 (MQ=11)

                                                                                   gCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATc                                                                 >  2:1544647/1‑150 (MQ=34)

                                                                                     cGGTAGCGCCAACTCTCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcg                                                               <  2:1388918/150‑1 (MQ=21)

                                                                                       gTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGGACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGGCGGCATTGATCACCAGGTCGGTCGGGGCGGCAGGGGCGGTGATATCcggg                                                             >  2:3928991/1‑149 (MQ=255)

                                                                                         aGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcggtgc                                                           <  2:6843099/150‑1 (MQ=34)

                                                                                          gcgcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcggtgcg                                                          <  1:4138838/150‑1 (MQ=34)

                                                                                          gcgcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCCTACG‑GTGCTGCCCGGGTCGGCG‑CGGCCGGTGACCACGCAGCCGGCGGCCTTGATCACCAGTGCGGTCGGGGCGGCAGGGGCGGTGATATCcggtgcg                                                          >  1:2165874/1‑150 (MQ=11)

                                                                                           cgcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcggtgcgg                                                         <  1:5154359/150‑1 (MQ=34)

                                                                                            gcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGcc                                                                                                                                                >  2:882478/1‑63 (MQ=35)

                                                                                            gcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGcc                                                                                                                                                <  1:882478/63‑1 (MQ=35)

                                                                                            gcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcggtgcggt                                                        >  2:476662/1‑150 (MQ=34)

                                                                                             ccAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACCCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACCCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATAcggtgcggtg                                                       >  1:5642368/1‑150 (MQ=255)

                                                                                              cAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGa                                                      <  2:6456386/150‑1 (MQ=34)

                                                                                                  tATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCg                                                                                                                               <  1:3826373/73‑1 (MQ=39)

                                                                                                  tATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcggtgcggtgtgttg                                                  >  2:3384978/1‑145 (MQ=21)

                                                                                                   aTCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgg                                                 <  2:6594127/150‑1 (MQ=33)

                                                                                                   aTCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgg                                                 <  2:3614733/150‑1 (MQ=33)

                                                                                                      aCGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCG‑CCCGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGGGCCCACGCTGCCGTCGGCCGTGATCACCAGTTCGGTCGGGGCGGCAGGGGCGGTGAGATCCGGTGCGGGGATTtgggtg                                              >  2:4367010/1‑150 (MQ=11)

                                                                                                      aCGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAC‑CACCCGTACG‑GTGCTGCCCCGGTCGGCG‑CGGCCCGTTGCCACCCAGCCGTCGGCCTTGATCACCAGTTCGGTCGGTGCGGCACGGGCGGTGAAATGCGGTGCGGTGATTtgggtg                                              >  2:2528247/1‑150 (MQ=11)

                                                                                                       cGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCG‑CCCGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCACTTCGGTCGGTGCGGCAGGGGCGGGGATATACGGTGCGGTGATTtgggtgg                                             >  1:2817910/1‑150 (MQ=11)

                                                                                                       cGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtgg                                             <  2:5105364/150‑1 (MQ=32)

                                                                                                       cGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtgg                                             <  2:3308710/150‑1 (MQ=32)

                                                                                                       cGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtgg                                             <  1:84534/150‑1 (MQ=32)

                                                                                                              cAGCACAGTGGTGCCGTCGGGGGCCCC‑CACCCGGAAC‑GAGCTGCCCGGGTCGGGGCCGGCCGG‑GACCACCCAGCCGGGGGCCGTGATAACCCGTACGGTCGGGGCGGGGGGGGCGGGGATATCCGGTGCGGGGAGTTGGGTGGgcggcgg                                      >  2:6406485/1‑149 (MQ=2)

                                                                                                              cAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGa                                      <  1:996787/150‑1 (MQ=32)

                                                                                                              cAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGa                                      <  1:4671528/150‑1 (MQ=32)

                                                                                                               aGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAg                                     <  2:5177542/150‑1 (MQ=32)

                                                                                                                 cacaGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCGCGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGa                                   <  1:2903898/150‑1 (MQ=17)

                                                                                                                 cacaGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGa                                   <  2:4812079/150‑1 (MQ=32)

                                                                                                                 cacaGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGa                                   <  1:5890969/150‑1 (MQ=32)

                                                                                                                    aGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAac                                >  2:5041110/1‑150 (MQ=32)

                                                                                                                     gtggtgCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAaca                               <  1:5140548/150‑1 (MQ=32)

                                                                                                                      tggtgCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAacac                              <  1:3992390/150‑1 (MQ=32)

                                                                                                                         tgCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACAccgc                           <  2:3866424/150‑1 (MQ=32)

                                                                                                                             gTCGGCGGCCAT‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGTGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGtc                       <  2:409113/150‑1 (MQ=17)

                                                                                                                                ggcggcCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGtcatc                    >  1:6282556/1‑150 (MQ=32)

                                                                                                                                  cgggCAG‑CAC‑CGTACGCGTGCTGGCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCaa                  >  1:2097511/1‑150 (MQ=11)

                                                                                                                                   gggcAG‑CAC‑CGTACGCGTGCTGGCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAAt                 >  1:3571733/4‑150 (MQ=11)

                                                                                                                                   gggcAG‑CAC‑CGTACGCGTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGAt                                                     >  1:5013856/4‑114 (MQ=11)

                                                                                                                                    tccAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATg                <  2:2197255/149‑1 (MQ=17)

                                                                                                                                    gcCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATg                <  1:4324517/150‑1 (MQ=32)

                                                                                                                                     ccAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATgg               <  1:6542789/150‑1 (MQ=32)

                                                                                                                                      cgg‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGc              <  2:222607/148‑1 (MQ=17)

                                                                                                                                        g‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGcca            >  1:6867015/1‑150 (MQ=32)

                                                                                                                                             gcgTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCcaccac        >  1:2349115/1‑150 (MQ=32)

                                                                                                                                                 aCG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTcc    >  1:6010331/1‑150 (MQ=32)

                                                                                                                                                  cG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCca   <  1:5566733/150‑1 (MQ=32)

                                                                                                                                                   g‑gTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCcag  >  1:6455311/1‑150 (MQ=32)

GAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAG‑CACGCGTACG‑GTGCTGCCCGGTTCGGCG‑CGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCCAG  >  NC_002947/939782‑940077

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`