BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,675,147	A→C	9.0%	M166L (ATG→CTG) ‡	PP_2340 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,675,147	0	A	C	9.0%	450.1 / 30.4	181	M166L (ATG→CTG) ‡	PP_2340	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (91/66); new base C (7/10); total (98/83)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

Only 100 of 181 total aligned reads displayed.

GACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGCTG  >  NC_002947/2675000‑2675295

 aCCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATgg                                                                                                                                                   >  2:2902320/1‑150 (MQ=255)

  ccTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATggg                                                                                                                                                  <  1:585513/150‑1 (MQ=255)

  ccTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATggg                                                                                                                                                  >  2:4255073/1‑150 (MQ=255)

    tGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCa                                                                                                                                                >  2:5386990/1‑150 (MQ=255)

    tGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCa                                                                                                                                                <  1:5408208/150‑1 (MQ=255)

    tGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCa                                                                                                                                                <  1:1927643/150‑1 (MQ=255)

     ggAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAg                                                                                                                                               >  1:3685972/1‑150 (MQ=255)

       aGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAgcg                                                                                                                                             <  1:5114852/150‑1 (MQ=255)

         gCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAgcgcg                                                                                                                                           >  2:1552576/1‑150 (MQ=255)

         gCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAgcgcg                                                                                                                                           <  2:3753324/150‑1 (MQ=255)

         gCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCAGCGCGATGCCATGGGCAgcgcg                                                                                                                                           >  2:4588554/1‑150 (MQ=255)

           cTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCa                                                                                                                                         >  1:4361251/1‑150 (MQ=255)

            tATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCCGCGCGCAg                                                                                                                                        >  1:469224/1‑150 (MQ=255)

            tATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAg                                                                                                                                        >  2:263598/1‑150 (MQ=255)

            tATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAg                                                                                                                                        >  1:5441495/1‑150 (MQ=255)

            tATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAg                                                                                                                                        >  2:190274/1‑150 (MQ=255)

            tATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAg                                                                                                                                        >  1:4279796/1‑150 (MQ=255)

                     gtgtTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGa                                                                                                                               <  2:377948/150‑1 (MQ=255)

                       gtTGCTGCCGGTGTTCGGTGACTGGGAACTGCCCCTATTGCGCGATGCGATCG‑ATCGCGAGTTCTGCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCTTGGGCAGGGCGCAGGGCCAGTGAGc                                                                                                                             <  2:1380479/150‑1 (MQ=255)

                       gtTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGc                                                                                                                             >  2:3156615/1‑150 (MQ=255)

                        ttGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCCGCGCGCAGGGCCAGTGAGcc                                                                                                                            >  2:1441959/1‑150 (MQ=255)

                         tgctgcCGGTGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCTATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCg                                                                                                                           <  1:1300614/150‑1 (MQ=255)

                            tgcCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCACAGGGCCAGTGAGCCGGGc                                                                                                                        >  2:2448551/1‑150 (MQ=255)

                               cGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGaa                                                                                                                     <  2:2191366/150‑1 (MQ=255)

                                ggAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTCCCTGCGCGCTGCCCATGGCATCGCGCa                                                                                                                                         >  1:2950530/1‑129 (MQ=255)

                                 gAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAACCGGGCGAATc                                                                                                                   <  2:4692281/150‑1 (MQ=255)

                                 gAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAACCGGGCGAATc                                                                                                                   <  2:4692359/150‑1 (MQ=255)

                                      ccGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATAGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCg                                                                                                              >  2:2654307/1‑150 (MQ=255)

                                       cGTGACGGGGACCTGCCCCAATTGCGCGAGGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCCGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGc                                                                                                             <  2:1614493/150‑1 (MQ=255)

                                         tGACTGGGACCTGCCCCAATTGCGCGATGCTATCG‑ATCGCGTGTTCTTCCTGCTTTCAGAGGCGCACTCCCAGCGCTATATCGCCCCCCGCTTCCTGCGCGTTGCCTGGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCg                                                                                                           <  1:737503/150‑1 (MQ=255)

                                            cTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGgtg                                                                                                        >  2:4150731/1‑150 (MQ=255)

                                             tGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGgtgt                                                                                                       >  2:301075/1‑150 (MQ=255)

                                               ggAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGCGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCCGGCGACTGATTCGcccggctca                                                                                                     <  1:4465030/150‑8 (MQ=255)

                                               ggAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCGTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTAt                                                                                                     <  2:4263260/150‑1 (MQ=255)

                                                 aaCTGCCCCAATGGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTc                                                                                                   <  2:777262/150‑1 (MQ=255)

                                                        ccAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCAgg                                                                                                                                       >  2:3216062/1‑107 (MQ=38)

                                                           aTTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTg                                                                                         <  2:4440920/150‑1 (MQ=255)

                                                            ttGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGTTGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTgg                                                                                        <  2:537337/150‑1 (MQ=255)

                                                             ctcgcgAT‑CGATCGCATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGt                                                                                       >  2:125322/3‑150 (MQ=255)

                                                              gcgcGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTg                                                                                      >  1:3921458/1‑150 (MQ=255)

                                                                gcgatgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGTGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCa                                                                                    <  1:3601537/150‑1 (MQ=255)

                                                                gcgatgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCa                                                                                    >  2:3492200/1‑150 (MQ=255)

                                                                   atgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGt                                                                                 >  2:702538/1‑150 (MQ=255)

                                                                   atgcgatCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCACCCAGCTTCCTGCGCGATGCCATGGGCCGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGt                                                                                 >  2:2584878/1‑150 (MQ=255)

                                                                    tgcgatCG‑ATCGCGAGTTCTTCCTGCTAACAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGGGTATTCGCAACCGCTGGTGCAGGtt                                                                                >  2:1469508/1‑150 (MQ=255)

                                                                      cgatcg‑atCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTc                                                                              >  1:2400521/1‑150 (MQ=255)

                                                                       gatcg‑atcGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCCGCGCGCAGGGCCAGGGAGCCGGGCGAAGCAGTCGCCGGTGGCTTCGCAACCGCTGGTGCCGGTTTCg                                                                             >  1:4854696/1‑150 (MQ=255)

                                                                       gatcg‑atcGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCg                                                                             >  1:4442768/1‑150 (MQ=255)

                                                                        ttcg‑atcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGt                                                                            <  2:4976207/149‑1 (MQ=255)

                                                                        atct‑atcgGGAGTTCTTCCTGCTTTCAGAGGCGCACTACCCGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCTTGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGt                                                                            <  1:3470850/150‑1 (MQ=255)

                                                                        atcg‑atcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGt                                                                            <  2:2950304/150‑1 (MQ=255)

                                                                        atcg‑atcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGAAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGt                                                                            <  1:1292729/150‑1 (MQ=255)

                                                                         tcg‑atcgCGTGTTCTTCCTGCTTTCAGCGGCGCACTCCCCGCGCTATAGCGCCCCCCGCTTCCTGCGCGATGCCCTGGGCCGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCACAACCGCTGGTGCAGGTTTCGTg                                                                           <  1:4120708/150‑1 (MQ=255)

                                                                         tcg‑atcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTg                                                                           <  2:4361251/150‑1 (MQ=255)

                                                                             ctcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGAt                                                                        <  2:3360279/149‑1 (MQ=255)

                                                                             atcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGAt                                                                        >  2:1113390/1‑150 (MQ=255)

                                                                              tcgCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATg                                                                       >  1:4206281/1‑150 (MQ=255)

                                                                                 cgAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGa                                                                    >  1:2901516/1‑150 (MQ=255)

                                                                                  gAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCa                                                                                                                                         >  1:5021956/1‑80 (MQ=37)

                                                                                  gAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGCTGCCCATGGCATCGCGCa                                                                                                                                         <  2:5021956/80‑1 (MQ=37)

                                                                                   tgTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGATATATCGCCCCCAGCTTCCTGCGCTATGCCTTGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAAc                                                                  <  1:4190375/149‑1 (MQ=255)

                                                                                     ttcttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCt                                                                >  1:2993990/1‑150 (MQ=255)

                                                                                       cttcCTGCTTTCAGAGGCGCACTACCCGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                              <  1:2172536/150‑1 (MQ=255)

                                                                                        ttcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTgg                                                             >  1:814273/1‑150 (MQ=255)

                                                                                         tcCTGCTTTCCGCGGCGCACCACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCATTGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGa                                                            <  1:3329003/150‑1 (MQ=255)

                                                                                         tcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCCGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGa                                                            <  1:3329104/150‑1 (MQ=255)

                                                                                         tcCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGa                                                            <  2:2613437/150‑1 (MQ=255)

                                                                                           cTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACa                                                          >  2:5433459/1‑150 (MQ=255)

                                                                                           cTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACa                                                          >  1:1006668/1‑150 (MQ=255)

                                                                                           cTGCTTTCAGAGGCGCACCACCAGCGCTATACCGCCCCCCGCTTCCTGCGCGAGGCCAGGGGCAGCGCCCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACa                                                          <  1:674602/150‑1 (MQ=255)

                                                                                            tGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAg                                                         >  2:1197528/1‑150 (MQ=255)

                                                                                            tGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAg                                                         >  2:4675660/1‑150 (MQ=255)

                                                                                             gCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGc                                                        >  1:2907150/1‑150 (MQ=255)

                                                                                                   agagGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGc                                                  >  2:2521396/1‑150 (MQ=255)

                                                                                                    gagGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCa                                                 >  1:112046/1‑150 (MQ=255)

                                                                                                     agGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                              <  2:4535821/136‑1 (MQ=255)

                                                                                                     agGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTg                                                              >  1:4535821/1‑136 (MQ=255)

                                                                                                       gcgcACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCa                                              >  2:2183462/1‑150 (MQ=255)

                                                                                                       gcgcACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCa                                              <  2:4272292/150‑1 (MQ=255)

                                                                                                       gcgcACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCa                                              >  1:5300980/1‑150 (MQ=255)

                                                                                                        cgcACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAg                                             >  2:3862798/1‑150 (MQ=255)

                                                                                                         gcACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGc                                            >  1:5160647/1‑150 (MQ=255)

                                                                                                         gcACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGc                                            <  2:580681/150‑1 (MQ=255)

                                                                                                         gcACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATACCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATc                                               >  2:2537010/1‑147 (MQ=255)

                                                                                                            cTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGc                                         >  2:4528579/1‑150 (MQ=255)

                                                                                                               ccAGCGCTATATCGCCCACAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTg                                      >  1:2223226/1‑150 (MQ=255)

                                                                                                                cAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGa                                     <  2:3685972/150‑1 (MQ=255)

                                                                                                                 aGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGAc                                    >  2:1736454/1‑150 (MQ=255)

                                                                                                                         atCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTccc                            >  1:994443/1‑150 (MQ=255)

                                                                                                                           cGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAg                          >  2:2823282/1‑150 (MQ=255)

                                                                                                                              ccccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCgg                                                                                                          >  1:3153617/1‑67 (MQ=255)

                                                                                                                                ccAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCa                     >  2:3998165/1‑150 (MQ=255)

                                                                                                                                  aGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACt                   >  1:2720690/1‑150 (MQ=255)

                                                                                                                                        cTGCGCGCTGCCCATGGCATCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGa             >  2:4760770/1‑150 (MQ=255)

                                                                                                                                         tGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGAt            >  2:1891978/1‑150 (MQ=255)

                                                                                                                                             cgATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACt        >  2:5204466/1‑150 (MQ=255)

                                                                                                                                               ttGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTc      <  1:3243685/149‑1 (MQ=255)

                                                                                                                                                   cATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGAt            >  2:4907954/1‑140 (MQ=255)

                                                                                                                                                   cATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGAt            <  1:4907954/140‑1 (MQ=255)

                                                                                                                                                   cATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGCtg  >  2:407251/1‑150 (MQ=255)

GACCTGGAGGCCTATGGCATGGTGTTGCTGCCGGAGTTCCGTGACTGGGAACTGCCCCAATTGCGCGATGCGATCG‑ATCGCGAGTTCTTCCTGCTTTCAGAGGCGCACTACCAGCGCTATATCGCCCCCAGCTTCCTGCGCGATGCCATGGGCAGCGCGCAGGGCCAGTGAGCCGGGCGAATCAGTCGCCGGTGTATTCGCAACCGCTGGTGCAGGTTTCGTGGATGCGAACCTTGGACAGCTCCGGCATCAGCGGCTTGACCTGGTCCCAGATCCACTTGGCGATCACTTCGCTG  >  NC_002947/2675000‑2675295

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: `‑`