Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,826,287	Δ1 bp	6.1%	coding (4903/12954 nt)	pvdL ←	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,826,287	0	C	.	6.1%	208.6 / 6.5	66	coding (4903/12954 nt)	pvdL	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base . (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

GAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATG  >  NC_002947/4826141‑4826429                                                                                                                                                    |                                                                                                                                                   gAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCcg                                                                                                                                                  >  5:282960/1‑149 (MQ=255)       ggCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCt                                                                                                                                            >  6:235375/1‑149 (MQ=255)        gCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCtg                                                                                                                                           <  1:288215/149‑1 (MQ=255)              aGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGCCCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGttgttg                                                                                                                                     <  8:219496/149‑1 (MQ=255)                  gCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGc                                                                                                                                 >  7:105504/1‑149 (MQ=255)                     gCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCg                                                                                                                              <  5:235375/149‑1 (MQ=255)                      cTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCgg                                                                                                                             >  3:151666/1‑149 (MQ=255)                      cTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCgg                                                                                                                             >  4:174042/1‑149 (MQ=255)                         gACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGTCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCTGTTc                                                                                                                          >  6:20492/1‑149 (MQ=255)                         gACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGtt                                                                                                                           <  3:151827/148‑1 (MQ=255)                          aCATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCg                                                                                                                         >  1:78788/1‑149 (MQ=255)                            aTCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTg                                                                                                                       <  4:346200/149‑1 (MQ=255)                             tCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTg‑‑‑g                                                                                                                   <  3:174042/149‑1 (MQ=255)                             tCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTg‑‑‑g                                                                                                                   <  8:316659/149‑1 (MQ=255)                             tCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATACAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTg‑‑‑g                                                                                                                   >  2:199038/1‑149 (MQ=255)                               gCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCt                                                                                                                 <  8:195632/149‑1 (MQ=255)                                 aGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCg                                                                                                               <  2:78788/149‑1 (MQ=255)                                       ggCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGc                                                                                                         >  1:240094/1‑149 (MQ=255)                                           ggTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCAc                                                                                                     >  6:185332/1‑149 (MQ=255)                                             tGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACgc                                                                                                   >  7:116212/1‑149 (MQ=255)                                                 cgAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACgctg‑ct                                                                                              <  5:204516/149‑1 (MQ=255)                                                 cgAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACgctg‑ct                                                                                              <  5:185332/149‑1 (MQ=255)                                                  gAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGctg‑ctg                                                                                             >  7:133259/1‑149 (MQ=255)                                                       gATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GATCGACATGCTCACGCTG‑CTGcccaa                                                                                        >  7:196171/1‑146 (MQ=255)                                                           tGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTgg                                                                                    <  4:133885/149‑1 (MQ=255)                                                             ggTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑Gc                                                                                                                  <  3:324515/118‑1 (MQ=255)                                                             ggTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑Gc                                                                                                                  >  4:324515/1‑118 (MQ=255)                                                             ggTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCacac                                                                                       <  8:116212/144‑1 (MQ=255)                                                                   cAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTGCCTGCCCCAGTGTggcagca                                                                             <  8:105504/149‑3 (MQ=255)                                                                   cAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGggcaggca                                                                            <  6:282960/149‑1 (MQ=255)                                                                      cAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCcgc                                                                                                                                                 >  6:294958/1‑81 (MQ=255)                                                                      cAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGAAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTGCCTGCCCCAGTGTggcagcagcg                                                                          <  8:245886/149‑6 (MQ=255)                                                                       aaCGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGc                                                                        <  8:49838/149‑1 (MQ=255)                                                                            gTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGc                                                                   >  6:258599/1‑149 (MQ=255)                                                                                  accGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGa                                                             <  6:206403/149‑1 (MQ=255)                                                                                      aGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCtt                                                         >  3:270319/1‑149 (MQ=255)                                                                                        ggCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGc                                                       >  2:46703/1‑149 (MQ=255)                                                                                           gaagaaGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGc                                                    >  7:33904/1‑149 (MQ=255)                                                                                            aagaagTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCt                                                   >  4:352490/1‑149 (MQ=255)                                                                                            aagaagTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGGCGACCTTGCCGCt                                                   <  1:199038/149‑1 (MQ=255)                                                                                              gaagTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTgg                                                 >  1:31324/1‑149 (MQ=255)                                                                                                     tcCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCaa                                          >  6:353677/1‑149 (MQ=255)                                                                                                       ctgctgGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATg                                        <  4:270319/149‑1 (MQ=255)                                                                                                        tGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATgg                                       >  4:305278/1‑149 (MQ=255)                                                                                                         gCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATgg                                       >  8:202900/1‑148 (MQ=255)                                                                                                                aaTGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCtgtg                               >  6:86545/1‑149 (MQ=255)                                                                                                                aaTGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGAGGATGGTCGACATCCTCACGCTC‑CTGCCACACTCGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGCCATCt                                  >  1:289846/1‑149 (MQ=255)                                                                                                                    cgcgGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCa                            >  4:310706/1‑148 (MQ=255)                                                                                                                      cgGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAg                         <  2:31324/149‑1 (MQ=255)                                                                                                                          aGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGAt                     <  1:46703/149‑1 (MQ=255)                                                                                                                          aGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGAt                     <  5:176544/149‑1 (MQ=255)                                                                                                                             ttCTACACGTAGC‑GCCAGATTGCCGCGACCTGTTGTTGCAGCTCGTTTCGTT‑‑‑GCTCGACATGCTCACGCTG‑CTGACACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAATGGATGa                   <  8:196171/148‑1 (MQ=255)                                                                                                                             ttCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATc                                   <  1:300311/132‑1 (MQ=255)                                                                                                                             ttCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATc                                   >  2:300311/1‑132 (MQ=255)                                                                                                                             ttCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAg                  <  3:310706/149‑1 (MQ=255)                                                                                                                             ttCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAg                  >  5:354830/1‑149 (MQ=255)                                                                                                                              acAACAGGTCGCGGCCATCTGG‑CGCGACGTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTggg                                                >  3:152443/2‑118 (MQ=37)                                                                                                                              acAACAGGTCGCGGCCATCTGG‑CGCGACGTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTggg                                                <  4:152443/117‑1 (MQ=37)                                                                                                                               cAACAGGTCGCGGCCATCTGG‑CGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCt                                  <  7:164505/131‑1 (MQ=255)                                                                                                                               cAACAGGTCGCGGCCATCTGG‑CGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCt                                  >  8:164505/1‑131 (MQ=255)                                                                                                                                aaCACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGtt                <  4:151666/148‑1 (MQ=255)                                                                                                                                aaCACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTg               >  2:294309/1‑149 (MQ=255)                                                                                                                                   acGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTgcgc            <  2:240094/149‑1 (MQ=255)                                                                                                                                        gc‑gcCAGATGGCGGCGACCTGTTGGTGCAGCGCGGTGCGTG‑‑‑GCTCGCCTTGCTCTCGGTG‑GTGCCACACTGGGGCTGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCa       <  2:289846/149‑1 (MQ=255)                                                                                                                                             cAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCAt   >  1:341388/1‑149 (MQ=255)                                                                                                                                              aGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATg  >  7:3573/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGGCCAATGCGCGGCAGATTCAACACGTCGC‑GCCAGATGGCCGCGACCTGTTGTTGCAGCTCGGTTCGTG‑‑‑GCTCGACATGCTCACGCTG‑CTGCCACACTGGGGCAGGCAGTGCCTTGCGGTCGACCTTGCCGCTGGGGCCCAATGGCATCTGTGCCAAAGGGATGAGTTGCGCCGGCACCATG  >  NC_002947/4826141‑4826429

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