Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,618,115	C→G	8.8%	R179R (CGG→CGC)	wzy ←	O‑antigen polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,618,115	0	C	G	8.8%	204.5 / 9.9	68	R179R (CGG→CGC)	wzy	O‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base C (28/34);  new base G (2/4);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

CCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAG  >  NC_002947/5617971‑5618263                                                                                                                                                  |                                                                                                                                                      ccGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAg                                                                                                                                                   >  7:54934/1‑149 (MQ=255)    gCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCCAGGCAAGCCAGCGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGt                                                                                                                                                 <  1:250313/149‑1 (MQ=255)     cTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTgg                                                                                                                                               >  5:56859/1‑149 (MQ=255)           cGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTa                                                                                                                                         <  3:86466/149‑1 (MQ=255)           cGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACCTCGGCGACGTGGCttgca                                                                                                                                         <  7:141027/149‑5 (MQ=255)            ggTCACGGAACCACAACCGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACCTCGGCGACGTGGCttgcag                                                                                                                                        <  7:245893/149‑6 (MQ=255)               cACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTagag                                                                                                                                      >  4:307679/1‑148 (MQ=255)                 cGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CAt                                                                                                                                  <  1:228240/149‑1 (MQ=255)                        acaACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGc                                                                                                                           <  5:230714/149‑1 (MQ=255)                                gcCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCg                                                                                                                    >  7:57620/1‑148 (MQ=255)                                gcCATCACCACGGTAATCACCAGTGCCAGAACCCCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGa                                                                                                                   <  6:112027/149‑1 (MQ=255)                                  cATCTCCACGGTAATCTCCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGGATGGCTAAGG‑CCCCCAGTCAAGCCAGTGCAGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGATGAAGGGCGAAc                                                                                                                 <  6:273661/149‑1 (MQ=255)                                        cacGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGa                                                                                                           >  7:122478/1‑149 (MQ=255)                                        cacGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACCTCGGCGACGTGGc                                                                                                                                              >  8:82094/1‑115 (MQ=255)                                        cacGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACCTCGGCGACGTGGc                                                                                                                                              <  7:82094/115‑1 (MQ=255)                                         acGGTAATCACCAGTGCCAGAACCGCCCCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCACGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGa                                                                                                           >  6:33303/1‑148 (MQ=255)                                            gTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACAt                                                                                                        >  2:39232/1‑148 (MQ=255)                                                tCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCCATCAAGAGCAAGGCCGAACCAATGACATAAGc                                                                                                    <  1:39232/149‑1 (MQ=255)                                                  aCCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGccc                                                                                                  <  2:167981/148‑1 (MQ=255)                                                    cAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCa                                                                                               <  3:246199/149‑1 (MQ=255)                                                          cAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCaa                                                                                                                               <  1:220410/111‑1 (MQ=255)                                                          cAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCaa                                                                                                                               >  2:220410/1‑111 (MQ=255)                                                           aGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATggg                                                                                        >  1:77129/1‑149 (MQ=255)                                                               ccGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTgg                                                                                    <  6:56859/149‑1 (MQ=255)                                                                cGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTggg                                                                                   >  2:78195/1‑149 (MQ=255)                                                                cGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTACCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATCACATAAGCCCCCAGGATGGGGTggg                                                                                   >  3:9743/1‑149 (MQ=255)                                                                 gctccacGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCACGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGa                                                                                  <  5:33303/146‑1 (MQ=255)                                                                  caccacGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGat                                                                                 <  1:238416/149‑1 (MQ=255)                                                                  caccacGGCTCTGGCTGAGCATAGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGat                                                                                 <  4:138557/149‑1 (MQ=255)                                                                        ggCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCg                                                                           >  4:68911/1‑149 (MQ=255)                                                                          ctctGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGcc                                                                         <  3:307679/149‑1 (MQ=255)                                                                               gCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGc                                                                          >  8:73290/1‑143 (MQ=255)                                                                               gCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGc                                                                          <  7:73290/143‑1 (MQ=255)                                                                               gCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATAc                                                                    >  5:41053/1‑149 (MQ=255)                                                                                        ggCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGatat                                                                               <  2:160612/129‑1 (MQ=255)                                                                                        ggCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGatat                                                                               >  1:160612/1‑129 (MQ=255)                                                                                        ggCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGcc                                                           >  7:118528/1‑149 (MQ=255)                                                                                         gCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGcc                                                           <  8:96280/148‑1 (MQ=255)                                                                                         gCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCt                                                          >  6:51858/1‑149 (MQ=255)                                                                                            aaGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGaa                                                       >  1:169937/1‑149 (MQ=255)                                                                                            aaGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGaa                                                       >  1:252472/1‑149 (MQ=255)                                                                                                  ccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACAt                                                                                                        >  6:16076/1‑95 (MQ=255)                                                                                                  ccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACAt                                                                                                        <  5:16076/95‑1 (MQ=255)                                                                                                   cccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCa                                                  >  7:273325/1‑148 (MQ=255)                                                                                                    ccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAAt                                                                                                             <  7:151589/88‑1 (MQ=255)                                                                                                    ccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAAt                                                                                                             >  8:151589/1‑88 (MQ=255)                                                                                                       ggCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGgc                                             >  3:163730/1‑149 (MQ=255)                                                                                                          aaGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCa                                          >  7:235173/1‑149 (MQ=255)                                                                                                             ccAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGc                                       <  8:122478/149‑1 (MQ=255)                                                                                                                 tGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCa                                   <  7:95192/149‑1 (MQ=255)                                                                                                                      ccagccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCa                              <  8:54934/149‑1 (MQ=255)                                                                                                                        agccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATa                            <  3:68911/149‑1 (MQ=255)                                                                                                                           cacaACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGaa                         >  6:205189/1‑149 (MQ=255)                                                                                                                            acaacaACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAAc                        <  4:125892/149‑1 (MQ=255)                                                                                                                               acaacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtg                     <  5:51858/149‑1 (MQ=255)                                                                                                                                caacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtga                    >  8:10917/1‑149 (MQ=255)                                                                                                                                 aacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtgat                   <  6:41053/149‑1 (MQ=255)                                                                                                                                 aacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtgat                   <  8:235173/149‑1 (MQ=255)                                                                                                                                  acTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtgatg                  <  1:78195/149‑1 (MQ=255)                                                                                                                                   cTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCCCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCa                              <  3:23528/136‑1 (MQ=255)                                                                                                                                   cTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCCCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCa                              >  4:23528/1‑136 (MQ=255)                                                                                                                                    tacgaGCCACCTCGGCGACGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCAc                                  <  8:52104/127‑1 (MQ=255)                                                                                                                                    tacgaGCCACCTCGGCGACGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCAc                                  >  7:52104/5‑131 (MQ=255)                                                                                                                                       aaGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGga             <  2:77129/149‑1 (MQ=255)                                                                                                                                          ccACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATggg                                                                                        <  8:252646/71‑1 (MQ=255)                                                                                                                                          ccACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATggg                                                                                        >  7:252646/1‑71 (MQ=255)                                                                                                                                                 gCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCa   <  8:273325/149‑1 (MQ=255)                                                                                                                                                 gCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCa   <  7:10917/149‑1 (MQ=255)                                                                                                                                                  ccGAGGTTGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCATGATGGGGGGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAg  <  6:183009/149‑1 (MQ=255)                                                                                                                                                  ccGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAg  >  2:203835/1‑149 (MQ=255)                                                                                                                                                  ccGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAg  >  3:79128/1‑149 (MQ=255)                                                                                                                                                  ccGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAg  <  1:317954/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                      CCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAG  >  NC_002947/5617971‑5618263

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