BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,613,990	C→G	5.7%	R337R (CGC→CGG) ‡	PP_4931 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,613,990	0	C	G	5.7%	649.9 / 21.0	194	R337R (CGC→CGG) ‡	PP_4931	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (102/81); new base G (8/3); total (110/84)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 203 total aligned reads displayed.

ATTGCCGCAGCGCAGAAGGAAATCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCTGGGCGTACCGGCCTGGGA  >  NC_002947/5613841‑5614139

aTTGCCGCAGCGCAGAAGGAAATCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAgcgc                                                                                                                                                       <  2:376467/150‑1 (MQ=255)

 ttGTCGCAGCGCAGAAGGAAATCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGcc                                                                                                                                                      <  1:2386164/150‑1 (MQ=255)

    ccGCAGCGCAGAAGGAAATCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCctgc                                                                                                                                                   <  1:2786464/150‑1 (MQ=255)

        agcgcagAAGGAAATCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCtggt                                                                                                                                               >  2:4942661/1‑150 (MQ=255)

            cagAAGGAAATCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGc                                                                                                                                           >  1:2150773/1‑150 (MQ=255)

             agAAGGAAATCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGca                                                                                                                                          <  1:3322158/150‑2 (MQ=255)

              gAAGGAAATCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTg                                                                                                                                         >  2:3326945/1‑150 (MQ=255)

              gAAGGAAATCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTg                                                                                                                                         >  2:2867141/1‑150 (MQ=255)

               aaGGAAATCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTgg                                                                                                                                        <  1:2718835/150‑1 (MQ=255)

                 ggAAATCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGcc                                                                                                                                      >  1:1760428/1‑150 (MQ=255)

                 ggAAATCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGcc                                                                                                                                      >  1:2896146/1‑150 (MQ=255)

                   aaaTCGACAACCTGCTGCCATGGGTCGTTCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGGGCCTGCACAATCCGGCCTGGCACGCCCGGACAGCGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGa                                                                                                                                    >  2:2366398/1‑150 (MQ=255)

                   aaaTCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGa                                                                                                                                    <  1:5179539/150‑1 (MQ=255)

                    aaTCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGAc                                                                                                                                   <  2:3687592/150‑1 (MQ=255)

                        gACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCaaa                                                                                                                               <  1:96873/150‑1 (MQ=255)

                         aCAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAac                                                                                                                              >  2:1119073/1‑150 (MQ=255)

                          cAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAact                                                                                                                             <  2:1224124/150‑1 (MQ=255)

                               tgctgcCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGGCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAactggcac                                                                                                                        >  2:1929154/1‑150 (MQ=255)

                                  tgcCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAactggcactgg                                                                                                                     <  1:4813156/150‑1 (MQ=255)

                                  tgcCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAactggcactgg                                                                                                                     <  2:3248772/150‑1 (MQ=255)

                                  tgcCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTAGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAactggcactgg                                                                                                                     >  2:819479/1‑150 (MQ=255)

                                    ccATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAActggcactggca                                                                                                                   >  1:132404/1‑150 (MQ=255)

                                     cATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACtggcactggcac                                                                                                                  <  1:1645173/150‑1 (MQ=255)

                                      aTGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCAcc                                                                                                                 >  1:4609027/1‑150 (MQ=255)

                                         ggTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGa                                                                                                              >  2:1988905/1‑150 (MQ=255)

                                              aTCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCa                                                                                                         <  2:1733750/150‑1 (MQ=255)

                                              aTCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCa                                                                                                         <  1:4860417/150‑1 (MQ=255)

                                              aTCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCa                                                                                                         <  2:2948486/150‑1 (MQ=255)

                                                  gagCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGa                                                                                                     >  2:894726/1‑150 (MQ=255)

                                                      cTGGTGCGCTGGGCGACAGGGCGCGTCGTCCGCGCTGAGCCTGCACATTCCGGCGTGGCAGGCCCGTACTTCGCCTTCATTGCCGACCTGCTGCGCCTGCTGGTGGGCTTGCGGACCTAACTGGCACTGGCACCGGACTTCAGCGACCgg                                                                                                 <  2:2410147/150‑1 (MQ=255)

                                                      cTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCgg                                                                                                 >  1:1736906/1‑150 (MQ=255)

                                                      cTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCgg                                                                                                 <  1:4942661/150‑1 (MQ=255)

                                                       gggTGCGCGGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCggg                                                                                                <  2:4255913/149‑1 (MQ=255)

                                                         gTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCGCAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTa                                                                                              >  1:640819/1‑150 (MQ=255)

                                                          tGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTaa                                                                                             <  1:3670574/150‑1 (MQ=255)

                                                              cTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAg                                                                                         >  2:5168706/1‑150 (MQ=255)

                                                               tGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGc                                                                                        >  2:1756876/1‑150 (MQ=255)

                                                               tGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGc                                                                                        <  1:751390/150‑1 (MQ=255)

                                                                  gCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCAc                                                                                     >  2:5278808/1‑150 (MQ=255)

                                                                    cacaCTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCt                                                                                   >  1:107463/1‑150 (MQ=255)

                                                                    cacaCTGCGCGTCGACCGCGCTGAGCCTACACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCtggtg                                                                                                                                              <  1:1377280/91‑1 (MQ=255)

                                                                     acagTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGCCCAGCCGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTg                                                                                  <  1:742952/150‑1 (MQ=255)

                                                                       acTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAactggca                                                                                                                         >  1:1025376/1‑109 (MQ=255)

                                                                       acTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAactggca                                                                                                                         <  2:1025376/109‑1 (MQ=255)

                                                                          gcgcGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAactggcact                                                                                                                       >  1:302424/1‑108 (MQ=255)

                                                                          gcgcGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAgcg                                                                             >  1:4384096/1‑150 (MQ=255)

                                                                            gcgTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGcgcg                                                                           >  2:2610059/1‑150 (MQ=255)

                                                                            gcgTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGcgcg                                                                           >  2:1841829/1‑150 (MQ=255)

                                                                                  aCCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCa                                                                     >  1:12818/1‑150 (MQ=255)

                                                                                  aCCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCa                                                                     <  2:2292147/150‑1 (MQ=255)

                                                                                    cgcgCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCAtt                                                                   >  2:1652054/1‑150 (MQ=255)

                                                                                    cgcgCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCAtg                                                                   >  1:3273609/1‑149 (MQ=255)

                                                                                        cTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTc                                                               <  1:3936574/150‑1 (MQ=255)

                                                                                         tGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTcc                                                              >  2:4733279/1‑150 (MQ=255)

                                                                                             ccTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACaa                                                          >  2:2915388/1‑150 (MQ=255)

                                                                                             ccTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACaa                                                          >  2:3087395/1‑150 (MQ=255)

                                                                                             ccTCCACACTCCGGCCTGCCACGCCCGGACACCGCCTCCCGGGCCGCCCCGCAGCGCCTGCGGGGGGGCGGGCCCACCAACCTGGCCCTGGCACCGCACGCCAGCGACCGGGTAACCAGCCACCTGGAGCGCGACGGCATTCGTCCACaa                                                          <  1:1988905/150‑1 (MQ=255)

                                                                                               tGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGc                                                        >  1:187645/1‑150 (MQ=255)

                                                                                                  acacaccGCCCTGGCACGCCCGCACCACCCCTGCCGTCCCCCCCAGCCGCGCCTCCGGGGGGGCGGGCCGACCAAACTGCCACTGGCACCGGACTTCAGCGACCGGGCAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGcaca                                                     <  1:2337953/145‑1 (MQ=255)

                                                                                                  acaATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGcaca                                                     >  2:1855979/1‑150 (MQ=255)

                                                                                                  acaATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGcaca                                                     >  2:4031681/1‑150 (MQ=255)

                                                                                                      tCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGcc                                                 >  2:829755/1‑150 (MQ=255)

                                                                                                      tCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGcc                                                 >  1:3156012/1‑150 (MQ=255)

                                                                                                        cGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCa                                                                                                         >  2:84279/1‑92 (MQ=255)

                                                                                                        cGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGa                                               >  1:4778627/1‑150 (MQ=255)

                                                                                                        cGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGa                                               >  2:4997467/1‑150 (MQ=255)

                                                                                                         ggCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGAc                                              <  2:3156012/150‑1 (MQ=255)

                                                                                                             tGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTga                                          <  1:5396909/150‑1 (MQ=255)

                                                                                                              ggCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTgac                                         >  1:2990345/1‑150 (MQ=255)

                                                                                                                 aCGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTgaccga                                      >  2:1895366/1‑150 (MQ=255)

                                                                                                                  aggccGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTTGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTgaccgac                                     >  2:195434/4‑150 (MQ=255)

                                                                                                                     ccGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGa                                                                                                              <  2:4592347/74‑1 (MQ=255)

                                                                                                                     ccGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTg                                  >  2:749181/1‑150 (MQ=255)

                                                                                                                     ccGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTg                                  >  2:3398497/1‑150 (MQ=255)

                                                                                                                        gACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCa                               >  2:4638915/1‑150 (MQ=255)

                                                                                                                           aaCGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGc                            >  1:4692544/1‑150 (MQ=255)

                                                                                                                           aaCGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGc                            >  1:436507/1‑150 (MQ=255)

                                                                                                                           aaCGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACATGCCACGc                            >  2:2101999/1‑150 (MQ=255)

                                                                                                                            aCGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGTACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGcc                           <  1:3411842/150‑1 (MQ=255)

                                                                                                                               ccttccttGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTc                        >  1:4507439/1‑150 (MQ=255)

                                                                                                                                 ttccttGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCg                      >  1:5009810/1‑150 (MQ=255)

                                                                                                                                   ccttGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCt                    >  1:1461469/1‑150 (MQ=255)

                                                                                                                                    cttGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCTg                   >  2:4720266/1‑150 (MQ=255)

                                                                                                                                     ttGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCTgg                  <  1:4635605/150‑1 (MQ=255)

                                                                                                                                      tGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGcc                           <  1:934930/140‑1 (MQ=255)

                                                                                                                                        ccGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCTGGGCg               >  1:4771741/1‑150 (MQ=255)

                                                                                                                                        ccGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCTGGGCg               >  1:179055/1‑150 (MQ=255)

                                                                                                                                         ccaCCAGCAGGCGCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCTGGGCGt              >  2:3033905/3‑150 (MQ=255)

                                                                                                                                         ccaCCAGCAGGCGCTGCTGGTCGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCTGGGCGt              >  2:142751/3‑150 (MQ=255)

                                                                                                                                         cGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGGACCGGGCTTCAGCGACCGGGTAATCTGCCACCTGGAGCGCGTCGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCGGCCACGCCCTCCGCTGGGCGt              >  2:1593520/1‑150 (MQ=255)

                                                                                                                                          gACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCTGGGCGTa             >  1:4182740/1‑150 (MQ=255)

                                                                                                                                          caCCAGCAGGCGCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCTgg                  <  1:1844508/144‑1 (MQ=255)

                                                                                                                                          caCCAGCAGGCGCTGCTGGTCGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTga                                          >  1:1802805/2‑121 (MQ=255)

                                                                                                                                          caCCAGCAGGCGCTGCTGGTCGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTga                                          <  2:1802032/120‑1 (MQ=255)

                                                                                                                                           aCCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCATTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTgac                                         <  2:3617683/121‑1 (MQ=255)

                                                                                                                                           aCCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCATTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTgac                                         >  1:3617683/1‑121 (MQ=255)

                                                                                                                                            ccAGCAGCGCCGGCTGGTGGGCTGGCCGACCCACCTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCAGCGGAGGGCGtggc           <  2:4552726/150‑4 (MQ=255)

                                                                                                                                                cagcGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCTGGGCGTACCGGcc       <  2:1461469/150‑1 (MQ=255)

                                                                                                                                                   cgcCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCTGGGCGTACCGGCCTgg    <  1:1119073/150‑1 (MQ=255)

                                                                                                                                                     ccTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCTGGGCGTACCGGCCTGGGa  <  2:4384096/150‑1 (MQ=255)

ATTGCCGCAGCGCAGAAGGAAATCGACAACCTGCTGCCATGGGTCGATCAGAGCCTGGTGCGCTGGGCCACACTGCGCGTCGACCGCGCTGAGCCTGCACAATCCGGCCTGGCACGCCCGGACAACGCCTTCCTTGCCGACCAGCAGCGCCTGCTGGTGGGCTGGCCGACCAAACTGGCACTGGCACCGGACTTCAGCGACCGGGTAATCAGCCACCTGGAGCGCGACGGCATTCGTCCACAAGCACAGGCCGACCTGACCGACCTGCCACGCCCTCCGCTGGGCGTACCGGCCTGGGA  >  NC_002947/5613841‑5614139

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