BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,637,088	A→G	6.6%	E284G (GAG→GGG)	lepA →	elongation factor 4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,637,088	0	A	G	6.6%	765.0 / 27.2	243	E284G (GAG→GGG)	lepA	elongation factor 4
Reads supporting (aligned to +/- strand): ref base A (141/85); new base G (7/9); total (149/94)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 247 total aligned reads displayed.

GCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCT  >  NC_002947/1636939‑1637237

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

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

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

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

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

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

          gACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGgtg                                                                                                                                             >  1:5300810/1‑150 (MQ=255)

          gACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGgtg                                                                                                                                             >  1:116376/1‑150 (MQ=255)

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

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

            cAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCCTCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGACCTCGACCTCAGGGGTCGAggacag                                                                                                                                           <  2:3601252/150‑5 (MQ=255)

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

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

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

                gTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCa                                                                                                                                       <  1:1229208/150‑1 (MQ=255)

                 tCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAg                                                                                                                                      >  2:318352/1‑150 (MQ=255)

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

                        cTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaa                                                                                                                               <  2:2609883/150‑1 (MQ=255)

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

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

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

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

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

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

                          tCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaaga                                                                                                                             >  1:3725198/1‑150 (MQ=255)

                          tCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaaga                                                                                                                             <  1:3924248/150‑1 (MQ=255)

                           cACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaagaa                                                                                                                            <  1:3354544/150‑1 (MQ=255)

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

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

                              cccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGAt                                                                                                                         <  2:1711050/150‑1 (MQ=255)

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

                                    gCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGcc                                                                                                                   >  1:1245602/1‑150 (MQ=255)

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

                                      acacCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGACCTCGACCTCAGGGGTCGAggacaggg                                                                                                                                         <  2:1680663/126‑7 (MQ=255)

                                       ccaggacaCCGCTGATCTGAAAGCCGGCGAAGTAGGCGTCATCATCGCCAGCCTCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGCGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCa                                                                                                                >  2:550272/8‑150 (MQ=255)

                                       cacCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCa                                                                                                                >  1:4296498/1‑150 (MQ=255)

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

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

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

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

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

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

                                                     aTCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGcc                                                                                                  >  1:3168184/1‑150 (MQ=255)

                                                     aTCTGAAAGCCGGCGAAGTAGGCGTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGcc                                                                                                  <  2:284434/150‑1 (MQ=255)

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

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

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

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

                                                          aaaGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTc                                                                                             <  2:447019/150‑1 (MQ=255)

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

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

                                                            aGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTccc                                                                                           >  1:2925732/1‑150 (MQ=255)

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

                                                                   gAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGc                                                                                    >  2:3128823/1‑150 (MQ=255)

                                                                    caggaGGCTTCATCATCGCCAGCCTCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGAGTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCt                                                                                   <  1:5169827/146‑1 (MQ=255)

                                                                      gTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTcc                                                                                 <  1:4016210/150‑1 (MQ=255)

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

                                                                         ggCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCcgac                                                                              >  1:355114/1‑150 (MQ=255)

                                                                         ggCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCcgac                                                                              >  2:2182176/1‑150 (MQ=255)

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

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

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

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

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

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

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

                                                                                             cAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTCTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGa                                                          <  2:386374/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                       cACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAg                                                <  2:1435414/150‑1 (MQ=255)

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

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

                                                                                                              ccccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAGGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGa                                         >  1:1858485/1‑150 (MQ=255)

                                                                                                                ccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGAcc                                       >  1:4483951/1‑150 (MQ=255)

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

                                                                                                                          acacCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGAct                             >  2:398527/1‑150 (MQ=255)

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

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

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

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

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

                                                                                                                                    cctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCa                   <  1:5060906/150‑1 (MQ=255)

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

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

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

                                                                                                                                     aacaccTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATCTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAg                  >  1:5008110/5‑150 (MQ=255)

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

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

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

                                                                                                                                          cTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACAt             >  2:570506/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)

                                                                                                                                                 ccTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGaa      >  1:1368503/1‑150 (MQ=255)

                                                                                                                                                   tGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAg    >  1:3466499/1‑150 (MQ=255)

                                                                                                                                                   tGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAg    >  2:1837038/1‑150 (MQ=255)

GCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGCT  >  NC_002947/1636939‑1637237

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