BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,189,953	G→A	7.8%	R163H (CGC→CAC) ‡	PP_4567 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,189,953	0	G	A	7.8%	896.7 / 44.6	297	R163H (CGC→CAC) ‡	PP_4567	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (129/144); new base A (11/12); total (140/157)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 317 total aligned reads displayed.

CTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCC  >  NC_002947/5189805‑5190102

cTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTgcgc                                                                                                                                                      >  2:812614/1‑150 (MQ=255)

cTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTgcgc                                                                                                                                                      >  2:6263219/1‑150 (MQ=255)

       tCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGcgccacgcca                                                                                                                                               >  2:1673498/1‑150 (MQ=255)

           cTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATccc                                                                                                                                           >  2:5745008/1‑150 (MQ=255)

            tGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATcccc                                                                                                                                          <  2:7956912/150‑1 (MQ=255)

               gcCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGaa                                                                                                                                       <  2:305057/150‑1 (MQ=38)

               gcCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGaa                                                                                                                                       <  1:966414/150‑1 (MQ=38)

               gcCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGaa                                                                                                                                       <  1:7714015/150‑1 (MQ=38)

               gcCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGaa                                                                                                                                       <  1:5525138/150‑1 (MQ=255)

                ccTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAag                                                                                                                                      <  2:3568185/150‑1 (MQ=255)

                ccTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAag                                                                                                                                      >  1:1013119/1‑150 (MQ=255)

                 cTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAaga                                                                                                                                     <  1:7122363/150‑1 (MQ=255)

                       aaCCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCg                                                                                                                               >  2:1697421/1‑150 (MQ=255)

                             aCAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAggg                                                                                                                         >  1:4656709/1‑150 (MQ=255)

                               aaGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGcc                                                                                                                       >  1:7553704/1‑150 (MQ=255)

                                 gCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTa                                                                                                                     <  2:5585010/150‑1 (MQ=255)

                                 gCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTa                                                                                                                     >  1:7601349/1‑150 (MQ=255)

                                     tGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCAGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCt                                                                                                                 >  1:2096048/1‑150 (MQ=255)

                                       aaGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCGCCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTc                                                                                                               >  2:6089100/1‑150 (MQ=255)

                                       aaGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTc                                                                                                               >  1:6237801/1‑150 (MQ=255)

                                          aTGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGc                                                                                                            <  1:8489883/150‑1 (MQ=255)

                                             ttGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGATGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAg                                                                                                         <  1:8800532/150‑1 (MQ=255)

                                             ttGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAg                                                                                                         <  2:7899857/150‑1 (MQ=255)

                                               ggCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGgc                                                                                                       >  1:6667758/1‑150 (MQ=255)

                                               ggCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGgc                                                                                                       >  1:4867808/1‑150 (MQ=255)

                                                    tCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGcgcgcg                                                                                                  >  1:3542306/1‑150 (MQ=255)

                                                     cGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCgcgg                                                                                                 <  1:7963169/150‑1 (MQ=255)

                                                      gAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCgcggc                                                                                                <  2:7702767/150‑1 (MQ=255)

                                                      gAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCgcggc                                                                                                <  1:1651539/150‑1 (MQ=255)

                                                      gAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCgcggc                                                                                                >  2:5967983/1‑150 (MQ=255)

                                                       aGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCgcggcg                                                                                               <  1:6547468/150‑1 (MQ=255)

                                                       aGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCgcggcg                                                                                               <  1:6721986/150‑1 (MQ=255)

                                                        ggTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGcggcgg                                                                                              <  1:7051629/150‑1 (MQ=255)

                                                         gTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCggg                                                                                             >  2:1484403/1‑150 (MQ=255)

                                                         gTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCggg                                                                                             >  1:272325/1‑150 (MQ=255)

                                                         gTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCggg                                                                                             >  2:1376583/1‑150 (MQ=255)

                                                          tCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGc                                                                                            >  1:1127195/1‑150 (MQ=255)

                                                            ccgccgCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCtt                                                                                          >  1:297275/1‑150 (MQ=255)

                                                                  cAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGaaa                                                                                    <  1:3190847/150‑1 (MQ=255)

                                                                    aCGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAaca                                                                                  <  1:4468517/150‑1 (MQ=255)

                                                                    aCGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAaca                                                                                  <  2:3886815/150‑1 (MQ=255)

                                                                        tcgctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGaa                                                                              >  1:1727482/1‑150 (MQ=255)

                                                                          gctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGCAACACGAAca                                                                            >  2:8759470/1‑150 (MQ=255)

                                                                           ctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAcag                                                                           >  2:5466191/1‑150 (MQ=255)

                                                                            tcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAcagc                                                                          >  2:1808402/1‑150 (MQ=255)

                                                                            tcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAcagc                                                                          >  2:2583970/1‑150 (MQ=255)

                                                                            tcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAcagc                                                                          >  1:2376504/1‑150 (MQ=255)

                                                                              gATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAcagcag                                                                        <  2:8478761/150‑1 (MQ=255)

                                                                               aTACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACagcagc                                                                       <  2:1076002/150‑1 (MQ=255)

                                                                                  cccTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCt                                                                    <  2:8660905/150‑1 (MQ=255)

                                                                                   ccTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTa                                                                   >  2:355411/1‑150 (MQ=255)

                                                                                   ccTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTACGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTa                                                                   >  1:693030/1‑150 (MQ=255)

                                                                                     tGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTAcc                                                                 >  2:2208291/1‑150 (MQ=255)

                                                                                       ctctACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCAc                                                               >  1:7857798/1‑150 (MQ=255)

                                                                                          tactacCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGAt                                                            >  2:3126410/1‑150 (MQ=255)

                                                                                          tactacCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGAt                                                            >  1:3297332/1‑150 (MQ=255)

                                                                                             tacCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAg                                                         >  1:3783565/1‑150 (MQ=255)

                                                                                              acCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGc                                                        <  1:1327939/150‑1 (MQ=255)

                                                                                              acCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGc                                                        <  1:2209244/150‑1 (MQ=255)

                                                                                               ccTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACagcagc                                                                       <  2:4813226/134‑1 (MQ=255)

                                                                                               ccTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGcc                                                       >  2:6246325/1‑150 (MQ=255)

                                                                                               ccTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGcc                                                       <  2:8296854/150‑1 (MQ=255)

                                                                                                 tGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGcccc                                                     >  2:3159700/1‑150 (MQ=255)

                                                                                                  gCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCa                                                    <  2:4656709/150‑1 (MQ=255)

                                                                                                  gCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCa                                                    <  2:4919256/150‑1 (MQ=255)

                                                                                                    tCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATg                                                  >  1:3095895/1‑150 (MQ=255)

                                                                                                     cACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGc                                                 >  1:1166114/1‑150 (MQ=255)

                                                                                                        ccAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCca                                              >  2:2272648/1‑150 (MQ=255)

                                                                                                         cAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCcag                                             >  2:4291303/1‑150 (MQ=255)

                                                                                                         cAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCcag                                             >  2:8893899/1‑150 (MQ=255)

                                                                                                          aGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCcagc                                            >  2:7129535/1‑150 (MQ=255)

                                                                                                          aGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCcagc                                            >  2:3786611/1‑150 (MQ=255)

                                                                                                            cgcgAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCcagcag                                          >  2:657787/1‑150 (MQ=255)

                                                                                                               gAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCt                                       >  1:8607670/1‑150 (MQ=255)

                                                                                                                      gCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGt                                <  2:2876936/150‑1 (MQ=255)

                                                                                                                      gCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGt                                <  2:2169307/150‑1 (MQ=255)

                                                                                                                      gCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGGTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGt                                <  1:4981058/150‑1 (MQ=255)

                                                                                                                        ctgctgGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGc                              <  2:4121779/150‑1 (MQ=255)

                                                                                                                        ctgctgGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGc                              <  1:4738694/150‑1 (MQ=255)

                                                                                                                          ggaactcTTCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCa                            >  1:5906305/7‑150 (MQ=37)

                                                                                                                          gctgGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCa                            <  1:2146970/150‑1 (MQ=255)

                                                                                                                                 ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACATGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGtt                     >  1:4918879/1‑150 (MQ=255)

                                                                                                                                 ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCCGGCGGtt                     >  2:8484601/1‑150 (MQ=255)

                                                                                                                                 ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGtt                     >  1:4299524/1‑150 (MQ=255)

                                                                                                                                 ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGtt                     >  1:4075376/1‑150 (MQ=255)

                                                                                                                                 ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGtt                     >  1:3586374/1‑150 (MQ=255)

                                                                                                                                 ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGtt                     >  1:821944/1‑150 (MQ=255)

                                                                                                                                 ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGtt                     >  1:1873468/1‑150 (MQ=255)

                                                                                                                                   cGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTAc                   <  2:781647/150‑1 (MQ=255)

                                                                                                                                   cGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTAc                   <  2:6837111/150‑1 (MQ=255)

                                                                                                                                    gAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTAcc                  >  2:7213775/1‑150 (MQ=255)

                                                                                                                                       gAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGaa               <  1:5660506/150‑1 (MQ=255)

                                                                                                                                       gAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGaa               <  1:5834639/150‑1 (MQ=255)

                                                                                                                                       gAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGaa               <  1:9136033/150‑1 (MQ=255)

                                                                                                                                          ggCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGAt            >  2:6999246/1‑150 (MQ=255)

                                                                                                                                            cGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCa          <  1:5222560/150‑1 (MQ=255)

                                                                                                                                            cGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCa          <  1:1376604/150‑1 (MQ=255)

                                                                                                                                             gtgtTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAg         <  1:8083037/150‑1 (MQ=255)

                                                                                                                                                ttGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCaa      <  1:3013020/150‑1 (MQ=255)

                                                                                                                                                  gcgcCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGt    >  1:5449251/1‑150 (MQ=255)

CTGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCC  >  NC_002947/5189805‑5190102

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