BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	433,850	C→G	7.1%	V451V (GTG→GTC)	glcB ←	malate synthase G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	433,850	0	C	G	7.1%	691.2 / 32.2	210	V451V (GTG→GTC)	glcB	malate synthase G
Reads supporting (aligned to +/- strand): ref base C (97/98); new base G (10/5); total (107/103)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 221 total aligned reads displayed.

GGCCAGGCCAACGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTC  >  NC_002947/433703‑433999

 gccaggccaACGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCAcg                                                                                                                                                    >  1:659807/1‑150 (MQ=255)

    aggccaACGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACgcgc                                                                                                                                                 >  1:1146725/1‑150 (MQ=255)

      gccaACGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTc                                                                                                                                               >  1:2770906/1‑150 (MQ=255)

       ccaACGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTcg                                                                                                                                              <  1:1587019/150‑1 (MQ=255)

          aCGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTcggcg                                                                                                                                           <  2:878499/150‑1 (MQ=255)

            gTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCggcggc                                                                                                                                         >  2:3897818/1‑150 (MQ=255)

            gTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCggcggc                                                                                                                                         <  1:4931595/150‑1 (MQ=255)

            gTCGACGTTGTTGTGTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCggcggc                                                                                                                                         >  2:3380417/1‑150 (MQ=255)

                 cGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCt                                                                                                                                    <  1:3981639/150‑1 (MQ=255)

                    tgttgttTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGa                                                                                                                                 >  1:3191725/1‑150 (MQ=255)

                    tgttgttTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGa                                                                                                                                 >  1:3610532/1‑150 (MQ=255)

                    tgttgttTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGa                                                                                                                                 <  2:321171/150‑1 (MQ=255)

                        gttTTCGTAGGCGCCGTTCGACTTCTCGTTCTTCATGGCTCCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTTGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTTTTGATGAACACCACGCGCTCGGCTGCTGCCTTGATGCa                                                                                                                             <  1:4384162/150‑1 (MQ=255)

                        gttTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCa                                                                                                                             <  1:5434282/150‑1 (MQ=255)

                         ttttCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAg                                                                                                                            <  1:1283181/150‑1 (MQ=255)

                               tAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTg                                                                                                                      <  1:2654988/150‑1 (MQ=255)

                               tAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTg                                                                                                                      <  1:299451/150‑1 (MQ=255)

                                aGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGa                                                                                                                     <  2:1120886/150‑1 (MQ=255)

                                  gcgcCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCtt                                                                                                                                   >  1:1079679/1‑134 (MQ=255)

                                   cgcCGATCCACGTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGt                                                                                                                  <  2:4265317/150‑1 (MQ=255)

                                    gcCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGa                                                                                                                     >  1:1001194/1‑146 (MQ=255)

                                            aCTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTgg                                                                                                         >  1:4825029/1‑150 (MQ=255)

                                            aCTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTgg                                                                                                         >  1:4549856/1‑150 (MQ=255)

                                            aCTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTgg                                                                                                         >  2:5251070/1‑150 (MQ=255)

                                            aCTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATAGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTgg                                                                                                         <  1:2265192/150‑1 (MQ=255)

                                               tctcGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTAc                                                                                                      <  1:1238501/150‑1 (MQ=255)

                                                ctcGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACg                                                                                                     >  1:1409251/1‑150 (MQ=255)

                                                ctcGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACg                                                                                                     <  1:1200111/150‑1 (MQ=255)

                                                ctcGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACg                                                                                                     >  2:401808/1‑150 (MQ=255)

                                                 tcGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACgg                                                                                                    >  2:3634091/1‑150 (MQ=255)

                                                    ttcttcATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGt                                                                                                 >  2:3062710/1‑150 (MQ=255)

                                                    ttattcATGGCACCTGTGCGCACCACGGCGCCGGCTTCCGGCGAGGTGTGGCTTTCGTCGGCAGTGAGGTCAAGGAGGCCGGTGGGGGTGAACACCACGCGCTCGGCGGCGGCCTTGGTGGAGGACGTGAGGTGGGCCGGGGTGGGGCGt                                                                                                 >  2:4453984/4‑150 (MQ=255)

                                                         cATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGa                                                                                                                                 <  1:970854/113‑1 (MQ=255)

                                                          aTGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCg                                                                                           >  2:5433892/1‑150 (MQ=255)

                                                              cACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCa                                                                                       <  2:1520211/150‑1 (MQ=255)

                                                               aCCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCat                                                                                      >  2:1931543/1‑150 (MQ=255)

                                                                  tttGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCatgat                                                                                   >  2:4944688/1‑150 (MQ=255)

                                                                     gcgcACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGcc                                                                                <  1:4944688/150‑1 (MQ=255)

                                                                     gcgcACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGcc                                                                                <  2:1653539/150‑1 (MQ=255)

                                                                       gcACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGa                                                                              >  2:1118171/1‑150 (MQ=255)

                                                                         accacGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGAcc                                                                            <  2:1388522/150‑1 (MQ=255)

                                                                         accacGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGAcc                                                                            >  1:4087711/1‑150 (MQ=255)

                                                                          ccacGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCt                                                                           <  2:1822831/150‑1 (MQ=255)

                                                                                cgcCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACgg                                                                     <  1:2028319/150‑1 (MQ=255)

                                                                                 gcCAGCTTCCATCGAGGTATGGATTTCGTCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTc                                                                                               >  2:4385899/1‑123 (MQ=255)

                                                                                  ccGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGgtg                                                                   <  1:5251049/150‑1 (MQ=255)

                                                                                   cGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGgtgt                                                                  >  2:4939896/1‑150 (MQ=255)

                                                                                    ggCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGtt                                                                 <  2:944505/150‑1 (MQ=255)

                                                                                     gCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTg                                                                >  1:3860075/1‑150 (MQ=255)

                                                                                          cATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCgg                                                           <  2:611512/150‑1 (MQ=255)

                                                                                          cATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCgg                                                           >  2:5415529/1‑150 (MQ=255)

                                                                                          cATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCgg                                                           >  1:3625430/1‑150 (MQ=255)

                                                                                           aTCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGc                                                          >  2:2776986/1‑150 (MQ=255)

                                                                                             cGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCAt                                                        <  2:2099308/150‑1 (MQ=255)

                                                                                              gAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATg                                                       >  2:4205425/1‑150 (MQ=255)

                                                                                               aGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGc                                                      <  2:1379755/150‑1 (MQ=255)

                                                                                                 gtgtGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGcc                                                     <  1:533572/149‑1 (MQ=255)

                                                                                                 gtgtGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGcc                                                     >  2:533678/1‑149 (MQ=255)

                                                                                                  tgtgGATTTCATCGCCAGTGCGGTCAAGGTAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGTGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCca                                                   >  2:385502/1‑150 (MQ=255)

                                                                                                  tgtgGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCa                                                                                                                             <  1:2113791/76‑1 (MQ=255)

                                                                                                   gtgGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCcag                                                  >  2:2492352/1‑150 (MQ=255)

                                                                                                     ggATTTCGTCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTCCTCATCCATGATGCCGATTTTCAGTGTGTTGGGCTTCATACCaagca                                                <  1:5369012/150‑1 (MQ=255)

                                                                                                     ggATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCcagca                                                <  2:515348/150‑1 (MQ=255)

                                                                                                      gATTTCGTCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTCCTCATCCATGATGCCGATTTTCAGTGTGTTGGGCTTCATACCaagcag                                               >  1:160502/1‑150 (MQ=255)

                                                                                                      gATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCcagcag                                               >  1:3923005/1‑150 (MQ=255)

                                                                                                       aTTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAgg                                              >  2:476326/1‑150 (MQ=255)

                                                                                                         ttCGTCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTc                                                                                               >  1:3288647/1‑99 (MQ=255)

                                                                                                           cgtCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTCCTCATCCATGATGCCGATTTTCAGTGTGTTGGGCTTCATACCAAGCAGCTCtt                                          <  1:1825641/148‑1 (MQ=255)

                                                                                                             tCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCg                                        >  1:2956846/1‑150 (MQ=255)

                                                                                                             tCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTCCTCATCCATGATGCCGATTTTCAGTGTGTTGGGCTTCATACCAAGCAGCTCTTCg                                        >  1:545625/1‑150 (MQ=255)

                                                                                                               gCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTCCTCATCCATGATGCCGATTTTCAGTGTGTTGGGCTTCATACCAAGCAGCTCTTCggt                                      >  2:2345917/1‑148 (MQ=255)

                                                                                                                ccAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTCCTCATCCATGATGCCGATTTTCAGTGTGTTGGGCTTCATACCAAGCAGCTCTTCggtg                                     >  1:5025322/1‑147 (MQ=255)

                                                                                                                   gTGCGGTCAAGGAAGCCGGGGTTGCTGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTgg                                  <  1:2581345/150‑1 (MQ=255)

                                                                                                                     gCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCt                                <  1:169951/150‑1 (MQ=255)

                                                                                                                     gCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCt                                <  1:3478712/150‑1 (MQ=255)

                                                                                                                         tCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAg                            <  2:1705585/150‑1 (MQ=255)

                                                                                                                         tCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAg                            >  1:4388050/1‑150 (MQ=255)

                                                                                                                         tCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAg                            >  1:4128625/1‑150 (MQ=255)

                                                                                                                          cAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGCCTTCGACCTGGCTGAAGa                           <  1:1120062/150‑1 (MQ=255)

                                                                                                                           aaGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCatgat                                                                                   >  1:3813256/1‑93 (MQ=255)

                                                                                                                           aaGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGAt                          <  1:89473/150‑1 (MQ=255)

                                                                                                                             ggAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGtt                                                                 <  2:1001494/109‑1 (MQ=255)

                                                                                                                             ggAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGAtct                        >  1:2716422/1‑150 (MQ=255)

                                                                                                                             ggAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGAtct                        >  2:796865/1‑150 (MQ=255)

                                                                                                                              gAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGAtctc                       >  1:2410380/1‑150 (MQ=255)

                                                                                                                                 gCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGt                    >  2:3706390/1‑150 (MQ=255)

                                                                                                                                 gCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGt                    <  2:2754275/150‑1 (MQ=255)

                                                                                                                                   cGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACg                                                                      >  2:1768523/1‑98 (MQ=255)

                                                                                                                                   cGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACg                                                                      <  1:1768523/98‑1 (MQ=255)

                                                                                                                                    ggTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGc                 <  1:4004692/150‑1 (MQ=255)

                                                                                                                                     gtgtTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCg                >  2:5007949/1‑150 (MQ=255)

                                                                                                                                          gatgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGc           <  1:1058041/150‑1 (MQ=255)

                                                                                                                                          gatgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGc           >  2:2810356/1‑150 (MQ=255)

                                                                                                                                           atgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGcc          >  2:4730033/1‑150 (MQ=255)

                                                                                                                                           atgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGcc          <  1:4606711/150‑1 (MQ=255)

                                                                                                                                              aaCACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACt       >  1:2537826/1‑150 (MQ=255)

                                                                                                                                              aaCACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACt       <  1:1466788/150‑1 (MQ=255)

                                                                                                                                               acacCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACtt      >  2:969283/1‑150 (MQ=255)

                                                                                                                                                caccacGCTCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTc     >  1:4300007/1‑150 (MQ=255)

                                                                                                                                                caccacGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTc     <  2:3191725/150‑1 (MQ=255)

GGCCAGGCCAACGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTCCTC  >  NC_002947/433703‑433999

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