BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,378,737	T→G	9.6%	A266A (GCA→GCC)	dnaK ←	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,378,737	0	T	G	9.6%	784.4 / 54.5	251	A266A (GCA→GCC)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand): ref base T (112/115); new base G (16/8); total (128/123)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 254 total aligned reads displayed.

GTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTC  >  NC_002947/5378588‑5378886

gTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCt                                                                                                                                                       <  2:8903241/150‑1 (MQ=255)

 tACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTg                                                                                                                                                      <  2:5326039/150‑1 (MQ=255)

 tACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTg                                                                                                                                                      >  2:8716791/1‑150 (MQ=255)

   cGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCg                                                                                                                                                    <  2:6125369/150‑1 (MQ=255)

   cGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCg                                                                                                                                                    <  2:5525701/150‑1 (MQ=255)

    gCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCgg                                                                                                                                                   >  2:4362623/1‑150 (MQ=255)

     cTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGc                                                                                                                                                  >  1:7800853/1‑150 (MQ=255)

        gACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGcttc                                                                                                                                               <  2:8103798/150‑1 (MQ=255)

             ggTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCa                                                                                                                                          <  1:638766/150‑1 (MQ=255)

             ggTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCa                                                                                                                                          <  2:6061108/150‑1 (MQ=255)

             ggTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCa                                                                                                                                          <  1:9120133/150‑1 (MQ=255)

              gTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAg                                                                                                                                         >  2:2466711/1‑150 (MQ=255)

               tCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGa                                                                                                                                        <  2:1731234/150‑1 (MQ=255)

                    gACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTg                                                                                                                                   >  2:4134214/1‑150 (MQ=255)

                    gACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTg                                                                                                                                   >  1:4326347/1‑150 (MQ=255)

                             cTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCag                                                                                                                          >  2:644052/1‑150 (MQ=255)

                              tCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCaga                                                                                                                         <  1:1964832/150‑1 (MQ=255)

                              tCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCaga                                                                                                                         >  2:246017/1‑150 (MQ=255)

                              tCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTATTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCaga                                                                                                                         <  2:8823606/150‑1 (MQ=255)

                              gccAGCGTGGCGCGGGAGATCGTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCaga                                                                                                                         <  2:3261311/149‑1 (MQ=255)

                               ccAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCagag                                                                                                                        >  1:2597095/1‑150 (MQ=255)

                                cAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAgg                                                                                                                       >  1:2120242/1‑150 (MQ=255)

                                cAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAgg                                                                                                                       >  1:4430754/1‑150 (MQ=255)

                                 aGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGTAGGGCCAGAGGa                                                                                                                      <  1:3479367/150‑1 (MQ=255)

                                 aGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGa                                                                                                                      <  1:5347011/150‑1 (MQ=255)

                                 aGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGa                                                                                                                      >  1:174949/1‑150 (MQ=255)

                                 aGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGa                                                                                                                      >  2:1540882/1‑150 (MQ=255)

                                        cgcgGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGttct                                                                                                               <  1:1046654/150‑1 (MQ=255)

                                         gcgGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGttctt                                                                                                              >  1:4932303/1‑150 (MQ=255)

                                          cgGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGttcttc                                                                                                             <  2:2743023/150‑1 (MQ=255)

                                             gagaTCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAgg                                                                                                          >  2:6229905/1‑150 (MQ=255)

                                             gagaTCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAgg                                                                                                          >  2:1629562/1‑150 (MQ=255)

                                                aTCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTcc                                                                                                       >  1:396177/1‑150 (MQ=255)

                                                aTCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTcc                                                                                                       <  2:5715918/150‑1 (MQ=255)

                                                   ttCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATg                                                                                                    <  1:5181214/150‑1 (MQ=255)

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

                                                     cACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGcc                                                                                                  >  2:4262028/1‑150 (MQ=255)

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

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

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

                                                          tCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGAct                                                                                             <  2:6461481/150‑1 (MQ=255)

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

                                                               tGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGTTTGACGTCGGTCGACTGAGCGGACGACATCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttc                                                                                        <  1:767436/150‑1 (MQ=255)

                                                               tGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttc                                                                                        <  2:1148816/150‑1 (MQ=255)

                                                               tGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttc                                                                                        <  1:6921719/150‑1 (MQ=255)

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

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

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

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

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

                                                                    aGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGaa                                                                                   <  1:1317955/150‑1 (MQ=255)

                                                                     ggACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAAc                                                                                  >  1:556197/1‑150 (MQ=255)

                                                                        ccGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcg                                                                               >  2:1437794/1‑150 (MQ=255)

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

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

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

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

                                                                                  cTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGt                                                                     >  1:1501505/1‑150 (MQ=255)

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

                                                                                     cAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGt                                                                  <  1:4065648/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                              aCGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCATGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTAGAACTCGTCGACAAGGTAGACGATCAGgc                                                         >  2:3877363/1‑150 (MQ=255)

                                                                                               cGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGgcg                                                        <  2:4627368/150‑1 (MQ=255)

                                                                                               cGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGgcg                                                        <  1:9059641/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                  gACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATc                                                                                                                    >  1:3669369/1‑71 (MQ=255)

                                                                                                                  gACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCa                                     <  2:5442057/150‑1 (MQ=255)

                                                                                                                    cTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCATGCGCATGTCGAAGTCTTCGCCAcc                                   >  2:7538570/1‑150 (MQ=255)

                                                                                                                     tGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGccg                                  >  2:4083978/1‑149 (MQ=255)

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

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

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

                                                                                                                          ggACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGac                             >  2:1992170/1‑149 (MQ=255)

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

                                                                                                                              gaaaGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGAACgtg                         <  1:6166721/147‑1 (MQ=255)

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

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

                                                                                                                                   cTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGcc                    >  2:4197559/1‑150 (MQ=255)

                                                                                                                                   cTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGcc                    >  2:5924192/1‑150 (MQ=255)

                                                                                                                                      gATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGtt                 >  1:8357220/1‑150 (MQ=255)

                                                                                                                                             gCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCTCCACCCAGGAAGGTGTCGCCGTTGGTAGcc          >  1:5986810/1‑150 (MQ=255)

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

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

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

GTACGCTTGACCAGGTCTTCGACCAGCGACTCCAGCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTC  >  NC_002947/5378588‑5378886

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