BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,378,815	A→C	8.9%	L240L (CTT→CTG)	dnaK ←	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,378,815	0	A	C	8.9%	394.8 / 20.9	129	L240L (CTT→CTG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand): ref base A (65/51); new base C (5/7); total (70/59)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

GCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGAT  >  NC_002947/5378666‑5378964

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

   tCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAAATCGTCGACAAgg                                                                                                                                                    <  2:23298/150‑1 (MQ=255)

      gCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAg                                                                                                                                                 <  1:1029398/150‑1 (MQ=255)

        aGTGATGTACGGCTGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTc                                                                                                                                               <  2:3627209/150‑1 (MQ=255)

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

         gTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCg                                                                                                                                              >  2:2617728/1‑150 (MQ=255)

            aTGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATc                                                                                                                                           <  2:2219493/150‑1 (MQ=255)

             tgaatggCAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCa                                                                                                                                          <  1:519848/145‑1 (MQ=255)

               tACGGCAGGTTGACGTCGGTCGACTGCGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAgg                                                                                                                                        <  2:4314436/150‑1 (MQ=255)

                atggCAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAgac                                                                                                                                       <  2:2458940/148‑3 (MQ=255)

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

                    cAGGTTGAGTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCAt                                                                                                                                   <  2:3046139/150‑1 (MQ=255)

                    cAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCAt                                                                                                                                   >  1:4874147/1‑150 (MQ=255)

                      ggTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGt                                                                                                                                 >  2:4785682/1‑150 (MQ=255)

                      ggTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGt                                                                                                                                 >  2:4785686/1‑150 (MQ=255)

                       gTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTc                                                                                                                                >  1:4123038/1‑150 (MQ=255)

                       gTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTc                                                                                                                                >  1:1249512/1‑150 (MQ=255)

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

                            tttcggtctGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGt                                                                                                                           >  2:2387316/3‑150 (MQ=255)

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

                               cggtcgACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCtt                                                                                                                        <  2:4321015/150‑1 (MQ=255)

                                ggtctgcTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTc                                                                                                                       >  2:3859707/7‑150 (MQ=255)

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

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

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

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

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

                                 gtcgACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCg                                                                                                                      <  1:4916886/150‑1 (MQ=255)

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

                                     gcTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCAccgc                                                                                                                  <  1:4859876/149‑3 (MQ=255)

                                     aCTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCAc                                                                                                                  >  2:4768506/1‑150 (MQ=255)

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

                                        gAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGccga                                                                                                               <  2:3635384/150‑3 (MQ=255)

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

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

                                             gacgacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAg                                                                                                          >  1:3157418/1‑150 (MQ=255)

                                               cgacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGgt                                                                                                        >  2:2272522/1‑150 (MQ=255)

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

                                                  cAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTc                                                                                                     <  1:3580490/150‑1 (MQ=255)

                                                  cAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTc                                                                                                     <  1:3602913/150‑1 (MQ=255)

                                                     cTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACCAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGcc                                                                                                  >  1:38534/1‑150 (MQ=255)

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

                                                         aTCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTg                                                                                              <  1:1072150/150‑1 (MQ=255)

                                                            ttGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTa                                                                                           >  2:3814998/1‑150 (MQ=255)

                                                            ttGGCATTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTa                                                                                           <  1:2989521/150‑1 (MQ=255)

                                                                       tGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTc                                                                                <  1:3094222/150‑1 (MQ=255)

                                                                       tGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTc                                                                                <  2:104775/150‑1 (MQ=255)

                                                                        gCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCg                                                                               <  1:1696829/150‑1 (MQ=255)

                                                                         cGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGa                                                                              <  2:4396404/150‑1 (MQ=255)

                                                                          ggCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGaa                                                                             <  1:1059373/150‑1 (MQ=255)

                                                                           gCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAAc                                                                            >  1:3357962/1‑150 (MQ=255)

                                                                            cttcttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACt                                                                           >  1:4163659/1‑150 (MQ=255)

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

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

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

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

                                                                              tcttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTAGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTgg                                                                         <  1:4355325/150‑1 (MQ=255)

                                                                               cttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGt                                                                        >  2:4476543/1‑150 (MQ=255)

                                                                                tttCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTa                                                                                    <  1:367927/137‑1 (MQ=255)

                                                                                tttCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTg                                                                       <  1:4460408/150‑1 (MQ=255)

                                                                                 ttCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACTAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGc                                                                      <  2:4509941/150‑1 (MQ=255)

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

                                                                                           gCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCAACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGa                                                            >  1:4514273/1‑150 (MQ=255)

                                                                                              gggCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACtt                                                         >  1:953223/1‑150 (MQ=255)

                                                                                              gggCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACtt                                                         >  1:5024946/1‑150 (MQ=255)

                                                                                                   agagGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCg                                                    >  2:4371694/1‑150 (MQ=255)

                                                                                                   agagGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCg                                                    >  1:4336907/1‑150 (MQ=255)

                                                                                                         tCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCg                                              <  1:2292565/150‑1 (MQ=255)

                                                                                                         tCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCg                                              >  1:421613/1‑150 (MQ=255)

                                                                                                           gTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGAt                                            >  1:1631777/1‑150 (MQ=255)

                                                                                                              cttcAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGAc                                         >  1:4978824/1‑150 (MQ=255)

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

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

                                                                                                                 cAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGa                                      >  1:2850226/1‑150 (MQ=255)

                                                                                                                 cAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGa                                      >  1:2850316/1‑150 (MQ=255)

                                                                                                                 cAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGa                                      >  1:4370899/1‑150 (MQ=255)

                                                                                                                    gTCCAGGCCCGACTCTTTCTTGAACTCGTCGACCAGGTAGTCGATCAGGCGCAGGTCGAAGTCTTCGCCTCCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGAGGACCGAAAc                                   <  1:2369038/150‑1 (MQ=255)

                                                                                                                         tGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGa                              <  1:1022115/150‑1 (MQ=255)

                                                                                                                          gCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGaa                             >  2:4881515/1‑150 (MQ=255)

                                                                                                                          gCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGaa                             <  1:2111529/150‑1 (MQ=255)

                                                                                                                          gCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGaa                             <  2:3309788/150‑1 (MQ=255)

                                                                                                                            ccGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAgg                           >  1:4590697/1‑150 (MQ=255)

                                                                                                                             cGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGt                          >  1:1233351/1‑150 (MQ=255)

                                                                                                                              gACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTa                         <  1:1275829/150‑1 (MQ=255)

                                                                                                                                  ctttcttGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACTGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTaccac                     <  1:275820/150‑1 (MQ=255)

                                                                                                                                       ttGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCAccc                >  1:3002031/1‑150 (MQ=255)

                                                                                                                                       ttGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCAccc                <  2:3094274/150‑1 (MQ=255)

                                                                                                                                         gAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAg              >  2:1393498/1‑150 (MQ=255)

                                                                                                                                          aaCTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAgg             >  2:3600582/1‑150 (MQ=255)

                                                                                                                                           aCTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGt            >  2:3881375/1‑150 (MQ=255)

                                                                                                                                           aCTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGt            >  2:4314731/1‑150 (MQ=255)

                                                                                                                                            cTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTc           <  1:4757423/150‑1 (MQ=255)

                                                                                                                                            cTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTc           <  1:4097614/150‑1 (MQ=255)

                                                                                                                                            cTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTc           <  1:4097592/150‑1 (MQ=255)

                                                                                                                                             tcgtcgACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCa          >  2:2624731/1‑150 (MQ=255)

                                                                                                                                             tcgtcgACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCa          >  2:545885/1‑150 (MQ=255)

                                                                                                                                                tcgACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATaa       >  1:2908950/1‑150 (MQ=255)

                                                                                                                                                  gACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAAc     >  2:2417433/1‑150 (MQ=255)

                                                                                                                                                   aCAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACg    >  1:423735/1‑150 (MQ=255)

                                                                                                                                                     aaGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGAt  <  1:1887840/150‑1 (MQ=255)

GCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCGGCGATTTCGATGACCGAAACGTCGAAGGTACCACCACCCAGGTCATAAACGAT  >  NC_002947/5378666‑5378964

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