BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,163,464	A→C	18.5%	S161R (AGC→CGC)	PP_3661 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,163,464	0	A	C	18.5%	420.9 / 19.7	229	S161R (AGC→CGC)	PP_3661	membrane protein
Reads supporting (aligned to +/- strand): ref base A (96/82); new base C (29/20); total (125/104)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

CGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAA  >  NC_002947/4163315‑4163613

cGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCa                                                                                                                                                       <  1:3522233/150‑1 (MQ=255)

    gCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcgc                                                                                                                                                   <  1:5112091/150‑1 (MQ=255)

      tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccc                                                                                                                                                 >  2:6581999/1‑150 (MQ=255)

      tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGccc                                                                                                                                                 >  2:4363288/1‑150 (MQ=255)

       ggTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccc                                                                                                                                                <  2:2923752/150‑1 (MQ=255)

        gtcgTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGATAGTGTTCTAGCTCGACCGCGATACCGCATTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTCGAATCGGGGGTGCGCAGCGccccc                                                                                                                                               <  2:6860875/147‑1 (MQ=255)

        gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccgat                                                                                                                                               <  1:6385647/150‑4 (MQ=255)

        gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccccc                                                                                                                                               <  2:3515192/150‑1 (MQ=255)

        gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccccc                                                                                                                                               <  2:7707144/150‑1 (MQ=255)

        gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGccccc                                                                                                                                               >  1:2440057/1‑150 (MQ=255)

         tAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccc                                                                                                                                              >  2:2008969/1‑150 (MQ=255)

            tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGcccccccccaa                                                                                                                                           >  2:1411719/1‑147 (MQ=255)

              aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGAATCGGCGCGGCGCCGCGcccccccgaaa                                                                                                                                         >  1:2708400/1‑145 (MQ=255)

               cctgTACCCTGCTGTTGGGCGTGGGGTGGGGCATGATTGTGTTCATGCTCGACCGCGATACCGCCTTGCTGACCGCCGGCGGCAGCGCCTTGTGGGGCGCGGCCGCCGTGCTGGCATTCGAATCGGGGCTGCGCTGCGCCCCCCATaaaa                                                                                                                                        <  2:4753158/147‑1 (MQ=255)

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

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

                cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCAAAAAg                                                                                                                                       >  2:3348956/1‑150 (MQ=255)

                   tACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCAACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAgcgc                                                                                                                                    <  2:5033659/150‑1 (MQ=255)

                     cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCa                                                                                                                                  >  2:2119443/1‑150 (MQ=255)

                       ctgctgTTGGGCCTGTGGTTTGGCTTGATAGTTTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCGGTGCTTGCATTCGAATCTGCGCTGCGAAGCGCCCCCCATAAAAGCGCCATg                                                                                                                                <  2:2703114/150‑1 (MQ=255)

                       ctgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAAGCGCCATg                                                                                                                                >  2:6789899/1‑150 (MQ=255)

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

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

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

                           tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCGAAAAGCGCCAAGGCAg                                                                                                                            >  2:2486332/1‑150 (MQ=255)

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

                                 gcgtgtggtgtgtATTTAATGTGTTAAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCGGGCTGCGCCTTGTGGGGCGCCGCGGCCGTGGTGGATTTCTTTTGGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                                      <  1:624263/147‑1 (MQ=255)

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

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

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

                                    tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCACTCGAATCGGCGGGGGGCAGCGCCCCCCCGAAAAGCGCCCAGGGAGGCGGAcgcg                                                                                                                   >  1:8262193/1‑150 (MQ=255)

                                     gtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCGCCCCCCAAAAAAGCGCCCAGGCCGACGGAACCgg                                                                                                                  >  1:5889909/1‑149 (MQ=255)

                                       ggtgtggCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAg                                                                                                                <  2:4584652/150‑1 (MQ=255)

                                         tgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGTAGTg                                                                                                              >  2:4654319/1‑150 (MQ=255)

                                             gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTgg                                                                                                          >  2:3183435/1‑149 (MQ=255)

                                              cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGtt                                                                                                         <  2:8116044/150‑1 (MQ=255)

                                              cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCCCCCCCCCAAAAAGCGCCAAGGCACTCGGAAGCGGAGTGCTGtt                                                                                                         >  1:7609527/1‑150 (MQ=255)

                                               aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTc                                                                                                        >  2:5893410/1‑150 (MQ=255)

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

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

                                                 gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCTAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCgg                                                                                                      >  2:5471920/1‑150 (MQ=255)

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

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

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

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

                                                          cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                             <  2:7292707/150‑1 (MQ=255)

                                                           aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCAAAAAGCGCCCAGGCAGGCGGAAGCGGAGTGGGGTTGGGTAGCGTGTCg                                                                                            >  1:4658097/1‑150 (MQ=255)

                                                           aaGCTCGACCGCGATACCGCCTGGCTGACCGCCGCCGGCAGCGCAATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCg                                                                                            <  2:8204800/150‑1 (MQ=255)

                                                             gCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                          >  2:5969860/1‑150 (MQ=255)

                                                                 gACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTc                                                                                      <  2:514790/150‑1 (MQ=255)

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

                                                                  aCCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTcc                                                                                     <  2:6182242/150‑1 (MQ=255)

                                                                   ccGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCGGTTGGGTAACTTGTCGATGTTACt                                                                                    >  1:4159193/1‑150 (MQ=255)

                                                                     gctaTACCGCCTGGCTGACCGCCGCGGGCAGCGCCATGTGCGGCGCCGCCGCCGTGCTGGCATTCGAATGGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCttt                                                                                  <  1:3551990/147‑1 (MQ=255)

                                                                      cgATACCGCCTTGCCGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTa                                                                                 <  2:6395767/150‑1 (MQ=255)

                                                                       gATCCCGCCTTGCTGACCGCCGCCGGCAGCGCCATCGGCGGCGCCGCCGCGGTGCTGGCATTCGACGCGGCGGTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAc                                                                                <  1:4363288/150‑1 (MQ=255)

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

                                                                             gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCCCCCCCCCTAAAAGCGCCCTGGCAGTCGGGAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCAcag                                                                          >  1:8574359/1‑148 (MQ=255)

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

                                                                                  gCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGAt                                                                     >  2:2090939/1‑150 (MQ=255)

                                                                                   cTGCCCGCCGCCGGCAGCGCCTTCGGGGGCGCCGCGGCGGTGCTGTATTTCGAATGGGCGGGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATc                                                                    <  1:3141193/150‑1 (MQ=255)

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

                                                                                     gACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCCTAAAAGCGCCATGGCCGTCGGAAGCGTAGTGCGGTTCGGTACCGTGTCGTTGTTCCTTTACCCACTGGCGATCaa                                                                  >  1:7152569/1‑150 (MQ=255)

                                                                                      aCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAt                                                                 >  2:2315016/1‑150 (MQ=255)

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

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

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

                                                                                          ccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTTTCGATGTTCCTTTAACCACAGGCGATCAAAGCCg                                                             >  2:4821869/1‑150 (MQ=255)

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

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

                                                                                                        aTCTGCGGCGCCGCCGCCGTGCTGTCATTCGACTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTg                                               <  1:6486076/150‑1 (MQ=255)

                                                                                                         tCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCATGTCGAAGTTCCTTTACCCACGGGCGATCAATGCCGGCTGGCTGCATCTgg                                              >  2:6453936/1‑150 (MQ=255)

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

                                                                                                               gcgcCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACAGGCGATCAATGCCGGCTGGCTGCATCTGGATACCa                                        >  2:5223855/1‑150 (MQ=255)

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

                                                                                                                cgccgccgcCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCGGTTCGGTACCTTGTCGATGTTCCTTTACCCACGGGCGAGCAAAGCCGGCTGGCTGCCACTGGCTACCaa                                       >  1:7905244/1‑149 (MQ=255)

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

                                                                                                                  ccgccgcggTGCTGGCTTTCGACGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                                     <  1:6843729/150‑1 (MQ=255)

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

                                                                                                                               ggaaTTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGGACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgct                        <  2:3859812/147‑1 (MQ=255)

                                                                                                                               ggCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgct                        >  2:2641051/1‑150 (MQ=255)

                                                                                                                                gatttCGACGCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctg                       <  1:6789899/147‑1 (MQ=255)

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

                                                                                                                                    tcgacccggcgcggcgCCGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg                   <  1:5515557/144‑1 (MQ=255)

                                                                                                                                    tCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATAACATGGGCGCCGGGCTGCTGCTgg                   >  2:6990457/1‑150 (MQ=255)

                                                                                                                                    tCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGGTCGGTACCTTTTCGATGTTCCTTTAACCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAAGGGCGCCGGGCTGCTGCTgg                   >  2:7539795/1‑150 (MQ=255)

                                                                                                                                      gattCGGCGCGGCCCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc                 <  1:1411719/147‑1 (MQ=255)

                                                                                                                                      gAATCGGCGCTGCGCCGCGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCAGTTCGGGACCTTGTCGATGTTCCATTACCCACAGGCGATCAAATCCCGCTGGCTGCATCAGGATACCAAGGGCGCCGGGCTGCTGCTGggc                 >  2:110622/1‑150 (MQ=255)

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

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

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

                                                                                                                                          gggCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGACGCGGAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCACTGGCGAGCAAAGCCGGCTGGCTGCATCTGGATACCAAGGGCGCCGGGCTGCTCCTGggcggcc             >  2:4390192/2‑149 (MQ=255)

                                                                                                                                          gggCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa             <  2:7841002/149‑1 (MQ=255)

                                                                                                                                            gcgcTGCGCCGCGCCCCCCCTAACAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt           <  1:3390357/150‑1 (MQ=255)

                                                                                                                                             cggtgcgctgcgcCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTg          <  1:8037114/147‑1 (MQ=255)

                                                                                                                                                tGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTcc       <  1:6581999/150‑1 (MQ=255)

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

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

                                                                                                                                                  cgcAGCGCCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCAc     <  2:1084943/150‑1 (MQ=255)

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

CGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAA  >  NC_002947/4163315‑4163613

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