BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,036,872	T→A	5.0%	L5* (TTG→TAG) ‡	PP_0889 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,036,872	0	T	A	5.0%	1284.1 / 36.3	371	L5* (TTG→TAG) ‡	PP_0889	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (144/208); new base A (7/12); total (151/220)
Fisher's exact test for biased strand distribution p-value = 8.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 382 total aligned reads displayed.

CAATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGACTGGCAG  >  NC_002947/1036723‑1037020

cAATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCtt                                                                                                                                                      >  1:7446687/1‑150 (MQ=255)

 aaTATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCCGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTCATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGACTCCGCTCCTGCTACTGAGAGCGCGGGCTTg                                                                                                                                                     >  1:204154/1‑150 (MQ=255)

 aaTATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTg                                                                                                                                                     >  1:6774267/1‑150 (MQ=255)

 aaTATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGACTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGTGCGCGGGCTTg                                                                                                                                                     >  1:5637729/1‑150 (MQ=255)

     tctcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcg                                                                                                                                                 >  1:4431819/1‑150 (MQ=255)

      ctcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgc                                                                                                                                                <  1:3050722/150‑1 (MQ=255)

       tcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcgcg                                                                                                                                               >  2:3950352/1‑150 (MQ=255)

       tcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcgcg                                                                                                                                               >  1:6700666/1‑150 (MQ=255)

       tcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcgcg                                                                                                                                               >  1:1721187/1‑150 (MQ=255)

       tcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcg                                                                                                                                               <  2:1769597/150‑1 (MQ=255)

        cGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCATAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgc                                                                                                                                              <  1:4054151/150‑2 (MQ=255)

        cGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgc                                                                                                                                              <  1:941366/150‑2 (MQ=255)

         ggCCACCGCCCGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGaa                                                                                                                                             <  1:2870458/150‑1 (MQ=255)

          gCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAg                                                                                                                                            >  2:5277968/1‑150 (MQ=255)

           ccAGCGCCCGTGCGGTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAgg                                                                                                                                           >  1:4690757/1‑150 (MQ=255)

            cACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAggg                                                                                                                                          <  2:2219677/150‑1 (MQ=255)

              ccGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCt                                                                                                                                        <  2:4812629/150‑1 (MQ=255)

               cGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTg                                                                                                                                       <  1:693747/150‑1 (MQ=255)

               cGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTg                                                                                                                                       >  2:416654/1‑150 (MQ=255)

                   cGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCaaa                                                                                                                                   <  2:4645808/150‑1 (MQ=255)

                      gcgcTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCa                                                                                                                                <  1:1038386/150‑1 (MQ=255)

                       cgcTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAg                                                                                                                               >  1:6838669/1‑150 (MQ=255)

                        gcTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGt                                                                                                                              >  2:169813/1‑150 (MQ=255)

                            ggCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCa                                                                                                                          <  2:2226512/150‑1 (MQ=255)

                               aTGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATc                                                                                                                       >  2:5202661/1‑150 (MQ=255)

                               aTGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATc                                                                                                                       <  1:3264274/150‑1 (MQ=255)

                                 gCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCAAAGGGCTGCAAAGCAGTCCCAATCaa                                                                                                                     <  1:6108722/150‑1 (MQ=255)

                                   aCCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAgg                                                                                                                   >  1:541061/1‑150 (MQ=255)

                                      gacggacgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGa                                                                                                                >  1:3526094/1‑150 (MQ=255)

                                           acgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCg                                                                                                           >  2:2474862/1‑150 (MQ=255)

                                            cgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGa                                                                                                          <  2:3173145/150‑1 (MQ=255)

                                            cgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGa                                                                                                          <  2:6912140/150‑1 (MQ=255)

                                             attGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGAt                                                                                                         <  2:5733997/149‑1 (MQ=255)

                                                gCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGaa                                                                                                      >  2:2907793/1‑150 (MQ=255)

                                                gCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGaa                                                                                                      >  1:7074832/1‑150 (MQ=255)

                                                 ccgcgcAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                                     <  1:7118637/148‑1 (MQ=255)

                                                 cAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                                     <  1:3362410/150‑1 (MQ=255)

                                                 cAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                                     <  1:6240794/150‑1 (MQ=255)

                                                 cAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                                     <  2:512128/150‑1 (MQ=255)

                                                 cAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                                     <  1:964063/150‑1 (MQ=255)

                                                 cAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                                     <  1:3665215/150‑1 (MQ=255)

                                                   gcgcAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGt                                                                                                   >  2:283468/1‑150 (MQ=255)

                                                       aaGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTg                                                                                               >  2:2108623/1‑150 (MQ=255)

                                                       aaGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTg                                                                                               <  1:5935007/150‑1 (MQ=255)

                                                        aGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGt                                                                                              <  1:4497566/150‑1 (MQ=255)

                                                        aGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGt                                                                                              >  1:3197138/1‑150 (MQ=255)

                                                          cTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCTCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAATTCGAGCTCGATGAAGGTACTGTTg                                                                                            <  2:4415864/150‑1 (MQ=255)

                                                          cTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                            <  2:6860757/150‑1 (MQ=255)

                                                          cTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                            <  1:5216415/150‑1 (MQ=255)

                                                           tGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTgg                                                                                           <  1:6399711/150‑1 (MQ=255)

                                                             cAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGaa                                                                                                                                             <  2:6422313/98‑1 (MQ=255)

                                                                aaGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATc                                                                                      >  2:4041915/1‑150 (MQ=255)

                                                                 aGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATcc                                                                                     <  2:5812243/150‑1 (MQ=255)

                                                                  gCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCt                                                                                    >  1:2234828/1‑150 (MQ=255)

                                                                    ggCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCCAGGCGAGCTCGATGAAGGTCCTGTTGGTCATCCTTg                                                                                  >  1:5814001/1‑150 (MQ=255)

                                                                      ccGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgcg                                                                                >  1:4864514/1‑150 (MQ=255)

                                                                         gTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTg                                                                             >  1:451525/1‑150 (MQ=255)

                                                                               cgcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCa                                                                       >  2:1213588/1‑150 (MQ=255)

                                                                                gcTGATTCCTTTTCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                                      <  1:2929219/150‑1 (MQ=255)

                                                                                gcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                                      <  2:1545900/150‑1 (MQ=255)

                                                                                gcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                                      >  2:7274677/1‑150 (MQ=255)

                                                                                   gATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTa                                                                   <  2:2723232/150‑1 (MQ=255)

                                                                                     ttCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAgc                                                                 <  1:1384574/150‑1 (MQ=255)

                                                                                      tCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAgcg                                                                >  1:2476241/1‑150 (MQ=255)

                                                                                         ttATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCa                                                             >  2:407693/1‑150 (MQ=255)

                                                                                         ttATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCa                                                             >  2:1360896/1‑150 (MQ=255)

                                                                                            tCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAg                                                          >  2:4566416/1‑150 (MQ=255)

                                                                                             cGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGc                                                         <  2:5431800/150‑1 (MQ=255)

                                                                                              gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCt                                                        >  2:431108/1‑150 (MQ=255)

                                                                                              gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCt                                                        <  2:766231/150‑1 (MQ=255)

                                                                                               atgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTg                                                       >  1:6789744/1‑150 (MQ=255)

                                                                                               atgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTg                                                       >  1:3924364/1‑150 (MQ=255)

                                                                                                tgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGGTGCAGGTAGCGCCAAAGCTGa                                                      <  2:6992658/150‑1 (MQ=255)

                                                                                                 gatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGaa                                                     <  2:1995499/150‑1 (MQ=255)

                                                                                                 gatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGaa                                                     <  1:4414132/150‑1 (MQ=255)

                                                                                                   tgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGa                                                   >  2:5525188/1‑150 (MQ=255)

                                                                                                      tGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTg                                                <  2:3406353/150‑1 (MQ=255)

                                                                                                       ggCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgct                                                                                                                                             <  1:5498026/56‑3 (MQ=255)

                                                                                                        gCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                              <  1:189080/150‑1 (MQ=255)

                                                                                                         ccTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGATCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGcc                                             >  1:7195926/1‑150 (MQ=255)

                                                                                                                gcgGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCg                                      >  1:2583693/1‑150 (MQ=255)

                                                                                                                  gggggACTCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCa                                    <  2:7171317/150‑1 (MQ=255)

                                                                                                                  gggggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCa                                    <  2:6011285/150‑1 (MQ=255)

                                                                                                                    gggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGt                                  <  1:3798854/150‑1 (MQ=255)

                                                                                                                     ggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTg                                 <  2:6565393/150‑1 (MQ=255)

                                                                                                                      gAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCGAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTgg                                <  1:3898415/150‑1 (MQ=255)

                                                                                                                       aaTCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGc                               <  1:1620933/150‑1 (MQ=255)

                                                                                                                         tCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAg                             >  2:591292/1‑150 (MQ=255)

                                                                                                                           cGCTCCTGCAACTGTGAACGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAg                                                                  >  1:1288885/1‑111 (MQ=255)

                                                                                                                            gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCaaa                                                           >  2:4920433/1‑117 (MQ=255)

                                                                                                                            gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCaaa                                                           <  1:4920433/117‑1 (MQ=255)

                                                                                                                            gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAAc                          >  2:3879002/1‑150 (MQ=255)

                                                                                                                                 tGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGAccc                     >  2:3701230/1‑150 (MQ=255)

                                                                                                                                  gCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCt                    >  1:3434982/1‑150 (MQ=255)

                                                                                                                                        gtgAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCaaa                                                           >  2:4250931/1‑105 (MQ=255)

                                                                                                                                        cccttcgcgGGCAAGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCaaa                                                           >  1:4325055/6‑105 (MQ=38)

                                                                                                                                           aGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTc           <  2:1847205/150‑1 (MQ=255)

                                                                                                                                            gcgcGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCg          >  1:4607958/1‑150 (MQ=255)

                                                                                                                                             cgcgGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGa         >  2:2354450/1‑150 (MQ=255)

                                                                                                                                               cgGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGACt       <  1:117936/150‑1 (MQ=255)

CAATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGACTGGCAG  >  NC_002947/1036723‑1037020

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