BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	593,390	C→G	7.2%	A375A (GCC→GCG)	PP_0503 →	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	593,390	0	C	G	7.2%	401.1 / 18.2	125	A375A (GCC→GCG)	PP_0503	MFS transporter
Reads supporting (aligned to +/- strand): ref base C (57/59); new base G (6/3); total (63/62)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

Only 100 of 126 total aligned reads displayed.

GCCGGACGTGGTGCGATGGCTAGCCGTTTGATCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAA  >  NC_002947/593241‑593535

               aTGGCTAGCCGTTTGATCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTAc                                                                                                                                    >  2:718003/1‑150 (MQ=255)

                tGGCTAGCCGTTTGATCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACt                                                                                                                                   >  1:948015/1‑150 (MQ=255)

                  gCTAGCCGTTTGATCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTgg                                                                                                                                 <  2:759465/150‑1 (MQ=255)

                  gCTAGCCGTTTGATCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTgg                                                                                                                                 <  1:1742714/150‑1 (MQ=255)

                   cTAGCCGTTTGATCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCGGCGGCACtgatgccta                                                                                                                                   >  2:2050523/1‑140 (MQ=255)

                   cTAGCCGTTTGATCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCGGCGGCACtgatgccta                                                                                                                                   <  1:2050523/147‑8 (MQ=255)

                   cTAGCCGTTTGATCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCGGCGGCACtgatgccca                                                                                                                                   >  2:2050735/1‑140 (MQ=255)

                   cTAGCCGTTTGATCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCGGCGGCACtgatgccca                                                                                                                                   <  1:2050735/147‑8 (MQ=255)

                     aGCCGTTTGATCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTa                                                                                                                              <  1:455667/150‑1 (MQ=255)

                     aGCCGTTTGATCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTa                                                                                                                              <  1:453948/150‑1 (MQ=255)

                      gCCGTTTGATCGGCGCATACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTAt                                                                                                                             >  2:215758/1‑150 (MQ=255)

                       ccGTTTGAGCGGCGCCTACATGGGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTCGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATTAGTGCCGCCGCACGGTGCTACTGGCTATg                                                                                                                            >  2:2121613/1‑150 (MQ=255)

                              aTCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGACTGGGTTGCGGTGTGACCGTCGACCAGGCGCAACAGTTCCGCCGCACTGTGCTACTGGCTATGGCTGCaa                                                                                                                     <  1:421707/150‑1 (MQ=255)

                               tCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAAc                                                                                                                    >  1:941754/1‑150 (MQ=255)

                               tCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑CGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAAc                                                                                                                    >  1:2289688/1‑150 (MQ=255)

                                 ggCGCCTACATGTGTTGCTATTCGTTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAAcgc                                                                                                                  <  1:1238724/150‑1 (MQ=255)

                                 ggCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTa                                                                                                                                     >  2:875573/1‑131 (MQ=255)

                                 ggCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTa                                                                                                                                     <  1:875573/131‑1 (MQ=255)

                                  gcgcCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGAGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAAcgcg                                                                                                                 >  1:748112/1‑150 (MQ=255)

                                    gcCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAAcgcgcg                                                                                                               >  2:2232776/1‑150 (MQ=255)

                                       tACATGTGTTGCTATTCGCTGGGCAGCGGTTTGGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAg                                                                                                             >  1:117478/1‑150 (MQ=255)

                                        aCATGTGTTGCTATTCGCTGGGCAGCGGTTTGGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAgg                                                                                                            <  2:117478/150‑1 (MQ=255)

                                        aCATGTGTTGCTATCCGCTGGGCAGCGGGTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAggg                                                                                                           >  2:1956010/1‑150 (MQ=255)

                                         cATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGggc                                                                                                          <  1:718003/150‑1 (MQ=255)

                                         cATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGggc                                                                                                          >  1:2336305/1‑150 (MQ=255)

                                            gtgtTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGggcggc                                                                                                       <  1:843126/150‑1 (MQ=255)

                                             tgtTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGAGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGcc                                                                                                      >  1:252298/1‑150 (MQ=255)

                                               ttGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGa                                                                                                    >  2:112953/1‑150 (MQ=255)

                                               ttGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGa                                                                                                    >  1:1511342/1‑150 (MQ=255)

                                                tGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAg                                                                                                   <  2:2156378/150‑1 (MQ=255)

                                                 gCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGa                                                                                                              >  1:1930137/1‑138 (MQ=255)

                                                 gCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGa                                                                                                              <  1:752905/138‑1 (MQ=255)

                                                 gCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGa                                                                                                              >  2:752905/1‑138 (MQ=255)

                                                  cTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGAACTGATTTTGCACACCGACCCCGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTAGTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgc                                                                                                 >  2:1289444/1‑150 (MQ=255)

                                                   tATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACCCGCGAGGGCGGCCGAGgcc                                                                                                >  2:244805/1‑150 (MQ=255)

                                                    aTTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgccg                                                                                               >  2:2406892/1‑150 (MQ=255)

                                                     ttCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgccgc                                                                                              <  2:723380/150‑1 (MQ=255)

                                                     ttCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgccgc                                                                                              >  1:1850886/1‑150 (MQ=255)

                                                       cGCTGGGCAGCGGTTT‑GGGGGAGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCAGCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGccc                                                                                            <  1:106745/150‑1 (MQ=255)

                                                        gCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCt                                                                                           <  1:112953/150‑1 (MQ=255)

                                                         cTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTg                                                                                          >  1:1472994/1‑150 (MQ=255)

                                                            ggCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCaa                                                                                       >  1:274882/1‑150 (MQ=255)

                                                                gCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGtt                                                                                   <  1:496042/150‑1 (MQ=255)

                                                                gCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGtt                                                                                   <  2:1175784/150‑1 (MQ=255)

                                                                  ggTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCa                                                                                 <  2:1472994/150‑1 (MQ=255)

                                                                   gTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAg                                                                                <  2:1579519/150‑1 (MQ=255)

                                                                   gTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAg                                                                                <  2:1838901/150‑1 (MQ=255)

                                                                   gTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGAGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAg                                                                                <  1:899291/150‑1 (MQ=255)

                                                                     tt‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGTCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGt                                                                              <  1:2040868/150‑1 (MQ=255)

                                                                            gCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAAc                                                                        >  1:872788/1‑150 (MQ=255)

                                                                              ggTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACtt                                                                      <  2:2001560/150‑1 (MQ=255)

                                                                              ggTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACtt                                                                      <  2:2289688/150‑1 (MQ=255)

                                                                              ggTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACtt                                                                      <  1:215758/150‑1 (MQ=255)

                                                                                      cGAGCTGGGTGTGTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCaa                                                              <  2:2301283/150‑1 (MQ=255)

                                                                                           tGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAgaccg                                                         >  1:2200597/1‑150 (MQ=255)

                                                                                             ggTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAgaccgac                                                       <  2:1850886/150‑1 (MQ=255)

                                                                                             ggTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAgaccgac                                                       >  1:2106675/1‑150 (MQ=255)

                                                                                             ggTGTTTGCACACTGGGGCTGGGCTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAgaccgac                                                       >  1:327490/1‑150 (MQ=255)

                                                                                               tgtTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCa                                                     >  2:1994785/1‑150 (MQ=255)

                                                                                               tgtTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCa                                                     >  2:2086480/1‑150 (MQ=255)

                                                                                                gtTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCAt                                                    <  2:1511342/150‑1 (MQ=255)

                                                                                                 tttGCACACTGGGGCTGGGTTGCGGTGTGGGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATg                                                   <  2:252298/150‑1 (MQ=255)

                                                                                                      acacTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCTCGGCCGccc                                                                                            <  1:2005361/104‑1 (MQ=255)

                                                                                                       cacTGGGGCTGGGTTGCGGTGTGTGGGCTCGGCATCAGCATCAGTGCCCCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATAcc                                             <  2:2200597/150‑1 (MQ=255)

                                                                                                        acTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATAccc                                            >  2:1567481/1‑150 (MQ=255)

                                                                                                        acTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATAGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATAccc                                            >  2:1004727/1‑150 (MQ=255)

                                                                                                           ggggCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACg                                         <  1:2232776/150‑1 (MQ=255)

                                                                                                             ggCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGTCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAg                                                                                <  2:353772/109‑1 (MQ=255)

                                                                                                             ggCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGTCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAg                                                                                >  1:353772/1‑109 (MQ=255)

                                                                                                               cTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                     <  1:1569665/150‑1 (MQ=255)

                                                                                                               cTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                     <  2:2005148/150‑1 (MQ=255)

                                                                                                                tGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGcc                                                                <  2:181077/122‑1 (MQ=255)

                                                                                                                    ttGCGGTGTGTGGGCTTGGCATGGGCATCAGTACCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgt                                >  2:1327513/1‑150 (MQ=255)

                                                                                                                            tgtgGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgtcg                        >  1:1544196/1‑150 (MQ=255)

                                                                                                                              tgGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACggg                                       >  1:1019046/1‑133 (MQ=255)

                                                                                                                              tgGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACggg                                       <  2:1019046/133‑1 (MQ=255)

                                                                                                                              tgGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgtcgac                      >  2:1287430/1‑150 (MQ=255)

                                                                                                                                 gCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCTCGGCCGCCCtcg                                                                                         <  2:1251663/80‑3 (MQ=39)

                                                                                                                                 gCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCTCGGCCGCCCtcg                                                                                         >  1:1251663/1‑78 (MQ=39)

                                                                                                                                 gCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCtt                   >  1:2119673/1‑150 (MQ=255)

                                                                                                                                  cTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTg                  <  2:2119673/150‑1 (MQ=255)

                                                                                                                                     ggCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGAcccgtggg               <  2:375900/150‑7 (MQ=255)

                                                                                                                                       cATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCtt             >  1:1224498/1‑150 (MQ=255)

                                                                                                                                        aTGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTg            >  1:915780/1‑150 (MQ=255)

                                                                                                                                         tGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                     <  2:1922939/124‑1 (MQ=255)

                                                                                                                                         tGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                     >  1:1922939/1‑124 (MQ=255)

                                                                                                                                         tGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGc           <  1:1415552/150‑1 (MQ=255)

                                                                                                                                         tGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATCGGTAAACCGCGAGTGGGTTGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGc           >  1:1623127/1‑150 (MQ=255)

                                                                                                                                          gggCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCt          >  1:1800424/1‑150 (MQ=255)

                                                                                                                                           ggCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATAccc                                            <  1:958778/115‑1 (MQ=255)

                                                                                                                                           ggCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATAccc                                            >  2:958778/1‑115 (MQ=255)

                                                                                                                                             cATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGTGCGATGCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTa       <  1:1956010/150‑1 (MQ=255)

                                                                                                                                              cacAGTGCGGCGGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAg      >  2:383401/3‑150 (MQ=255)

                                                                                                                                              cacAGTGCGGCGGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGGCGACCTTGGTCTTGCTGTAg      >  2:560089/3‑150 (MQ=255)

                                                                                                                                              aTCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAg      <  2:948015/150‑1 (MQ=255)

                                                                                                                                               acAGTGCGGCGGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTa                                                                   >  2:1018728/2‑88 (MQ=255)

                                                                                                                                               acAGTGCGGCGGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTa                                                                   <  1:1018728/87‑1 (MQ=255)

                                                                                                                                               acAGTGCGGCGGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGc     >  2:2367677/2‑150 (MQ=255)

                                                                                                                                                 aGTGCCGCCGGACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCa   >  2:1355649/1‑150 (MQ=255)

                                                                                                                                                  gTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCaa  >  2:1030843/1‑150 (MQ=255)

GCCGGACGTGGTGCGATGGCTAGCCGTTTGATCGGCGCCTACATGTGTTGCTATTCGCTGGGCAGCGGTTT‑GGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAA  >  NC_002947/593241‑593535

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`