BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,163,474	A→C	15.5%	H164P (CAT→CCT)	PP_3661 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,163,474	0	A	C	15.5%	525.9 / 36.3	231	H164P (CAT→CCT)	PP_3661	membrane protein
Reads supporting (aligned to +/- strand): ref base A (92/98); new base C (37/4); total (129/102)
Fisher's exact test for biased strand distribution p-value = 4.51e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

GTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGGCGCAGG  >  NC_002947/4163326‑4163623

gTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAt                                                                                                                                                      <  2:1291910/150‑1 (MQ=255)

 tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGGATGTAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGTCCGCCGCCGGCAGCGCGATCTGCGGCGCCGCCGCCGTGCTGGCATTCGTTTCGGGGCTGGGGAGGGCCCCCCATa                                                                                                                                                     >  2:5834858/1‑150 (MQ=255)

     cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAAg                                                                                                                                                 >  2:1881735/1‑150 (MQ=255)

     cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAAg                                                                                                                                                 >  2:1502705/1‑150 (MQ=255)

     cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCAAAAAg                                                                                                                                                 >  2:7779945/1‑150 (MQ=255)

           ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGAATCGGCGCTGCGCCGCGCCCCCCCACAAAGCGccct                                                                                                                                           >  2:3286371/1‑148 (MQ=255)

            ctgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATg                                                                                                                                          <  2:2702607/150‑1 (MQ=255)

            ctgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCTAAAAGCGCCaag                                                                                                                                          >  2:7371009/1‑148 (MQ=255)

             tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATCAAAGCGCCATgg                                                                                                                                         >  2:6236487/1‑150 (MQ=255)

             tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCAAgg                                                                                                                                         >  2:5536963/1‑150 (MQ=255)

               ctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCa                                                                                                                                       >  2:888152/1‑150 (MQ=255)

                  aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCCGCAAAGCGCCCTGGCAGTc                                                                                                                                    >  1:6027627/1‑150 (MQ=255)

                  aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTc                                                                                                                                    >  2:2853672/1‑150 (MQ=255)

                       ccTGTGGTGTGGCATGAAGGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCCGCTGCGGCGCCGCCGCCGTGCGGGCTGGCGCGTCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCGTGGCAGGCGGTAg                                                                                                                               >  2:3362347/1‑150 (MQ=255)

                         tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGTGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCCCCCCCCATAAAAGCGCCAAGGCACTCGGAAGcc                                                                                                                             >  2:833284/1‑149 (MQ=255)

                         tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCGACAAGCGCCAAGGCAGTCGGAAGCg                                                                                                                             >  2:5185092/1‑150 (MQ=255)

                         tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCg                                                                                                                             >  2:4272980/1‑150 (MQ=255)

                           tggtgtggCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGGTGACCGCCGCGGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                                           <  2:2244105/150‑1 (MQ=255)

                                tggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCAAAAAGAGCCAAGGCAGTCGGAAGCGGAGTGCt                                                                                                                      >  1:7015737/1‑150 (MQ=255)

                                  gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGGAGTGgtgt                                                                                                                    >  2:370508/1‑150 (MQ=255)

                                    aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCACCGTGCTGGCATTCGAATCGGCGCTGCGCCCCGCCCCCCCCAAAAACCCCCAGGCAGGCGGAAGCGTCTGGCTGtta                                                                                                                  >  1:323636/1‑149 (MQ=255)

                                     tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCg                                                                                                                 <  2:7211509/150‑1 (MQ=255)

                                     tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGCGCAGCGCCCCCCACAAAAGCGCCAAGGCAGTCGGAAGCGGAGTGCTGTTCg                                                                                                                 >  2:4645590/1‑150 (MQ=255)

                                        aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAAGCGCCAGGGCCGTCGGAAGCGGAGTGGGGGTCGGTa                                                                                                              >  2:3958965/1‑150 (MQ=255)

                                        aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTa                                                                                                              >  2:7934418/1‑150 (MQ=255)

                                        aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCAAAAAGCGCCCCGGCAGTCGGAAGCCTAGTGATGTTCGGTa                                                                                                              >  1:4286011/1‑150 (MQ=255)

                                           tgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAACCGCCCTTGCAGTCGGAATCGTAGTGCAGTTAGGTAACt                                                                                                           >  2:735459/1‑150 (MQ=255)

                                           tgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCAAAAAAGAGCCCAGGCAGTCGGAAGCGTAGTGCTGTTCGGTAcca                                                                                                           >  1:3544418/1‑149 (MQ=255)

                                            gtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAAAAAAACGCCCAGGCCGTCGGAAGCGGAGTGCTGTTCGGGAccat                                                                                                          >  1:6391217/1‑148 (MQ=255)

                                             ttCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCATGGCCGTCGGAAGCGTAGTGCTGGTCGGTAccggg                                                                                                         >  2:1430638/1‑147 (MQ=255)

                                               cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCAAAAAAACCCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCATGTc                                                                                                       >  1:6970549/1‑150 (MQ=255)

                                                  gCTCGCCCGCTATACGGCTTTGCTGCCCGCCGCCGGCCGCGCCTTCTGGGGCGCCGCCGCGGTGCTGGATTTCGACGGGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                                    <  1:436838/150‑1 (MQ=255)

                                                  gCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCCTTCGAATCGGCGCGGCGCCGCGCCCCCCCGAAAAGCGCCCGGGCAGTCGGACGCGGCGTGCGGGTCGGTACCGTGTCGag                                                                                                    >  2:4886758/1‑149 (MQ=255)

                                                    tCGCCCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCAGCTGCGGCGCCGCCGCCGTGCTGGCATTCGACTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                                  <  2:2027601/150‑1 (MQ=255)

                                                       aCCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCCCAAAGCGCCCGGGCCCGCGGACGCGGCGTGCGGGTCGGTGCCAGGTCgcggttgc                                                                                               >  1:4759095/1‑143 (MQ=255)

                                                        cggcgATCCCGCCTTGCTGACCGCCGCGGGCACCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCt                                                                                              <  1:5076869/148‑1 (MQ=255)

                                                               aCCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGGGCCGCGCCCCCCCGAAAAGCCCCCAGGCAGGCGGGAGCGGAGTGCTGTTGGGTGCCAAGTGGATGTTGCGTTaaaca                                                                                       >  1:8202740/1‑146 (MQ=255)

                                                                ccGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCAc                                                                                      >  2:3059261/1‑150 (MQ=255)

                                                                 cGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTACTTTACCCACt                                                                                     >  2:6659418/1‑150 (MQ=255)

                                                                   ccTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCGAAAAGCGCCAGGGCAGTCGGAAGCGGAGTGCTGTTCGGGACCTTGTCGATGTTACTTTAACCACAgg                                                                                   >  2:360698/1‑150 (MQ=255)

                                                                        cTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCTCCGCGTTGGTGGATTTCTAATGGGCGATGGGCAGCGGCCCCCATATAAGCTCCATGTAAGTCGGAAGCGTAGTGCTGTTGTGTTCCTTGTCGATGTTCCTTTACCCTCTggggttg                                                                              <  2:8639167/150‑6 (MQ=255)

                                                                         tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCAAAAAGCGCCAGGGCAGTCGGCAGCGGAGTGCTGGTGGGTACCGTGTCGGTGTTGCGTTTCCCACTGGCGCTcc                                                                             >  1:7775980/1‑149 (MQ=255)

                                                                           aCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAt                                                                           <  2:6765148/150‑1 (MQ=255)

                                                                            ccgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGGTCGGTACCGTGTCGATGTTCCTTTACCCACTGGCGATCAATg                                                                          >  2:5377132/1‑150 (MQ=255)

                                                                            ccgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCCCCCCCCCAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTTCCGTGTCGTTGGTGCGTTAACCACAGGCGATCaaat                                                                          >  1:5569899/1‑148 (MQ=255)

                                                                            ccgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAGAAAACCGCCATGGCCGTCGGAAGCGTAGTGCGGTTCGGGACCGTGGCGATGTTCCGTTAACCACAGGCGATCAATg                                                                          >  2:2395615/1‑150 (MQ=255)

                                                                            cccccgccgGCCGCGCCAGCTGCGGCGCCGCCGCCGTGCGGGCAGTCGACTCGGCGCGGCGCCCCGCCCCCCCAAAAACCGCCATGGCAGGCGGAACCCGCCGGCGGGTCGGGCCCGGGTCGAGGGTCCGTTACCCACCGGCGAGCCATg                                                                          >  2:6009184/4‑150 (MQ=255)

                                                                             cgccgcggGCTGCGCCATGGGCGGCGCCGCCGCGGTGGTGGCTTTCGAAGCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGc                                                                         <  1:8516815/150‑1 (MQ=255)

                                                                             cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCCGCGCAGCGCCCCCCAAAAAAGCCCCAAGGCACTCCGCAGCGTAGTGCGGTTCGGTACCATGTCGATGTTCCTTTAACCACAGGCGATCaaatc                                                                         >  2:4806146/1‑147 (MQ=255)

                                                                                          gcCATCTGCGGCGCCGCCGCCGTGCTGGATTTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCAt                                                            <  1:4183732/150‑1 (MQ=255)

                                                                                            catggggggcgccggcgccGTGGTGGAATTCTAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                                                          <  1:6488893/143‑1 (MQ=255)

                                                                                            cATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCAAAAAAGCGCCAGGGCAGTCGGAAGCGGAGTGCGGTTCGGTACCGTGTCGGTGTTCCGTTACCCCCTGGCGATCAAAGCCGGCTGGCGGCctct                                                          >  1:2087003/1‑150 (MQ=255)

                                                                                            cATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                                                          <  2:1682133/150‑1 (MQ=255)

                                                                                               cggcggcgCCGCCGCCGGGGTGGCTTTCGAATGGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                                       <  1:5742896/148‑1 (MQ=255)

                                                                                               cTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                                       <  2:5958158/150‑1 (MQ=255)

                                                                                                tGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATCAAAGCGCCCTGGCAGTCGGACGCGTAGTGCGGGTCGGTACCTTGTCGATGGTCCTGTACCCACTGGCGCTCAATGCCGGCTGGCGGCATCTggct                                                      >  2:8132061/1‑148 (MQ=255)

                                                                                                   ggcgCAGCCGCCGTCCTGGCATTCGAATCGGCGAAGCCCCGCGCCCCCCATAAAATCGCCATGGCAGTCGGAAGCGTAGTGCTGTTAGGTCACTTGTCGATGTTCCTTTAACCACTGGCGATCAATGCCGGCAGGCTGCATCTGGATAcc                                                   >  2:6450711/1‑150 (MQ=255)

                                                                                                    gcgcCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAATCGCCATGGGAGTCGGAAGCGTAGTGCTGTTCGGTACCTTTTCGATGTTACTTTACCCACAGGCGATAAATGACGGCTGGCTGCATAAGGATACCa                                                  >  2:878322/1‑150 (MQ=255)

                                                                                                      gccgccgcgGTGGTGGATTTCGACGCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATg                                                <  1:6675098/150‑1 (MQ=255)

                                                                                                       ccgccgcggTGGTGGATTTCTCAGGGGCGCGGCCCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                                               <  1:3000097/150‑1 (MQ=255)

                                                                                                       ccgccgcggTGCTGTCATTCGACGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                                               <  1:6843652/150‑1 (MQ=255)

                                                                                                       ccgccgccgTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                                               <  2:4385851/150‑1 (MQ=255)

                                                                                                                   tGGCATTCTAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgc                                   <  2:3877509/150‑1 (MQ=255)

                                                                                                                   tGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACAGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgc                                   >  2:5062140/1‑150 (MQ=255)

                                                                                                                       aTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCATGTCGATGTTCCTTTACCCACTGGCGATCAAAGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgct                               >  2:263502/1‑150 (MQ=255)

                                                                                                                       aTTCGAATCGGCGCGGCGCAGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCACTGGCGATCAATTCCGGCTGGCTGCATCTGGATAACAAGGGCGCCGGgctgctgcg                               >  2:1104182/1‑149 (MQ=255)

                                                                                                                        ggcgacgcggcgcggcgCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg                              <  1:118272/143‑1 (MQ=255)

                                                                                                                           taaTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc                           <  1:5834858/149‑1 (MQ=255)

                                                                                                                           gactCGGCGCTGCCCAGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc                           <  1:370508/147‑1 (MQ=255)

                                                                                                                           gAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc                           <  2:844929/150‑1 (MQ=255)

                                                                                                                           gAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc                           <  1:2854134/150‑1 (MQ=255)

                                                                                                                            ccgcggcgcggcgCCGCGCCCCCCAGACCCGCGCCAGGGCAGGCGGAAGCGTAGTGCTGTTCGGGACCGTGTCGATGTTCCTTTACCCACTGGCGAGCACTGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg                          <  2:7598003/147‑1 (MQ=255)

                                                                                                                            attCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg                          <  2:4678891/148‑1 (MQ=255)

                                                                                                                             gggggCGCTGCGCAGCGCCCCCCAAAAAATCGCCAAGGCCGTCGGAAGCGTAGTTCTGTTCGGTACCTTGTCGATGTTCCTTTAACCACAGGCGAACAATTCCGGCATGCTGCATCAGGATAACAAGGGCGCCGGGCTGCAGCTGggcgg                         >  2:2103402/4‑150 (MQ=255)

                                                                                                                             aTCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcgg                         <  1:1354796/150‑1 (MQ=255)

                                                                                                                             aTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTACTTTACCCACTTGCGAACAAAGCCGGCTGGCTGCATCTGGATACCAAGGGCGCCGGGCTGCTGCAGggcgg                         >  1:6169553/1‑150 (MQ=255)

                                                                                                                               cGGCGCTGCGCAGTGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa                       <  2:6108182/150‑1 (MQ=255)

                                                                                                                                 gcgcgccgccgcCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt                     <  1:3390186/144‑1 (MQ=255)

                                                                                                                                   ggtgcgctgcgcCCCCAATTAAAGCGCCTTGGCAGTCGGTCGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACGGGCGATCTATGCCGGCTGGCTGCATCTGGATTCCATGGGCGCCGGGCTGCTGCTGGGCGTCACTGt                   <  2:729710/148‑1 (MQ=255)

                                                                                                                                    cTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTc                  >  2:5237019/1‑150 (MQ=255)

                                                                                                                                      gcgcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCa                <  1:2008969/150‑1 (MQ=255)

                                                                                                                                       cgccgcGCCCCCCATACAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCAc               <  1:3348956/150‑1 (MQ=255)

                                                                                                                                       cgcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCAc               <  1:7257827/150‑1 (MQ=255)

                                                                                                                                       cgcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCAc               >  2:6962704/1‑150 (MQ=255)

                                                                                                                                        gcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACg              >  1:2661276/1‑150 (MQ=255)

                                                                                                                                        gcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACg              >  2:6775634/1‑150 (MQ=255)

                                                                                                                                         caccgcCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGa             <  1:6261350/147‑1 (MQ=255)

                                                                                                                                           gcgcCCCCCATACACGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCCCTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAg           <  2:2521392/150‑1 (MQ=255)

                                                                                                                                            cgcCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGt          >  2:4638185/1‑150 (MQ=255)

                                                                                                                                             gcCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTg         <  1:3860783/150‑1 (MQ=255)

                                                                                                                                              ccccccATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCCCGAAGTgg        >  2:7905446/1‑150 (MQ=255)

                                                                                                                                              ccccccATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTgg        <  1:793381/150‑1 (MQ=255)

                                                                                                                                                ccccATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGgcg      <  2:7077765/150‑1 (MQ=255)

                                                                                                                                                ccccATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGgcg      >  1:8760230/1‑150 (MQ=255)

                                                                                                                                                ccccATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGgcg      >  1:1523295/1‑150 (MQ=255)

                                                                                                                                                 cccATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATTGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGgcgc     >  2:7478116/1‑150 (MQ=255)

                                                                                                                                                 cccATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGgcgc     >  2:650545/1‑150 (MQ=255)

                                                                                                                                                 cccATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGgcgc     >  1:2485336/1‑150 (MQ=255)

                                                                                                                                                  ccAAAAAAGCGCCATGGAAGTCGTAAGCGTAGTGCTGTTCGGTAACTAGTCGATGTTCCTTTACCCAATGGCGATCAATGCCTGCTGGCTGCATCTCGAAACCATGGGCGCCCGGCTGCTGCTGGCCGGCAAAGTCCACGAAGAGGAGCa    >  2:4434165/1‑150 (MQ=255)

                                                                                                                                                    aTAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGGCGCAgg  >  1:8001293/1‑150 (MQ=255)

GTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGGCGCAGG  >  NC_002947/4163326‑4163623

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