BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,351,634	T→A	5.8%	L20H (CTC→CAC) ‡	PP_4703 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,351,634	0	T	A	5.8%	1086.1 / 35.5	324	L20H (CTC→CAC) ‡	PP_4703	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (146/158); new base A (14/5); total (161/163)
Fisher's exact test for biased strand distribution p-value = 3.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 327 total aligned reads displayed.

CGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCCCCTGGCCCTGG  >  NC_002947/5351486‑5351783

cGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTc                                                                                                                                                      >  2:3691166/1‑150 (MQ=255)

 gACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCg                                                                                                                                                     <  1:1136808/150‑1 (MQ=255)

 gACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCg                                                                                                                                                     <  2:3910354/150‑1 (MQ=255)

   cccGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGc                                                                                                                                                   >  1:5287450/1‑150 (MQ=255)

    ccGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGcc                                                                                                                                                  >  2:1348750/1‑150 (MQ=255)

    ccGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGcc                                                                                                                                                  >  1:3112758/1‑150 (MQ=255)

       aGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGc                                                                                                                                               <  1:3597700/150‑1 (MQ=255)

       aGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGc                                                                                                                                               <  2:5181017/150‑1 (MQ=255)

         cGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTc                                                                                                                                             >  2:2279829/1‑150 (MQ=255)

          gTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCGAGGGCCAGCCCg                                                                                                                                            <  1:1291829/150‑1 (MQ=255)

           tGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGa                                                                                                                                           <  1:2762662/150‑1 (MQ=255)

            ggTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGAc                                                                                                                                          <  1:74680/150‑1 (MQ=255)

                cAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTgg                                                                                                                                      >  2:2482400/1‑150 (MQ=255)

                 aGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGa                                                                                                                                     >  1:693659/1‑150 (MQ=255)

                        gATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGccc                                                                                                                              <  2:1874038/150‑1 (MQ=255)

                            gACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAgg                                                                                                                          <  2:2009241/150‑1 (MQ=255)

                              cAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGc                                                                                                                        <  2:2214958/150‑1 (MQ=255)

                              cAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGc                                                                                                                        <  2:1849219/150‑1 (MQ=255)

                               aaaCGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCg                                                                                                                       >  2:544229/1‑150 (MQ=255)

                                 aCGCCTCAACCAGTAACCGGCGGCCGCACATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGAc                                                                                                                     >  2:4293486/1‑150 (MQ=255)

                                   gCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGACGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGAccc                                                                                                                   <  1:5516186/150‑1 (MQ=255)

                                   gCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCTTGGGCGattc                                                                                                                   <  2:3004670/150‑4 (MQ=255)

                                     cTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCg                                                                                                                 <  2:5866767/150‑1 (MQ=255)

                                      tCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCgg                                                                                                                >  2:5947894/1‑150 (MQ=255)

                                       cAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGc                                                                                                               <  1:4387174/150‑1 (MQ=255)

                                        aaCCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGcc                                                                                                              >  1:39315/1‑150 (MQ=255)

                                        aaCCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGcc                                                                                                              >  2:5319771/1‑150 (MQ=255)

                                           cAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGt                                                                                                           >  1:5848442/1‑150 (MQ=255)

                                             gTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAGTCGCCCTTGGGCGattccag                                                                                                                >  1:2307570/1‑137 (MQ=255)

                                                   ggcggcCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCg                                                                                                   >  2:6199577/1‑150 (MQ=255)

                                                        ccGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAg                                                                                              >  2:4649082/1‑150 (MQ=255)

                                                          gCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGaa                                                                                                                                    >  2:2839753/1‑110 (MQ=255)

                                                          gCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTa                                                                                            >  2:4610144/1‑150 (MQ=255)

                                                          gCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTa                                                                                            >  2:5002942/1‑150 (MQ=255)

                                                           cAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTAt                                                                                           <  1:4614254/150‑1 (MQ=255)

                                                           cAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTAt                                                                                           <  1:4870524/150‑1 (MQ=255)

                                                           cAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTAt                                                                                           >  2:1432262/1‑150 (MQ=255)

                                                            aaaTCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTAtt                                                                                          >  1:1691595/1‑150 (MQ=255)

                                                             aaTCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTa                                                                                         <  2:3600854/150‑1 (MQ=255)

                                                               tCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGCGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGc                                                                                       <  1:1925218/150‑1 (MQ=255)

                                                               tCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGc                                                                                       >  2:4170880/1‑150 (MQ=255)

                                                               tCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGc                                                                                       >  1:2216427/1‑150 (MQ=255)

                                                               tCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGc                                                                                       >  2:3965415/1‑150 (MQ=255)

                                                                cTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGca                                                                                                 <  1:4484359/139‑1 (MQ=255)

                                                                 tGGCCGCCCTGTACAGCAGGTCGCCATGGAAGTCCTTCGTCGCCGTATCGAAACCCCGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTg                                                                                     <  1:969089/150‑1 (MQ=255)

                                                                 tGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTg                                                                                     <  2:1302251/150‑1 (MQ=255)

                                                                  ggCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTgg                                                                                    <  1:3756067/150‑1 (MQ=255)

                                                                  ggCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTgg                                                                                    <  1:334452/150‑1 (MQ=255)

                                                                  ggCCGCCCTGGACCGCAGGTCGCCATGGAAGCCCTTCGTCGCCGGATCGAACCCCCGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTgg                                                                                    <  2:3809685/150‑1 (MQ=255)

                                                                   gccgccCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGgc                                                                                   <  1:4211936/150‑1 (MQ=255)

                                                                     cgccCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGgcgc                                                                                 >  2:1706623/1‑150 (MQ=255)

                                                                        ccTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTc                                                                              <  1:89553/150‑1 (MQ=255)

                                                                              cagcagGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGc                                                                        <  2:4436438/150‑1 (MQ=255)

                                                                              cagcagGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGc                                                                        <  1:5397139/150‑1 (MQ=255)

                                                                                gcagGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGa                                                                      <  2:4486011/150‑1 (MQ=255)

                                                                                 cagGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGAc                                                                     <  1:4802297/150‑1 (MQ=255)

                                                                                 cagGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGAc                                                                     <  1:6501987/150‑1 (MQ=255)

                                                                                      cGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTccc                                                                <  2:2629636/150‑1 (MQ=255)

                                                                                         cATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCg                                                                                                                                                     <  2:4887594/62‑1 (MQ=255)

                                                                                         cATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAg                                                             >  2:3505091/1‑150 (MQ=255)

                                                                                              aaGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGc                                                        <  2:1958937/150‑1 (MQ=255)

                                                                                                gCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTg                                                      <  1:5201175/150‑1 (MQ=255)

                                                                                                gCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTg                                                      <  1:4918889/150‑1 (MQ=255)

                                                                                                  ccTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGt                                                    >  1:3564228/1‑150 (MQ=255)

                                                                                                   cTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTc                                                   >  1:199912/1‑150 (MQ=255)

                                                                                                   cTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTc                                                   >  2:6642401/1‑150 (MQ=255)

                                                                                                    ttCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCg                                                  >  1:1610708/1‑150 (MQ=255)

                                                                                                     tcgtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCgg                                                 <  2:3231748/150‑1 (MQ=255)

                                                                                                       gtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTg                                               >  2:2172082/1‑150 (MQ=255)

                                                                                                       gtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTg                                               <  1:6015251/150‑1 (MQ=255)

                                                                                                           ccGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCAt                                           >  2:4976429/1‑150 (MQ=255)

                                                                                                             gTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCa                                         >  1:1779146/1‑150 (MQ=255)

                                                                                                                   aaaCCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTa                                                                                            <  1:6669135/93‑1 (MQ=255)

                                                                                                                   aaaCCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAgcgcgc                                   <  1:1256989/150‑1 (MQ=255)

                                                                                                                   aaaCCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCGGGTCGGTGGCATCAgcgcgc                                   <  1:5668916/150‑1 (MQ=255)

                                                                                                                     aCCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTg                                 >  1:1089816/1‑150 (MQ=255)

                                                                                                                     aCCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTg                                 <  2:760294/150‑1 (MQ=255)

                                                                                                                       ccAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGAt                               >  2:3154958/1‑150 (MQ=255)

                                                                                                                        cAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGAtg                              <  2:5685315/150‑1 (MQ=255)

                                                                                                                        cAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGAtg                              >  1:3531983/1‑150 (MQ=255)

                                                                                                                            tCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGAtgctac                          <  2:1089816/150‑3 (MQ=255)

                                                                                                                              aTGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAg                        >  2:4358763/1‑150 (MQ=255)

                                                                                                                              aTGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAg                        <  2:2753521/150‑1 (MQ=255)

                                                                                                                                  gCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAgttgt                    <  1:4976429/150‑1 (MQ=255)

                                                                                                                                   ccTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCCAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAgttgtt                   <  2:4802292/150‑1 (MQ=255)

                                                                                                                                    cTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCGGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGttgttg                  <  2:1551627/150‑1 (MQ=255)

                                                                                                                                       gtcgagCTGGCCGAGGGCCAGCCCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCAt               >  1:3413557/6‑150 (MQ=255)

                                                                                                                                       aCCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCAt               <  2:1002232/150‑1 (MQ=255)

                                                                                                                                          gggCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATccc            <  1:1706623/150‑1 (MQ=255)

                                                                                                                                           agCTGGCCGAGGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATcccc           >  2:3284463/2‑150 (MQ=255)

                                                                                                                                           agCTGGCCGAGGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATcccc           >  1:2400912/2‑150 (MQ=255)

                                                                                                                                           agCTGGCCGAGGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATcccc           >  1:5119418/2‑150 (MQ=255)

                                                                                                                                              tGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctgg        <  1:3882871/150‑1 (MQ=255)

                                                                                                                                                 cccTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccc     >  1:5790116/1‑150 (MQ=255)

                                                                                                                                                  ccTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccct    <  2:354603/150‑1 (MQ=255)

                                                                                                                                                  ccTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccct    >  1:6172322/1‑150 (MQ=255)

                                                                                                                                                  ccTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccct    >  1:837945/1‑150 (MQ=255)

                                                                                                                                                   cTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccctg   >  1:2463863/1‑150 (MQ=255)

                                                                                                                                                   cTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccctg   <  2:693659/150‑1 (MQ=255)

                                                                                                                                                   cTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccctg   >  2:860178/1‑150 (MQ=255)

CGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCCCCTGGCCCTGG  >  NC_002947/5351486‑5351783

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