BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,441,191	T→C	6.4%	L45P (CTG→CCG) ‡	gpI →	tail formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,441,191	0	T	C	6.4%	1070.1 / 39.2	343	L45P (CTG→CCG) ‡	gpI	tail formation protein
Reads supporting (aligned to +/- strand): ref base T (160/157); new base C (11/11); total (171/172)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 361 total aligned reads displayed.

CGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGG  >  NC_002947/3441042‑3441340

cGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCctggct                                                                                                                                                       <  1:702281/150‑1 (MQ=255)

   cGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCtggctggc                                                                                                                                                    <  1:214121/150‑1 (MQ=255)

    gAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGcc                                                                                                                                                   >  2:5970215/1‑150 (MQ=255)

      aGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTg                                                                                                                                                 <  2:679971/150‑1 (MQ=255)

       ggtggtGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTgg                                                                                                                                                <  2:899063/150‑1 (MQ=255)

       ggtggtGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTgg                                                                                                                                                <  2:5654045/150‑1 (MQ=255)

        gtggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTggg                                                                                                                                               <  1:4451093/150‑1 (MQ=255)

           gtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAt                                                                                                                                            >  1:5336353/1‑150 (MQ=255)

             gACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATgg                                                                                                                                          >  2:4582577/1‑150 (MQ=255)

               cATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGaa                                                                                                                                        <  1:2360420/150‑1 (MQ=255)

               cATGAGCCTGTGACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCAGGCCCAGGCCAGCCAGGGCaggaa                                                                                                                                         <  2:4904424/150‑5 (MQ=255)

                   aGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTc                                                                                                                                    <  1:2069188/150‑1 (MQ=255)

                    gCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGACCTTCCATGCCCAGGCCAGCCAGGGCA‑GGAAGTTCg                                                                                                                                   <  1:5979096/150‑1 (MQ=38)

                     ccTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGTCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGa                                                                                                                                  <  1:1525087/150‑1 (MQ=255)

                     ccTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGa                                                                                                                                  <  2:3150934/150‑1 (MQ=255)

                     ccTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGa                                                                                                                                  <  2:1134881/150‑1 (MQ=255)

                     ccTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGa                                                                                                                                  <  2:1783171/150‑1 (MQ=255)

                     ccTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCACGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGa                                                                                                                                  >  1:3566595/1‑150 (MQ=255)

                      cTGTTACCGCACACCGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTAGCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAg                                                                                                                                 <  1:218407/150‑1 (MQ=255)

                      cTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAg                                                                                                                                 <  2:2544193/150‑1 (MQ=255)

                       tGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAgg                                                                                                                                >  2:432228/1‑150 (MQ=255)

                           aCCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGACGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTg                                                                                                                            <  2:1374241/150‑1 (MQ=255)

                            ccGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACCGTCTCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCCTGGCCGGTCGAGGGCTgg                                                                                                                           >  2:6364980/1‑150 (MQ=255)

                            ccGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTgg                                                                                                                           >  1:4644229/1‑150 (MQ=255)

                             cGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCCTGGAAGGTCGAGGGCgggg                                                                                                                          >  1:82886/1‑150 (MQ=255)

                             cGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCCTGGAAGGTCGAGGGCTggg                                                                                                                          >  2:2745051/1‑150 (MQ=255)

                             cGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTggg                                                                                                                          >  1:3257553/1‑150 (MQ=255)

                               cacaACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAg                                                                                                                        >  2:3587886/1‑150 (MQ=255)

                                    cGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcggc                                                                                                                   <  2:2108567/150‑1 (MQ=255)

                                     gCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCt                                                                                                                  <  2:3895338/150‑1 (MQ=255)

                                        aCCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTac                                                                                                               <  1:6543724/150‑1 (MQ=255)

                                            aGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACcg                                                                                                           <  1:6693464/150‑1 (MQ=255)

                                            aGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACcg                                                                                                           <  2:4726504/150‑1 (MQ=255)

                                              ttGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACcgac                                                                                                         <  2:3821010/150‑1 (MQ=255)

                                                ggAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACcgacga                                                                                                       <  1:5483868/150‑1 (MQ=255)

                                                   gcgcgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCa                                                                                                    <  1:129514/150‑1 (MQ=255)

                                                   gcgcgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCa                                                                                                    >  2:5799603/1‑150 (MQ=255)

                                                       gcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTggg                                                                                                                          <  2:5444646/124‑1 (MQ=255)

                                                        cccTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGcgcc                                                                                               <  2:4760008/150‑1 (MQ=255)

                                                             gAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCcgcgcg                                                                                          >  1:6811455/1‑150 (MQ=255)

                                                                gcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGACGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTg                                                                                       >  2:3010154/1‑150 (MQ=255)

                                                                gcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGCGGCGGCTTACACCGACGACCAGCGCCGCGCGCTg                                                                                       >  2:6598674/1‑150 (MQ=255)

                                                                gcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTg                                                                                       >  1:2690194/1‑150 (MQ=255)

                                                                   gccgccGACCTGGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCGCGGCGCTGGTc                                                                                    <  2:1268016/150‑1 (MQ=255)

                                                                     cgccgACCTTGGCCTTGACCCGGACATTATCCGTGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcg                                                                                  >  1:1003204/1‑150 (MQ=255)

                                                                     cgccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcg                                                                                  >  1:5747833/1‑150 (MQ=255)

                                                                      gccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcgc                                                                                 >  1:187720/1‑150 (MQ=255)

                                                                         gACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAg                                                                              <  2:6145757/150‑1 (MQ=255)

                                                                          aCCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAgg                                                                             <  1:2612337/150‑1 (MQ=255)

                                                                             ttGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCa                                                                          <  1:5544573/150‑1 (MQ=255)

                                                                             ttGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCa                                                                          >  2:2285021/1‑150 (MQ=255)

                                                                             gtGGCCCTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACGTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCa                                                                          <  1:4859526/149‑1 (MQ=255)

                                                                              tGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCAt                                                                         >  2:5878211/1‑150 (MQ=255)

                                                                                  cTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCg                                                                     <  1:3738569/150‑1 (MQ=255)

                                                                                     gACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGtt                                                                  <  1:481484/150‑1 (MQ=255)

                                                                                       cccGGATAATGTCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                                >  2:5211749/1‑150 (MQ=255)

                                                                                        ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                               <  2:3103282/150‑1 (MQ=255)

                                                                                        ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                               >  2:4432363/1‑150 (MQ=255)

                                                                                         cGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcaca                                                              >  1:3123723/1‑150 (MQ=255)

                                                                                          ggACATTATCCGGGGCGTGGCCGACTCGTCACGCGGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACaa                                                             <  1:2943233/150‑1 (MQ=255)

                                                                                                tATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCaa                                                       <  2:3473255/150‑1 (MQ=255)

                                                                                                gaTCCGGGGCGTGGCCGACTCGGCACGCCGCCCGCCGAACTTCC‑CGCCCGGGCGGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCGTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCCGCCCGGGTCACAAGCCCaa                                                       <  2:2882868/149‑1 (MQ=255)

                                                                                                 aTCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑GGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAg                                                      <  1:4064939/150‑1 (MQ=255)

                                                                                                   ccGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAACGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAggg                                                    >  1:285493/1‑150 (MQ=255)

                                                                                                     ggggCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGa                                                                               >  2:4284965/1‑121 (MQ=255)

                                                                                                      gggCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCAc                                                 >  2:4123418/1‑150 (MQ=255)

                                                                                                         cGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGt                                              >  2:4855640/1‑150 (MQ=255)

                                                                                                             gCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgc                                          >  1:4784795/1‑150 (MQ=255)

                                                                                                              ccGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcg                                         >  1:6164949/1‑150 (MQ=255)

                                                                                                                gACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGcgcg                                       >  1:2884036/1‑150 (MQ=255)

                                                                                                                   tCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGAtt                                    >  1:1056696/1‑150 (MQ=255)

                                                                                                                   tCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGAtt                                    >  1:660571/1‑150 (MQ=255)

                                                                                                                     gTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcg                                  >  1:3813984/1‑150 (MQ=255)

                                                                                                                        aCGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcg                               >  2:1457086/1‑150 (MQ=255)

                                                                                                                              cccGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCt                         >  1:3172347/1‑150 (MQ=255)

                                                                                                                               ccgccgAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTc                        >  1:1756082/1‑150 (MQ=255)

                                                                                                                                cgccgAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCa                       >  2:3710659/1‑150 (MQ=255)

                                                                                                                                 gccgAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCaa                      >  1:346909/1‑150 (MQ=255)

                                                                                                                                 gccgAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCaa                      >  1:626200/1‑150 (MQ=255)

                                                                                                                                  ccgAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAg                     <  2:2420875/150‑1 (MQ=255)

                                                                                                                                  ccgAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAg                     <  2:1900231/150‑1 (MQ=255)

                                                                                                                                    gAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGc                   >  1:6744453/1‑150 (MQ=255)

                                                                                                                                     aaCGTCC‑GGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCg                  <  1:6364980/150‑1 (MQ=255)

                                                                                                                                       cttccatGCCCAGGCCAGCCAGGGCA‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGt                >  1:770539/7‑150 (MQ=255)

                                                                                                                                        ttccatGCCCAGGCCAGCCAGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGt                >  2:5570357/6‑150 (MQ=255)

                                                                                                                                          ccatGCCCAGGCCAGCCAGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTc              >  1:5072012/4‑150 (MQ=255)

                                                                                                                                          cc‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCg             <  2:3272432/150‑1 (MQ=255)

                                                                                                                                           c‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg            <  1:6564943/150‑1 (MQ=255)

                                                                                                                                           c‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg            <  2:314449/150‑1 (MQ=255)

                                                                                                                                            atGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg            >  1:4441187/2‑150 (MQ=255)

                                                                                                                                            atGCCCAGGCCAGCCAGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg            >  2:6523849/2‑150 (MQ=255)

                                                                                                                                             tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCggg           <  2:4956403/150‑1 (MQ=255)

                                                                                                                                             tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCggg           >  2:1041751/1‑150 (MQ=255)

                                                                                                                                              gCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGt          >  1:2430845/1‑150 (MQ=255)

                                                                                                                                               cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTc         >  1:5699999/1‑150 (MQ=255)

                                                                                                                                                ccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCa        >  2:5136259/1‑150 (MQ=255)

                                                                                                                                                    gctggcCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAAcgc    <  1:1807499/150‑1 (MQ=255)

                                                                                                                                                     ctggcCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAAcgcg   >  2:4008607/1‑150 (MQ=255)

                                                                                                                                                      tggcCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCgg  >  2:5068515/1‑150 (MQ=255)

                                                                                                                                                      tggcCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCgg  >  2:282761/1‑150 (MQ=255)

CGTCGAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGG  >  NC_002947/3441042‑3441340

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