Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,618,111	G→C	8.8%	R181G (CGT→GGT)	wzy ←	O‑antigen polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,618,111	0	G	C	8.8%	206.8 / 9.2	68	R181G (CGT→GGT)	wzy	O‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base C (2/4);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

AGAACACCCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGA  >  NC_002947/5617964‑5618252                                                                                                                                                     |                                                                                                                                               agaACACCCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGt                                                                                                                                               <  1:128217/149‑1 (MQ=255) agaACACCCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGt                                                                                                                                               >  4:246199/1‑149 (MQ=255)        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                                                                                                                                   <  7:82094/115‑1 (MQ=255)                                               cacGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACCTCGGCGACGTGGc                                                                                                                                   >  8:82094/1‑115 (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                                                                                                                    >  2:220410/1‑111 (MQ=255)                                                                 cAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCaa                                                                                                                    <  1:220410/111‑1 (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                                                               <  7:73290/143‑1 (MQ=255)                                                                                      gCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGc                                                               >  8:73290/1‑143 (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                                                                                                  >  8:151589/1‑88 (MQ=255)                                                                                                           ccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAAt                                                                                                  <  7:151589/88‑1 (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        <  8:235173/149‑1 (MQ=255)                                                                                                                                        aacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtgat        <  6:41053/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)                                                                                                                                                     |                                                                                                                                               AGAACACCCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGG‑CCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAG‑CATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGA  >  NC_002947/5617964‑5618252

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