BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,378,134	G→A	5.1%	R467R (CGC→CGT)	dnaK ←	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,378,134	0	G	A	5.1%	622.6 / 14.2	197	R467R (CGC→CGT)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand): ref base G (109/78); new base A (6/4); total (115/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 209 total aligned reads displayed.

GCACCATACGCTCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAT  >  NC_002947/5377986‑5378283

  aCCATACGCTCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCcgcgc                                                                                                                                                    >  2:130684/1‑150 (MQ=255)

         gCTCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGTATGCCGTTGGCGTCGATCTCGAAGGTCACATCGATTTGCGGCACGCCGCGCGGTGCTg                                                                                                                                             >  1:1591362/1‑150 (MQ=255)

         gCTCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTg                                                                                                                                             <  1:417635/150‑1 (MQ=255)

         gCTCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTg                                                                                                                                             >  2:559191/1‑150 (MQ=255)

           tCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGc                                                                                                                                           >  1:3935071/1‑150 (MQ=255)

                  ctTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGt                                                                                                                                    >  1:2693374/1‑150 (MQ=255)

                          gACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAgg                                                                                                                            >  2:383848/1‑150 (MQ=255)

                            cAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTc                                                                                                                          <  1:1511875/150‑1 (MQ=255)

                             aGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCg                                                                                                                         <  2:1255278/150‑1 (MQ=255)

                               gCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGaa                                                                                                                       >  2:5030180/1‑150 (MQ=255)

                                ccGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAAc                                                                                                                      >  1:1458800/1‑150 (MQ=255)

                                 cGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACt                                                                                                                     >  2:3031742/1‑150 (MQ=255)

                                   gAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTg                                                                                                                   <  2:2862283/150‑1 (MQ=255)

                                      ttGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGccc                                                                                                                >  1:636455/1‑150 (MQ=255)

                                       tGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCa                                                                                                               <  2:1756644/150‑1 (MQ=255)

                                           cTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGa                                                                                                           >  1:2329558/1‑150 (MQ=255)

                                            ttGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCcgcgc                                                                                                                                                    <  1:3170309/108‑1 (MQ=255)

                                            ttGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCcgcgc                                                                                                                                                    >  2:3170309/1‑108 (MQ=255)

                                             tGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACt                                                                                                         >  1:2134212/1‑150 (MQ=255)

                                              gATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACtt                                                                                                        <  1:1999422/150‑1 (MQ=255)

                                               aTCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTg                                                                                                       <  1:1340379/150‑1 (MQ=255)

                                               aTCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTg                                                                                                       <  2:1519438/150‑1 (MQ=255)

                                               aTCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTg                                                                                                       <  2:623984/150‑1 (MQ=255)

                                               aTCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTg                                                                                                       <  1:2465300/150‑1 (MQ=255)

                                                 cACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGtt                                                                                                     >  1:951746/1‑150 (MQ=255)

                                                    gATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTg                                                                                                  >  2:3924795/1‑150 (MQ=255)

                                                     ctCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTgc                                                                                                 <  1:3148498/149‑1 (MQ=255)

                                                      tCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTgcg                                                                                                >  1:195408/1‑150 (MQ=255)

                                                      tCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTgcg                                                                                                <  1:4081090/150‑1 (MQ=255)

                                                           tGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGcc                                                                                           >  2:3675436/1‑150 (MQ=255)

                                                           tGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGcc                                                                                           >  1:4053634/1‑150 (MQ=255)

                                                           tGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGcc                                                                                           >  1:2767857/1‑150 (MQ=255)

                                                               tCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCt                                                                                       >  2:4078807/1‑150 (MQ=255)

                                                                 ttGCCTGTGCCCTTGTCCTTTGCGGACACGTTCACGATACCGTTGGCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGccatttca                                                                                     >  2:2424572/1‑144 (MQ=255)

                                                                 ttGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTa                                                                                     >  1:1085434/1‑150 (MQ=255)

                                                                 ttGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCATGCTTa                                                                                     >  1:5072860/1‑150 (MQ=255)

                                                                 ttGCCGGTGGCCTTGCCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTa                                                                                     >  2:4462899/1‑150 (MQ=255)

                                                                   gCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACg                                                                                   <  1:4957218/150‑1 (MQ=255)

                                                                          gCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCt                                                                            >  2:2441368/1‑150 (MQ=255)

                                                                          gCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCt                                                                            <  1:383848/150‑1 (MQ=255)

                                                                            cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCa                                                                                                                                  <  1:4697699/94‑1 (MQ=255)

                                                                            cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCa                                                                                                                                  >  2:4697699/1‑94 (MQ=255)

                                                                            cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTg                                                                          >  1:100934/1‑150 (MQ=255)

                                                                            cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTg                                                                          >  1:4275496/1‑150 (MQ=255)

                                                                             ttGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgc                                                                         >  2:3584693/1‑150 (MQ=255)

                                                                             ttGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgc                                                                         >  1:98162/1‑150 (MQ=255)

                                                                              tGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                        <  1:2956362/150‑1 (MQ=255)

                                                                              tGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                        <  2:4981605/150‑1 (MQ=255)

                                                                                tCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGTTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgcagc                                                                      >  1:1062923/1‑150 (MQ=255)

                                                                                tCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgcagc                                                                      >  1:442630/1‑150 (MQ=255)

                                                                                  tttCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCAc                                                                    >  1:4869438/1‑150 (MQ=255)

                                                                                  tttCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCAc                                                                    >  1:2843011/1‑150 (MQ=255)

                                                                                  tttCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCAc                                                                    >  2:3444749/1‑150 (MQ=255)

                                                                                          cTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAAt                                                                                                                                      >  2:4390566/1‑76 (MQ=255)

                                                                                             cGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCa                                                         <  1:1836426/150‑1 (MQ=255)

                                                                                              gTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCAc                                                        >  1:3796675/1‑150 (MQ=255)

                                                                                               tGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACg                                                       >  1:4854576/1‑150 (MQ=255)

                                                                                               tGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACg                                                       <  2:1881961/150‑1 (MQ=255)

                                                                                                gCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGACTGGTCACgg                                                      >  2:4199764/1‑150 (MQ=255)

                                                                                                gCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACgg                                                      >  1:676058/1‑150 (MQ=255)

                                                                                                 cAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGgc                                                     >  2:5290270/1‑150 (MQ=255)

                                                                                                  acgATACCGTTGGCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGgcc                                                    >  1:4335881/3‑149 (MQ=255)

                                                                                                   cgATACCGTTGGCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGgccg                                                   >  2:974084/2‑148 (MQ=255)

                                                                                                    gATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGct                                                  >  2:2199384/1‑150 (MQ=255)

                                                                                                     aTGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGaa                                                                                                                       <  1:4117737/80‑1 (MQ=255)

                                                                                                     aTGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGaa                                                                                                                       >  2:4117737/1‑80 (MQ=255)

                                                                                                       gCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTgc                                                                                                 <  1:5303098/100‑1 (MQ=255)

                                                                                                       gCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTg                                               <  1:4984775/150‑1 (MQ=255)

                                                                                                         cGTCGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGt                                             <  1:3422753/150‑1 (MQ=255)

                                                                                                          gTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGtt                                            <  2:5406590/150‑1 (MQ=255)

                                                                                                             ggCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTgtc                                         <  2:399109/150‑1 (MQ=255)

                                                                                                             ggCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTgtc                                         <  2:3611848/150‑1 (MQ=255)

                                                                                                              gCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTgtcg                                        >  1:4152048/1‑150 (MQ=255)

                                                                                                              gCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTgtcg                                        >  1:754823/1‑150 (MQ=255)

                                                                                                                 tCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTgtcttca                                     <  2:1653132/150‑2 (MQ=255)

                                                                                                                       tCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                               <  1:3534671/150‑1 (MQ=255)

                                                                                                                       tCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                               <  2:2693374/150‑1 (MQ=255)

                                                                                                                       tCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                               <  2:609203/150‑1 (MQ=255)

                                                                                                                        cGAAGGTCACTTCGATGTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCg                              <  2:4352931/150‑1 (MQ=255)

                                                                                                                         gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                             >  2:4431275/1‑150 (MQ=255)

                                                                                                                         gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                             >  1:3326904/1‑150 (MQ=255)

                                                                                                                         gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                             >  2:1192355/1‑150 (MQ=255)

                                                                                                                         gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                             >  2:1031645/1‑150 (MQ=255)

                                                                                                                              tCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAaca                        >  1:5370355/1‑150 (MQ=255)

                                                                                                                              tCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAaca                        >  1:4489858/1‑150 (MQ=255)

                                                                                                                               cACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAacac                       >  1:4679099/1‑150 (MQ=255)

                                                                                                                                 cTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCt                     <  2:5190282/150‑1 (MQ=255)

                                                                                                                                  ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCTGGTCGAACTTGCCCAGCGAATAGTTCTCCGCCGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                    >  2:2763355/1‑150 (MQ=255)

                                                                                                                                  ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGACCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACATTGCAGCACGTGAATGGTCACGGGGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                    >  2:3672717/1‑150 (MQ=255)

                                                                                                                                   tCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGc                   <  1:1031645/150‑1 (MQ=255)

                                                                                                                                     gATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGa                 <  2:5249666/150‑1 (MQ=255)

                                                                                                                                      aTTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGAc                <  2:2134212/150‑1 (MQ=255)

                                                                                                                                         tGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGActtc             <  2:514602/150‑1 (MQ=255)

                                                                                                                                            ggCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTg          <  2:2245278/150‑1 (MQ=255)

                                                                                                                                            ggCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTg          >  1:398183/1‑150 (MQ=255)

                                                                                                                                             gCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                               <  2:11832/128‑1 (MQ=255)

                                                                                                                                                cgccgcGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCg      <  2:4770822/150‑1 (MQ=255)

                                                                                                                                                   cgcgcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGGCTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGa   <  2:2959103/150‑1 (MQ=255)

                                                                                                                                                   cgcgcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGa   >  2:1013088/1‑150 (MQ=255)

                                                                                                                                                    gcgcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAt  <  2:4053634/150‑1 (MQ=255)

GCACCATACGCTCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAT  >  NC_002947/5377986‑5378283

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