BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,624,905	A→C	7.2%	M269R (ATG→AGG) ‡	PP_4941 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,624,905	0	A	C	7.2%	522.9 / 22.1	163	M269R (ATG→AGG) ‡	PP_4941	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (74/77); new base C (6/6); total (80/83)
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.98e-01

Only 100 of 167 total aligned reads displayed.

GCAAGGCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGGTAG  >  NC_002947/5624756‑5625054

 cAAGGCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCAt                                                                                                                                                      >  2:2200484/1‑150 (MQ=255)

  aaGGCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATCTCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATg                                                                                                                                                     >  2:2324311/1‑150 (MQ=255)

    ggCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggc                                                                                                                                                   <  2:4448073/150‑1 (MQ=255)

     gCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggca                                                                                                                                                  <  1:3643846/150‑1 (MQ=255)

       tCGCAGGTCGGCCAGCCGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcagg                                                                                                                                                >  1:2354028/1‑150 (MQ=255)

       tCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcagg                                                                                                                                                >  1:3532322/1‑150 (MQ=255)

       tCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcagg                                                                                                                                                >  2:3399035/1‑150 (MQ=255)

       tCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcagg                                                                                                                                                <  1:1261386/150‑1 (MQ=255)

         gCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcaggca                                                                                                                                              >  1:1552574/1‑150 (MQ=255)

         gCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcaggca                                                                                                                                              <  1:2465201/150‑1 (MQ=255)

         gCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcaggca                                                                                                                                              <  1:2464994/150‑1 (MQ=255)

              tCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGc                                                                                                                                         <  1:4572786/150‑1 (MQ=255)

               cGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCa                                                                                                                                        <  1:654562/150‑1 (MQ=255)

               cGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCa                                                                                                                                        <  1:4806941/150‑1 (MQ=255)

                 gCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCCAtgcaggca                                                                                                                                              <  2:3704657/142‑1 (MQ=255)

                 gCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCCAtgcaggca                                                                                                                                              >  1:3704657/1‑142 (MQ=255)

                   cagcagGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCCATGCAGGCAGca                                                                                                                                           <  1:3439380/143‑2 (MQ=255)

                   cagcagGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCCATGCAGGCAGca                                                                                                                                           >  2:3439380/1‑142 (MQ=255)

                    agcagGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCaa                                                                                                                                   <  1:2870043/150‑1 (MQ=255)

                      cagGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCaaaa                                                                                                                                 <  2:796515/150‑1 (MQ=255)

                       agGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAg                                                                                                                                <  2:789258/150‑1 (MQ=255)

                         gTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTa                                                                                                                              <  1:1838025/150‑1 (MQ=255)

                         gTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTa                                                                                                                              >  2:952511/1‑150 (MQ=255)

                              gCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCa                                                                                                                         >  1:3338238/1‑150 (MQ=255)

                               ccAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAg                                                                                                                        >  2:2001030/1‑150 (MQ=255)

                                cAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGc                                                                                                                       >  1:267079/1‑150 (MQ=255)

                                 aGGTTGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGcg                                                                                                                      >  1:123006/1‑149 (MQ=255)

                                 aGGTTGGTGCGCCATCCCGGTCGCCAGGTAACTGGAAAGGGTTGAACTGGCGGAACACGCGGATAACGGCAAGGTCGGCATCACGAGAGGCGGCCGCCTGGTCCGGTGCAGCCTGCAAGGCAGGCAGCTGCAGGCAAAAGTAGCCCATCt                                                                                                                      >  1:794544/1‑150 (MQ=255)

                                 aGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCt                                                                                                                      >  1:1033994/1‑150 (MQ=255)

                                   gtgggtgCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGt                                                                                                                    >  2:748182/1‑150 (MQ=255)

                                   gtgggtgCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGt                                                                                                                    <  1:2522533/150‑1 (MQ=255)

                                         gcgcCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCaa                                                                                                              <  1:2509255/150‑1 (MQ=255)

                                           gcCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAAcg                                                                                                            >  1:1008897/1‑150 (MQ=255)

                                           gcCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAAcg                                                                                                            >  1:1405809/1‑150 (MQ=255)

                                           gcCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAAcg                                                                                                            >  2:4748516/1‑150 (MQ=255)

                                             cATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAAcgcg                                                                                                          >  1:1438669/1‑150 (MQ=255)

                                              aTCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGa                                                                                                         <  1:1504447/150‑1 (MQ=255)

                                                   ggTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCg                                                                                                    <  1:3789891/150‑1 (MQ=255)

                                                       gCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCAc                                                                                                >  1:4800497/1‑150 (MQ=255)

                                                                 gggaacgggTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCaa                                                                                      >  2:1287699/1‑150 (MQ=255)

                                                                  ggaacgggTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAAc                                                                                     >  2:2749622/1‑150 (MQ=255)

                                                                    aacgggTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAAcgc                                                                                   <  1:4475559/150‑1 (MQ=255)

                                                                     acgggTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAAcgcg                                                                                  <  1:4925874/150‑1 (MQ=255)

                                                                         gTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCa                                                                               >  1:90562/1‑149 (MQ=255)

                                                                         gTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCa                                                                               <  2:90615/149‑1 (MQ=255)

                                                                             aaCTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCACAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCacgacg                                                                          >  2:3173792/1‑150 (MQ=255)

                                                                              aCTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCAcgacga                                                                         <  1:165412/150‑1 (MQ=255)

                                                                                tGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATg                                                                       >  1:1568035/1‑150 (MQ=255)

                                                                                tGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATg                                                                       >  2:3226528/1‑150 (MQ=255)

                                                                                tGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATg                                                                       <  2:2736017/150‑1 (MQ=255)

                                                                                        cacGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCa                                                               <  2:2345295/150‑1 (MQ=255)

                                                                                         acGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAg                                                              <  2:476657/150‑1 (MQ=255)

                                                                                         acGCGGATATCGGCCAGGTCGGCTTCACGAGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAg                                                              <  2:178646/150‑1 (MQ=255)

                                                                                          cgcgGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAgg                                                             <  1:2838330/150‑1 (MQ=255)

                                                                                           gcgGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGc                                                            >  2:1086311/1‑150 (MQ=255)

                                                                                                tatCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCa                                                       <  2:2173330/150‑1 (MQ=255)

                                                                                                tatCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCa                                                       <  2:651988/150‑1 (MQ=255)

                                                                                                   cGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGgc                                                    <  2:2959556/150‑1 (MQ=255)

                                                                                                     gCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGgcgc                                                  >  2:3166498/1‑150 (MQ=255)

                                                                                                     gCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGgcgc                                                  >  1:2884826/1‑150 (MQ=255)

                                                                                                       cAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGccc                                                <  1:4750899/150‑1 (MQ=255)

                                                                                                          gTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCtt                                             <  2:674504/150‑1 (MQ=255)

                                                                                                            cGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCa                                           <  1:2274797/150‑1 (MQ=255)

                                                                                                             ggCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAg                                          <  1:3208252/150‑1 (MQ=255)

                                                                                                             ggCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAg                                          <  1:2946720/150‑1 (MQ=255)

                                                                                                                 tCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgcgc                                      >  2:4911970/1‑150 (MQ=255)

                                                                                                                    cGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCGCCTGAAGgggcgcc                                   <  1:1829336/150‑2 (MQ=255)

                                                                                                                     gTGAGGCGGCCGGCTGGTCGGGTACTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcgg                                  >  1:3136243/1‑150 (MQ=255)

                                                                                                                       gAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcc                                <  2:3225452/150‑1 (MQ=255)

                                                                                                                          gcggcCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcctcg                             <  2:267079/150‑1 (MQ=255)

                                                                                                                          gcggcCGGCTGGTCGGGTACTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcctcg                             <  2:3136243/150‑1 (MQ=255)

                                                                                                                            ggccggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcctcggc                           <  1:2511859/150‑1 (MQ=255)

                                                                                                                                ggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCATGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCttt                       <  1:4190574/150‑1 (MQ=255)

                                                                                                                                ggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCttt                       <  2:2651860/150‑1 (MQ=255)

                                                                                                                                ggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCttt                       <  2:1552574/150‑1 (MQ=255)

                                                                                                                                  cTGGTCGGGTGCTGCCTGCCTGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGc                     >  1:2124768/1‑150 (MQ=255)

                                                                                                                                    ggTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCa                   <  2:2574496/150‑1 (MQ=255)

                                                                                                                                    ggTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCa                   >  2:628711/1‑150 (MQ=255)

                                                                                                                                     gTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAg                  <  2:4963065/150‑1 (MQ=255)

                                                                                                                                          gTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCa             >  1:2227357/1‑150 (MQ=255)

                                                                                                                                             gcagctgcctgcCATGCAGCTGCAGGCAAAAGTAGCCCAGCAGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCGTGGCCAGGCCCATct          >  1:3867203/5‑150 (MQ=255)

                                                                                                                                             ctgcctgcCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCtt                        >  2:4750255/1‑136 (MQ=255)

                                                                                                                                             ctgcctgcCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCtt                        <  1:4750255/136‑1 (MQ=255)

                                                                                                                                             ctgcctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATct          >  1:1099114/1‑150 (MQ=255)

                                                                                                                                              tgcctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATctg         >  1:3810198/1‑150 (MQ=255)

                                                                                                                                               gcctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATctgc        >  1:2247215/1‑150 (MQ=255)

                                                                                                                                                cctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATctgct       <  1:2200484/150‑1 (MQ=255)

                                                                                                                                                cctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATctgct       <  1:2169159/150‑1 (MQ=255)

                                                                                                                                                 ctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGAt                                                                        <  1:2622802/84‑1 (MQ=255)

                                                                                                                                                 ctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGAt                                                                        >  2:2622802/1‑84 (MQ=255)

                                                                                                                                                  tgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATTCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTgg     <  1:1620150/150‑1 (MQ=255)

                                                                                                                                                  tgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTgg     >  2:4834142/1‑150 (MQ=255)

                                                                                                                                                  tgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTgg     <  2:3338238/150‑1 (MQ=255)

                                                                                                                                                  tgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTgg     >  1:4061271/1‑150 (MQ=255)

                                                                                                                                                  tgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTgg     >  1:2869470/1‑150 (MQ=255)

                                                                                                                                                   gcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGGt    >  1:3027275/1‑150 (MQ=255)

                                                                                                                                                   gcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGGt    >  1:472425/1‑150 (MQ=255)

                                                                                                                                                    cATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGGTa   <  2:4478527/150‑1 (MQ=255)

                                                                                                                                                     aTGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGGTAg  >  1:3156333/1‑150 (MQ=255)

                                                                                                                                                     aTGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGGTAg  >  1:2663665/1‑150 (MQ=255)

GCAAGGCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGGTAG  >  NC_002947/5624756‑5625054

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