BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,378,785	C→G	9.8%	M250I (ATG→ATC)	dnaK ←	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,378,785	0	C	G	9.8%	774.1 / 56.4	242	M250I (ATG→ATC)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand): ref base C (124/93); new base G (16/8); total (141/101)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

ATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTC  >  NC_002947/5378636‑5378934

   ttCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATg                                                                                                                                                    >  1:6237715/1‑150 (MQ=255)

   ttCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATg                                                                                                                                                    >  2:1462018/1‑150 (MQ=255)

    tCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGc                                                                                                                                                   >  2:6568592/1‑150 (MQ=255)

    tCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGc                                                                                                                                                   >  2:6320263/1‑150 (MQ=255)

      aCGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGccc                                                                                                                                                 <  1:7788008/150‑1 (MQ=255)

      aCGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGccc                                                                                                                                                 <  2:2521906/150‑1 (MQ=255)

        gTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGa                                                                                                                                               >  1:1762412/1‑150 (MQ=255)

        gTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGa                                                                                                                                               <  1:2471261/150‑1 (MQ=255)

        gTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGa                                                                                                                                               >  2:4614996/1‑150 (MQ=255)

        gTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGa                                                                                                                                               <  2:3559246/150‑1 (MQ=255)

          caaGGTGCTTAGGACCGGTCTGGTCAGCCGTGATGAATGGCAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGccctgtt                                                                                                                                             <  1:1601977/148‑5 (MQ=255)

                gCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttct                                                                                                                                       >  1:1095387/1‑150 (MQ=255)

                gCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttct                                                                                                                                       <  1:6737932/150‑1 (MQ=255)

                gCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttct                                                                                                                                       >  2:8610004/1‑150 (MQ=255)

                gCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttct                                                                                                                                       <  2:5832112/150‑1 (MQ=255)

                 cTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttctt                                                                                                                                      >  1:2846382/1‑150 (MQ=255)

                 cTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttctt                                                                                                                                      <  1:3938826/150‑1 (MQ=255)

                  ttAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTg                                                                                                                                     >  2:2741783/1‑150 (MQ=255)

                  ttAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTg                                                                                                                                     >  2:1080392/1‑150 (MQ=255)

                  ttAGGACCGGTCTGGTCAGCCGTGATGAATGGCAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCAGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTg                                                                                                                                     >  1:796676/1‑150 (MQ=255)

                      gACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCTGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACt                                                                                                                                 <  1:3962401/150‑1 (MQ=255)

                      gACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTg                                                                                                                                     >  2:569258/1‑146 (MQ=255)

                      gACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTg                                                                                                                                     <  1:5143285/146‑1 (MQ=255)

                         cGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcgt                                                                                                                              >  2:7189870/1‑150 (MQ=255)

                         cGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcgt                                                                                                                              <  1:1057863/150‑1 (MQ=255)

                          ggTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcgtc                                                                                                                             <  2:777677/150‑1 (MQ=255)

                          ggTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcgtc                                                                                                                             <  2:2115941/150‑1 (MQ=255)

                               cATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACaa                                                                                                                        <  2:385201/150‑1 (MQ=255)

                                aTCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAg                                                                                                                       >  1:6556726/1‑150 (MQ=255)

                                 tCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAgg                                                                                                                      >  1:3614139/1‑150 (MQ=255)

                                   tGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTa                                                                                                                    <  1:7629288/150‑1 (MQ=255)

                                   tGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTa                                                                                                                    >  1:7446571/1‑150 (MQ=255)

                                   tGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTa                                                                                                                    <  1:8880306/150‑1 (MQ=255)

                                    gCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAg                                                                                                                   >  2:2535195/1‑150 (MQ=255)

                                      aGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTc                                                                                                                 >  2:8349439/1‑150 (MQ=255)

                                      aGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTc                                                                                                                 >  1:6827978/1‑150 (MQ=255)

                                        tGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGa                                                                                                               <  2:4446767/150‑1 (MQ=255)

                                             tACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAgg                                                                                                          >  2:1279916/1‑150 (MQ=255)

                                             tACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAgg                                                                                                          <  2:8272122/150‑1 (MQ=255)

                                             tACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTAGAACTCGTCGACAAGGTAGTCGATCAgg                                                                                                          >  2:1681487/1‑150 (MQ=255)

                                              aCGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGgc                                                                                                         >  1:2574906/1‑150 (MQ=255)

                                                ggCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGgcgc                                                                                                       <  1:2261461/150‑1 (MQ=255)

                                                ggCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTGGGCTTTTTCTGCGGCTTCTTTCCGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGgcgc                                                                                                       <  2:5140726/150‑1 (MQ=255)

                                                 gCAGGTTGAATTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTATCTAAAGATCGATGCCCTGTTCCTTCTTGAACTCGCCAACCAGATATTCGACCAGACGCa                                                                                                      >  2:3631120/1‑150 (MQ=255)

                                                  cAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCAt                                                                                                     <  2:8757676/150‑1 (MQ=255)

                                                  cAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCATTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCAt                                                                                                     <  1:4554274/150‑1 (MQ=255)

                                                       tGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGa                                                                                                >  2:9036881/1‑150 (MQ=255)

                                                        gACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGaa                                                                                               <  1:1622417/150‑1 (MQ=255)

                                                        gACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGaa                                                                                               <  1:5841208/150‑1 (MQ=255)

                                                          tttcggtctGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGt                                                                                             >  1:4100431/3‑150 (MQ=255)

                                                           gtcggtcgACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTc                                                                                            >  2:4248822/1‑150 (MQ=255)

                                                              ggtcgACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTc                                                                                         <  1:678746/150‑1 (MQ=255)

                                                              ggtcgACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTc                                                                                         <  2:785675/150‑1 (MQ=255)

                                                              ggtcgACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTc                                                                                         <  1:9002043/150‑1 (MQ=255)

                                                                  gACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCa                                                                                     >  1:4025569/1‑150 (MQ=255)

                                                                   aCTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCAc                                                                                    <  2:1501505/150‑1 (MQ=255)

                                                                      gAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGccga                                                                                 <  1:1500389/150‑3 (MQ=255)

                                                                       aGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAg                                                                                <  2:5288737/150‑1 (MQ=255)

                                                                          ggACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGaa                                                                             <  2:4932303/150‑1 (MQ=255)

                                                                          ggACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGaa                                                                             <  1:246245/150‑1 (MQ=255)

                                                                          ggACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGaa                                                                             >  1:574914/1‑150 (MQ=255)

                                                                          ggACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGaa                                                                             >  2:5944857/1‑150 (MQ=255)

                                                                           gacgaaAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGaac                                                                            <  1:609862/150‑2 (MQ=255)

                                                                              gacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGgtg                                                                         >  2:2806260/1‑150 (MQ=255)

                                                                              gacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGgtg                                                                         >  2:5209913/1‑150 (MQ=255)

                                                                              gacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGgtg                                                                         >  1:377046/1‑150 (MQ=255)

                                                                              gacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGgtg                                                                         >  1:3152471/1‑150 (MQ=255)

                                                                              gaaaGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGAACgtg                                                                         >  1:3546664/4‑150 (MQ=255)

                                                                              gaaaGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGAACgtg                                                                         >  2:2054099/4‑150 (MQ=255)

                                                                              gaaaGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGAACgtg                                                                         >  1:5205259/4‑150 (MQ=255)

                                                                              gaaaGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGAACgtg                                                                         >  1:8220492/4‑150 (MQ=255)

                                                                                  gCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGc                                                                     >  1:8644945/1‑150 (MQ=255)

                                                                                  gCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGc                                                                     >  2:6495129/1‑150 (MQ=255)

                                                                                     cGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGt                                                                  <  2:4326347/150‑1 (MQ=255)

                                                                                        tCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTgg                                                               <  2:5463934/150‑1 (MQ=255)

                                                                                             gCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGcc                                                          >  1:6786658/1‑150 (MQ=255)

                                                                                                           ttcttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTg                                            >  2:8708390/1‑150 (MQ=255)

                                                                                                           ttcttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTg                                            >  1:7323784/1‑150 (MQ=255)

                                                                                                           ttcttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTg                                            <  1:1540882/150‑1 (MQ=255)

                                                                                                            tcttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTgg                                           >  1:7571027/1‑150 (MQ=255)

                                                                                                              tttCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTg                                         >  2:5907502/1‑150 (MQ=255)

                                                                                                                tCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCt                                       >  1:9126701/1‑150 (MQ=255)

                                                                                                                      gCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGt                                 >  2:8586871/1‑150 (MQ=255)

                                                                                                                      gCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCAAGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGt                                 >  2:1765015/1‑150 (MQ=255)

                                                                                                                        tGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCg                               >  2:3674670/1‑150 (MQ=255)

                                                                                                                           aGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACt                            >  2:3105507/1‑150 (MQ=255)

                                                                                                                               ccAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCgg                        >  2:4189204/1‑150 (MQ=255)

                                                                                                                               ccAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCgg                        >  2:4066282/1‑150 (MQ=255)

                                                                                                                               ccAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCgg                        >  1:4350209/1‑150 (MQ=255)

                                                                                                                               ccAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGGTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCgg                        >  2:9100293/1‑150 (MQ=255)

                                                                                                                                cAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGc                       >  1:221315/1‑150 (MQ=255)

                                                                                                                                cAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGc                       >  1:1270929/1‑150 (MQ=255)

                                                                                                                                   agGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTATTTCGAACTGGTGCTCACCGTCGACTTCGGCGAt                    >  1:1576919/1‑150 (MQ=255)

                                                                                                                                     gATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGAttt                  >  1:1924788/1‑150 (MQ=255)

                                                                                                                                        cGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGa               <  1:1462018/150‑1 (MQ=255)

                                                                                                                                        cGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGa               >  1:7067481/1‑150 (MQ=255)

                                                                                                                                        cGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGa               <  1:644052/150‑1 (MQ=255)

                                                                                                                                          ttcttcAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTCGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATg             >  1:8343273/1‑150 (MQ=255)

                                                                                                                                              tcAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCg         >  2:1345635/1‑150 (MQ=255)

                                                                                                                                                     cATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTc  >  1:4153683/1‑150 (MQ=255)

ATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTC  >  NC_002947/5378636‑5378934

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