Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,948,661	Δ1 bp	6.3%	coding (187/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,661	0	T	.	6.3%	196.5 / 7.4	63	coding (187/435 nt)	fliO	flagellar protein FliO
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base . (2/2);  total (32/31)
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.48e-01

CCAGTACGTGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGTG  >  NC_002947/4948514‑4948808                                                                                                                                                    |                                                                                                                                                     ccAGTACGTGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGctc                                                                                                                                                     >  1:51258/1‑149 (MQ=255)       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                                                                                           >  1:103977/1‑140 (MQ=255)                                                                   gCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAgg                                                                                           <  2:103977/140‑1 (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                                                                                          >  1:310922/1‑129 (MQ=255)                                                                               ggCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGc                                                                                          <  2:310922/129‑1 (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                                             <  5:174757/149‑1 (MQ=255)                                                                                                        gCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAg                                             <  8:41197/149‑1 (MQ=255)                                                                                                            ccGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAgccgc                                         >  3:173588/1‑149 (MQ=255)                                                                                                            ccGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAgccgc                                         >  5:41331/1‑149 (MQ=255)                                                                                                            ccGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAgccgc                                         <  7:95079/149‑1 (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                                              >  5:244430/1‑124 (MQ=255)                                                                                                                                cccACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGa                                              <  6:244430/124‑1 (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)                                                                                                                                                    |                                                                                                                                                     CCAGTACGTGCAGGGCCTCGATGGTGCCGGGGGTATGGCCGATCAGGATCTGCTCCTTGCCCACCTGCACCAGCAGCAGGCGGTCACGCGGGCCGATGGCACGGCTGCCGACAATCTCGATCACCTGCCCACCCTTCACTGCGGTGCTCTGCAG‑GCGCCGCAGCAGCCAGGCCAGAAAGAAGATCAGGCCCACCACCAGCAGCAGGCCCAAGACCATCTGCGCCAGCTGCCCGCCAAGACTGCCAGGTTGAGCCGCCGGGGTGGCTGCTGGCGTGGCCGCAGCCAGGCAGGCGTG  >  NC_002947/4948514‑4948808

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