BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,637,084	C→T	6.8%	P283S (CCT→TCT) ‡	lepA →	elongation factor 4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,637,084	0	C	T	6.8%	734.5 / 28.3	236	P283S (CCT→TCT) ‡	lepA	elongation factor 4
Reads supporting (aligned to +/- strand): ref base C (136/84); new base T (7/9); total (143/93)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 240 total aligned reads displayed.

GTGCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAA  >  NC_002947/1636937‑1637233

gTGCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCt                                                                                                                                                     <  1:2900793/150‑1 (MQ=255)

  gCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGa                                                                                                                                                   >  1:2739357/1‑150 (MQ=255)

  gCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGa                                                                                                                                                   <  1:4982848/150‑1 (MQ=255)

  gCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGa                                                                                                                                                   <  2:672892/150‑1 (MQ=255)

    aCTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAgg                                                                                                                                                 >  2:3022827/1‑150 (MQ=255)

        ggTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGa                                                                                                                                             <  2:2491211/150‑1 (MQ=255)

        ggTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGa                                                                                                                                             <  1:4232942/150‑1 (MQ=255)

        ggTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCACCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGa                                                                                                                                             <  1:4753487/150‑1 (MQ=255)

         gTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGAg                                                                                                                                            <  2:2131157/150‑1 (MQ=255)

         gTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAg                                                                                                                                            <  1:956106/150‑1 (MQ=255)

             aCAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGAggaca                                                                                                                                        <  2:2985757/150‑4 (MQ=255)

             aCAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGAggaca                                                                                                                                        <  1:2536585/150‑4 (MQ=255)

              cAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGtt                                                                                                                                       <  1:201475/150‑1 (MQ=255)

              cAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGtt                                                                                                                                       >  1:2609883/1‑150 (MQ=255)

                gCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGc                                                                                                                                     <  1:5468070/150‑1 (MQ=255)

                     ggTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGc                                                                                                                                <  2:2557993/150‑1 (MQ=255)

                       tgtCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCtt                                                                                                                              <  1:103305/150‑1 (MQ=255)

                         tCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCa                                                                                                                            <  2:3160439/150‑1 (MQ=255)

                           ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaag                                                                                                                          >  2:1940634/1‑150 (MQ=255)

                           ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaag                                                                                                                          >  1:5055007/1‑150 (MQ=255)

                           ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaag                                                                                                                          >  1:2233080/1‑150 (MQ=255)

                           ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaag                                                                                                                          >  1:1258829/1‑150 (MQ=255)

                           ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaag                                                                                                                          >  2:4325043/1‑150 (MQ=255)

                           ttCACCCCGAAGCACACCCAGACCGCGGATCTGAAAGCCGGCGAAGTAGGCTTCCTCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGGCGAGGTGTTGCCAGGCTGCaag                                                                                                                          >  1:4573192/1‑150 (MQ=255)

                              aCCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaagaag                                                                                                                       >  2:1543669/1‑150 (MQ=255)

                              aCCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaagaag                                                                                                                       >  2:310038/1‑150 (MQ=255)

                               ccccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAagaaga                                                                                                                      <  2:27844/150‑1 (MQ=255)

                                cccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGAt                                                                                                                     >  2:2005812/1‑150 (MQ=255)

                                cccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGAt                                                                                                                     >  1:3270232/1‑150 (MQ=255)

                                cccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGAt                                                                                                                     >  2:4424221/1‑150 (MQ=255)

                                cccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGAt                                                                                                                     >  2:449418/1‑150 (MQ=255)

                                   gAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCa                                                                                                                  <  1:5520702/150‑1 (MQ=255)

                                        acacCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGACCTCGACCTCAGGGGTCGAggacaggg                                                                                                                                     >  1:1680663/1‑120 (MQ=255)

                                          acCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAg                                                                                                           <  1:5172404/150‑1 (MQ=255)

                                            ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGt                                                                                                         >  2:3416128/1‑150 (MQ=255)

                                            ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGt                                                                                                         <  1:3391994/150‑1 (MQ=255)

                                            ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGt                                                                                                         >  2:2597905/1‑150 (MQ=255)

                                            ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGt                                                                                                         >  1:465588/1‑150 (MQ=255)

                                               gACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTa                                                                                                      >  2:4997202/1‑150 (MQ=255)

                                                      gATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGc                                                                                               <  2:1481548/150‑1 (MQ=255)

                                                      gATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGc                                                                                               <  1:3075091/150‑1 (MQ=255)

                                                        tCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCt                                                                                             <  2:883815/150‑1 (MQ=255)

                                                        tCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCATGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCt                                                                                             <  1:1153166/150‑1 (MQ=255)

                                                          tGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGt                                                                                           >  1:3529802/1‑150 (MQ=255)

                                                          tGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGt                                                                                           >  2:764134/1‑150 (MQ=255)

                                                            aaaGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTc                                                                                         >  1:4459894/1‑150 (MQ=255)

                                                             aaGCCGGCGAAGTGGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTcc                                                                                        >  2:1173641/1‑150 (MQ=255)

                                                             aaGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTcc                                                                                        >  2:5512620/1‑150 (MQ=255)

                                                             aaGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTcc                                                                                        <  2:5409077/150‑1 (MQ=255)

                                                             aaGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTcc                                                                                        <  2:1360566/150‑1 (MQ=255)

                                                              aGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTccc                                                                                       <  2:358908/150‑1 (MQ=255)

                                                               gCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCCGTTCCCg                                                                                      <  2:1691146/150‑1 (MQ=255)

                                                                  ggCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTc                                                                                   <  2:5118816/150‑1 (MQ=255)

                                                                   gCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCa                                                                                  <  1:1410142/150‑1 (MQ=255)

                                                                           ggCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCcgac                                                                          <  1:2442280/150‑1 (MQ=255)

                                                                              ttCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCgacgac                                                                       <  2:1863857/150‑1 (MQ=255)

                                                                                  tcatcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACACTGACCCTGACCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTACAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGCCTTCg                                                                   >  1:5085367/1‑150 (MQ=255)

                                                                                   catcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                  <  1:5593476/150‑1 (MQ=255)

                                                                                    atcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAg                                                                 <  2:3453923/150‑1 (MQ=255)

                                                                                         cAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGAttt                                                            >  2:3619259/1‑150 (MQ=255)

                                                                                          agaaTCAAGGACATCCACGGAGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTACCCGGTCAGCTCCGACGACTTCGAGGATTTc                                                           >  1:3957287/4‑150 (MQ=255)

                                                                                          aGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTc                                                           >  2:3908311/1‑150 (MQ=255)

                                                                                              tCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCcgcg                                                       >  1:2985186/1‑150 (MQ=255)

                                                                                               cAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGa                                                      >  2:4603787/1‑150 (MQ=255)

                                                                                               cAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGa                                                      <  2:386223/150‑1 (MQ=255)

                                                                                                aaGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGAc                                                     >  2:17679/1‑150 (MQ=255)

                                                                                                      aTCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTg                                               <  2:4514254/150‑1 (MQ=255)

                                                                                                      aTCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTg                                               >  1:2348313/1‑150 (MQ=255)

                                                                                                       tCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGc                                              >  1:132253/1‑150 (MQ=255)

                                                                                                       tCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGc                                              >  1:1973648/1‑150 (MQ=255)

                                                                                                       tCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGc                                              >  1:857911/1‑150 (MQ=255)

                                                                                                         cACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAg                                            >  2:2356365/1‑150 (MQ=255)

                                                                                                         cACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAg                                            >  1:5346742/1‑150 (MQ=255)

                                                                                                         cACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAg                                            >  1:1517598/1‑150 (MQ=255)

                                                                                                                ccccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGa                                     >  1:1858571/1‑150 (MQ=255)

                                                                                                                 cccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGAc                                    >  1:1233766/1‑150 (MQ=255)

                                                                                                                       tGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGaa                              <  1:318352/150‑1 (MQ=255)

                                                                                                                        ggTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAAc                             <  2:4744549/150‑1 (MQ=255)

                                                                                                                          tGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGa                           >  1:3903821/1‑150 (MQ=255)

                                                                                                                              accctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActcc                       >  1:4915897/1‑150 (MQ=255)

                                                                                                                                cctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActcctc                     >  1:5342309/1‑150 (MQ=255)

                                                                                                                                   gaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCt                  <  2:531155/150‑1 (MQ=255)

                                                                                                                                    accctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTg                 <  1:3672171/150‑1 (MQ=255)

                                                                                                                                      cctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCa               >  2:3723228/1‑150 (MQ=255)

                                                                                                                                      cctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCa               >  1:1528682/1‑150 (MQ=255)

                                                                                                                                      cctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCa               >  1:4346690/1‑150 (MQ=255)

                                                                                                                                      cctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCa               >  1:5244124/1‑150 (MQ=255)

                                                                                                                                       aacaccTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAg              >  2:2899587/5‑150 (MQ=255)

                                                                                                                                         gTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTa            <  1:3723261/150‑1 (MQ=255)

                                                                                                                                         gTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTa            <  1:4023987/150‑1 (MQ=255)

                                                                                                                                          accTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTAc           >  1:2156396/2‑150 (MQ=255)

                                                                                                                                           ccTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACa          <  2:4788640/150‑1 (MQ=255)

                                                                                                                                            cTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACAt         >  2:570506/1‑150 (MQ=255)

                                                                                                                                             tCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATg        >  1:358311/1‑150 (MQ=255)

                                                                                                                                               gACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACTCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATgac      >  2:3665797/1‑148 (MQ=255)

                                                                                                                                               gACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTAGATGcc      >  1:1187468/1‑150 (MQ=255)

                                                                                                                                                aCCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCg     >  1:4999094/1‑150 (MQ=255)

                                                                                                                                                  cccTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGa   >  2:2675870/1‑150 (MQ=255)

                                                                                                                                                   ccTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGaa  <  2:5055007/150‑1 (MQ=255)

                                                                                                                                                   ccTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGaa  <  2:2226206/150‑1 (MQ=255)

GTGCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAA  >  NC_002947/1636937‑1637233

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