Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,948,667:1	+A	5.6%	coding (181/435 nt)	fliO ←	flagellar protein FliO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,948,667	1	.	A	5.6%	292.5 / 11.4	72	coding (181/435 nt)	fliO	flagellar protein FliO
Reads supporting (aligned to +/- strand):  ref base . (35/33);  new base A (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CGTGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGTGGCTG  >  NC_002947/4948520‑4948812                                                                                                                                                     |                                                                                                                                                  cGTGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAg                                                                                                                                                    <  8:277154/148‑1 (MQ=255)     cAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑Gcgcc                                                                                                                                              >  2:312503/1‑149 (MQ=255)        ggCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGCGC‑CTGCAGAGCACCgca                                                                                                                                           <  8:195060/149‑1 (MQ=255)        ggCCTAGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAC‑GCGCCgca                                                                                                                                           >  3:45043/1‑149 (MQ=255)         gCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGCGC‑CTGCAGAGCACCgcag                                                                                                                                          >  5:166326/1‑149 (MQ=255)         gCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGCGC‑CTGCAGAGCACCgcag                                                                                                                                          <  6:166326/149‑1 (MQ=255)          ccTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCgcagc                                                                                                                                         <  1:299056/149‑1 (MQ=255)           cTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCgcag                                                                                                                                          >  5:313630/1‑147 (MQ=255)           cTCGATGGGGCCGGGGGTATGGCCGATCAGGATGTGCTCCTTGCCCTCCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCgcagca                                                                                                                                        <  4:45043/149‑1 (MQ=255)             cGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCagcagc                                                                                                                                      >  5:282042/1‑149 (MQ=255)                  gTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAgccaggc                                                                                                                                 <  6:313630/149‑1 (MQ=255)                           gTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCagaaagaa                                                                                                                        <  2:51258/149‑1 (MQ=255)                           gTATGGCCGATCAGCCTCTCCGCCGTGCCCACCCCCCCCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAAGCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCagaaagaa                                                                                                                        <  4:263413/149‑1 (MQ=255)                                gCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGAt                                                                                                                     >  6:174757/1‑147 (MQ=255)                                  cGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAgg                                                                                                                 <  5:67827/149‑1 (MQ=255)                                         gATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCccacca                                                                                                          >  5:203641/1‑149 (MQ=255)                                               cTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCccaccac                                                                                                         <  8:20027/144‑1 (MQ=255)                                                  cTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACcagcag                                                                                                   <  1:224195/147‑1 (MQ=255)                                                     gCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCGTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAgg                                                                                               <  4:241974/148‑1 (MQ=255)                                                      cccACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCCCCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGcc                                                                                             >  2:106903/1‑149 (MQ=255)                                                       ccACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGccc                                                                                            <  2:120555/149‑1 (MQ=255)                                                         aCCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGccc                                                                                            >  7:148872/1‑147 (MQ=255)                                                           cTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGa                                                                                        <  8:258221/149‑1 (MQ=255)                                                             gCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAgg                                                                                               <  2:103977/140‑1 (MQ=255)                                                             gCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAgg                                                                                               >  1:103977/1‑140 (MQ=255)                                                               aCCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACATGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCACCAGCAGGCCCAAGACCa                                                                                     >  5:239565/1‑148 (MQ=255)                                                                ccAGCAGCAGGCGGTCACGCGGGCCGATGGCTCGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATc                                                                                   <  4:27555/149‑1 (MQ=255)                                                                 cagcagcagGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGAAGCAGGCCCAAGACCATCt                                                                                  >  7:252710/1‑149 (MQ=255)                                                                   gcagcagGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTgc                                                                                >  8:95079/1‑149 (MQ=255)                                                                    cagcagGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTgcg                                                                               >  8:310843/1‑149 (MQ=255)                                                                     agcagGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTgcgc                                                                              >  7:623/1‑149 (MQ=255)                                                                         ggCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGc                                                                                              <  2:310922/129‑1 (MQ=255)                                                                         ggCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGc                                                                                              >  1:310922/1‑129 (MQ=255)                                                                            ggTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGc                                                                       >  2:199964/1‑149 (MQ=255)                                                                                    gggCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGcc                                                                  >  8:241584/1‑146 (MQ=255)                                                                                     ggCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGa                                                              <  6:203641/149‑1 (MQ=255)                                                                                               aCGGCTGCCGACAATCTCGATCACCTGCCCCCCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGtt                                                    >  3:52544/1‑149 (MQ=255)                                                                                                 ggCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGtt                                                    >  7:16294/1‑147 (MQ=255)                                                                                                 ggCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGa                                                  <  8:623/149‑1 (MQ=255)                                                                                                  gCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAg                                                 <  8:41197/149‑1 (MQ=255)                                                                                                  gCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAg                                                 <  5:174757/149‑1 (MQ=255)                                                                                                      ccGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAgccgc                                             >  3:173588/1‑149 (MQ=255)                                                                                                      ccGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAgccgc                                             <  7:95079/149‑1 (MQ=255)                                                                                                      ccGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAgccgc                                             >  5:41331/1‑149 (MQ=255)                                                                                                       cGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAgccgcc                                            >  5:51507/1‑149 (MQ=255)                                                                                                        gACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGccgccg                                           <  8:252710/149‑1 (MQ=255)                                                                                                          cAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCggg                                         <  6:51507/149‑1 (MQ=255)                                                                                                          cAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCgg                                          <  8:148872/148‑1 (MQ=255)                                                                                                             tctcGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTg                                      >  7:256285/1‑149 (MQ=255)                                                                                                              ctcGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCACCCAGGCCAGAAAGAAGGTCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCACCGGGGTgg                                     >  3:261733/1‑149 (MQ=255)                                                                                                                 gATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGgctg                                  <  1:32080/149‑1 (MQ=255)                                                                                                                     aCCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTgg                              >  4:40521/1‑149 (MQ=255)                                                                                                                          cccACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGa                                                  <  6:244430/124‑1 (MQ=255)                                                                                                                          cccACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGa                                                  >  5:244430/1‑124 (MQ=255)                                                                                                                              cccTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGc                     <  7:241584/149‑1 (MQ=255)                                                                                                                               ccTTCACTGCGGTGCTCTGCAG‑GCGCCGCGGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCg                      >  1:108197/1‑147 (MQ=255)                                                                                                                               ccTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCg                      >  5:120615/1‑147 (MQ=255)                                                                                                                                cTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAg                   >  4:69616/1‑149 (MQ=255)                                                                                                                                   tgcTGCGGCGC‑CTGCAGAGCACCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCACCAAGACTGCCAGGGTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGc                  >  7:62010/3‑147 (MQ=255)                                                                                                                                   cACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGcc                 <  6:41331/148‑1 (MQ=255)                                                                                                                                    aCTGCGGTGCTCTGCAG‑GCGCCGCTGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCca                <  2:108197/148‑1 (MQ=255)                                                                                                                                    aCTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCcag               <  4:52544/149‑1 (MQ=255)                                                                                                                                           tGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCg        >  8:218272/1‑149 (MQ=255)                                                                                                                                             ctctGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGtg      >  5:92623/1‑149 (MQ=255)                                                                                                                                             ctctGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGtg      >  4:214076/1‑149 (MQ=255)                                                                                                                                             ctctGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCAGCGGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGtg      >  1:108317/1‑149 (MQ=255)                                                                                                                                                tGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCaa                                                                >  5:152940/1‑88 (MQ=255)                                                                                                                                                tGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCaa                                                                <  6:152940/88‑1 (MQ=255)                                                                                                                                                tGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGtgg     <  4:173588/147‑1 (MQ=255)                                                                                                                                                  cAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCGGCCAGGCAGGCGtggct   <  8:256285/147‑1 (MQ=255)                                                                                                                                                  cAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGtggc    <  1:199964/146‑1 (MQ=255)                                                                                                                                                   aG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGtggctg  <  2:332062/147‑1 (MQ=255)                                                                                                                                                   aG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGtggctg  >  1:332062/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CGTGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGTGGCTG  >  NC_002947/4948520‑4948812

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