BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,599,285	T→C	10.6%	Y175C (TAC→TGC) ‡	PP_4071 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,599,285	0	T	C	10.6%	247.5 / 16.9	94	Y175C (TAC→TGC) ‡	PP_4071	membrane protein
Reads supporting (aligned to +/- strand): ref base T (38/46); new base C (4/6); total (42/52)
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.39e-01

GTTTGGCTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCG  >  NC_002947/4599142‑4599429

gTTTGGCTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGc                                                                                                                                              >  1:311264/1‑148 (MQ=255)

 tttGGCTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGa                                                                                                                                            <  8:203339/149‑1 (MQ=255)

  ttGGCTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAg                                                                                                                                           >  5:50436/1‑149 (MQ=255)

     gCTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAgcgc                                                                                                                                        <  2:311264/149‑1 (MQ=255)

     gCTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAgcg                                                                                                                                         <  1:160895/148‑1 (MQ=255)

     gCTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAgcg                                                                                                                                         <  7:186580/148‑1 (MQ=255)

        tGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTg                                                                                                                                     >  4:421283/1‑149 (MQ=255)

             tCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACAGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTg                                                                                                                                <  5:350673/149‑1 (MQ=255)

               cAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggc                                                                                                                              <  6:346826/149‑1 (MQ=255)

               cAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTgg                                                                                                                               <  6:50436/148‑1 (MQ=255)

                agccGGAACCGAAGACTCAACGTCTGAGGAGTCACTGCCCAGCGACGCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGGCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggcc                                                                                                                             <  8:343117/146‑1 (MQ=255)

                aGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggcc                                                                                                                             >  6:270398/1‑149 (MQ=255)

                aGACGGTACCGAAGACTCAACGACAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTcgcc                                                                                                                             >  8:35316/1‑149 (MQ=255)

                 gACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggccg                                                                                                                            <  3:422620/149‑1 (MQ=255)

                   cGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTcgcctac                                                                                                                          >  3:77446/1‑146 (MQ=255)

                    ggTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTcgcctacc                                                                                                                         <  6:411848/149‑5 (MQ=255)

                     gTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTcgcctaccc                                                                                                                        <  6:168473/149‑6 (MQ=255)

                      tACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGcccc                                                                                                                       >  2:358462/1‑149 (MQ=255)

                        ccGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGt                                                                                                                     >  6:404695/1‑149 (MQ=255)

                         cGAAGACTCAACGTTAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTg                                                                                                                    >  8:83633/1‑149 (MQ=255)

                         cGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTg                                                                                                                    >  1:411301/1‑149 (MQ=255)

                          gAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGCTAAAGGCACCGGGCGGGGCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTgg                                                                                                                   <  4:282044/149‑2 (MQ=37)

                          gAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGc                                                                                                                   <  4:391401/149‑1 (MQ=255)

                           aaGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGcc                                                                                                                  >  6:374378/1‑149 (MQ=255)

                           aaGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTggg                                                                                                                  <  4:77446/149‑3 (MQ=38)

                            aGACTCAACGCCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTgggg                                                                                                                 <  7:35316/149‑4 (MQ=37)

                              aCTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                               >  8:8978/1‑149 (MQ=255)

                              aCTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                               >  8:2671/1‑149 (MQ=255)

                              aCTCAACGTCAGGGGAGGCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgc                                                                                                               <  3:151650/149‑1 (MQ=255)

                               cTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcg                                                                                                              <  2:35839/149‑1 (MQ=255)

                                tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                             >  7:20697/1‑149 (MQ=255)

                                tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                             >  5:201125/1‑149 (MQ=255)

                                tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTAcccc                                                                                                                       >  7:191986/1‑139 (MQ=38)

                                tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTAcccc                                                                                                                       <  8:191986/139‑1 (MQ=38)

                                 cAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGcc                                                                                                            <  7:2671/149‑1 (MQ=255)

                                 cAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGcc                                                                                                            <  7:8978/149‑1 (MQ=255)

                                  aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                             <  6:378962/147‑1 (MQ=255)

                                  aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                           >  5:332642/1‑149 (MQ=255)

                                  aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                           >  7:105538/1‑149 (MQ=255)

                                  aaCGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCt                                                                                                           >  5:52232/1‑149 (MQ=255)

                                            gAGTCACTGCCCAGCGAATCTGTGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCg                                                                                                 >  6:12035/1‑149 (MQ=255)

                                            gAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGc                                                                                                  <  6:81427/148‑1 (MQ=255)

                                               tCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGt                                                                                              >  8:51009/1‑149 (MQ=255)

                                               tCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGt                                                                                              <  5:54930/149‑1 (MQ=255)

                                                   tGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAAc                                                                                          >  7:214131/1‑149 (MQ=255)

                                                     cccAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAAc                                                                                          >  7:377068/1‑147 (MQ=255)

                                                       cAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgag                                                                                       >  7:132207/1‑148 (MQ=255)

                                                       cAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACAAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgag                                                                                       >  7:339325/1‑148 (MQ=255)

                                                        aGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                      >  5:122338/1‑148 (MQ=255)

                                                         gCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGaaa                                                                                    <  6:122338/149‑1 (MQ=255)

                                                             aTCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACttt                                                                                <  2:411301/149‑1 (MQ=255)

                                                                  cGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTg                                                                               <  5:441286/145‑1 (MQ=255)

                                                                  cGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATc                                                                            <  6:201125/148‑1 (MQ=255)

                                                                     ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                        >  2:69170/1‑149 (MQ=255)

                                                                     ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                        <  5:12035/149‑1 (MQ=255)

                                                                     ggCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                        <  6:52232/149‑1 (MQ=255)

                                                                      gCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTg                                                                        >  8:31085/1‑148 (MQ=255)

                                                                       ccGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCa                                                                      <  6:332642/149‑1 (MQ=255)

                                                                        cGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAg                                                                     <  7:51009/149‑1 (MQ=255)

                                                                          ttCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTc                                                                   <  7:83633/149‑1 (MQ=255)

                                                                              tgcagATAAACGCACCGGGGGGGTCGATGATCCCCGCATCATTGGTGATACACCCCTCCCTCGCGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCaa                                                                 <  1:358462/144‑1 (MQ=255)

                                                                                   aTAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGAt                                                                             >  5:234841/1‑130 (MQ=255)

                                                                                   aTAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGAt                                                                             <  6:234841/130‑1 (MQ=255)

                                                                                      aaGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggccg                                                                                                                            >  3:321757/1‑80 (MQ=255)

                                                                                      aaGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTggccg                                                                                                                            <  4:321757/80‑1 (MQ=255)

                                                                                      aaGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAgctgct                                                       <  3:212510/149‑1 (MQ=255)

                                                                                       aGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGctc                                                      <  1:231808/149‑1 (MQ=255)

                                                                                             cgggcgggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCg                                                <  7:31085/149‑1 (MQ=255)

                                                                                                   ggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgc                                          >  4:425101/1‑149 (MQ=255)

                                                                                                   ggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgc                                          <  6:325981/149‑1 (MQ=255)

                                                                                                   ggTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgc                                          <  5:374378/149‑1 (MQ=255)

                                                                                                         tgatCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcg                                    <  2:6542/149‑1 (MQ=255)

                                                                                                         tgatCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcg                                    <  8:20697/149‑1 (MQ=255)

                                                                                                          gatCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcgc                                   >  6:313327/1‑149 (MQ=255)

                                                                                                            tCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGcctgcctg                                       <  2:401858/143‑1 (MQ=255)

                                                                                                            tCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGcctgcctg                                       >  1:401858/1‑143 (MQ=255)

                                                                                                             cACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                      >  2:44751/1‑95 (MQ=255)

                                                                                                             cACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                      <  1:44751/95‑1 (MQ=255)

                                                                                                              aCCGCATCATTGCTGATCCACCAC‑CCCACGGGGTATGCGTGCGCCTTCGCTGGGCGGGCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCa                                <  2:299253/148‑1 (MQ=255)

                                                                                                              aCCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAg                               >  4:247087/1‑149 (MQ=255)

                                                                                                                    tcatTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATGAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAg                         <  8:48266/149‑1 (MQ=255)

                                                                                                                    tcatTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAg                         >  7:257281/1‑149 (MQ=255)

                                                                                                                     catTGCTGATCCACCAC‑CCCAGGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGGCAACGAGAACCTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGCgg                        >  8:386002/1‑149 (MQ=255)

                                                                                                                      atTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAgg                        <  3:421283/148‑1 (MQ=255)

                                                                                                                      atTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGt                       <  8:257281/149‑1 (MQ=255)

                                                                                                                          cTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCa                   >  7:300901/1‑149 (MQ=255)

                                                                                                                                  ccac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATc           <  5:270398/149‑1 (MQ=255)

                                                                                                                                   cac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCg          >  7:352522/1‑149 (MQ=255)

                                                                                                                                   cac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCg          >  4:430497/1‑149 (MQ=255)

                                                                                                                                    ac‑cCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCgg         >  6:135943/1‑149 (MQ=255)

                                                                                                                                     c‑cccACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATATAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGt        <  4:107898/149‑1 (MQ=255)

                                                                                                                                       cccACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTa       <  8:132207/149‑1 (MQ=255)

                                                                                                                                        gcACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTa       >  1:99909/2‑148 (MQ=38)

                                                                                                                                            ggggTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAgcg  <  8:339325/149‑1 (MQ=255)

GTTTGGCTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCAC‑CCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCG  >  NC_002947/4599142‑4599429

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