Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,660,182	G→T	11.1%	intergenic (+35/‑194)	metK → / → ligB	methionine adenosyltransferase/DNA ligase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,660,182	0	G	T	11.1%	178.3 / 12.1	63	intergenic (+35/‑194)	metK/ligB	methionine adenosyltransferase/DNA ligase B
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base T (4/3);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-01

CACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCC  >  NC_002947/5660037‑5660321                                                                                                                                                  |                                                                                                                                             cACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGa                                                                                                                                           <  6:405628/149‑1 (MQ=255)  aCTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCgct                                                                                                                                          >  5:358623/1‑147 (MQ=255)      cGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCgg                                                                                                                                       >  6:176724/1‑148 (MQ=255)       ggCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGGCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCgg                                                                                                                                       >  1:408833/1‑147 (MQ=255)         ccGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAgg                                                                                                                                   <  6:358623/149‑1 (MQ=255)            tACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAGGGGc                                                                                                                                <  2:280214/149‑1 (MQ=255)                 cGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGAt                                                                                                                                          >  4:359066/1‑134 (MQ=255)                 cGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGAt                                                                                                                                          <  3:359066/134‑1 (MQ=255)                  gcagcaGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCtttttt                                                                                                                          >  6:2886/1‑149 (MQ=255)                  gcagcaGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCtttttt                                                                                                                          >  3:44817/1‑149 (MQ=255)                    cgcagACGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGCGCGCACCGCCCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtgt                                                                                                                        <  2:408833/148‑1 (MQ=255)                                tCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACgg                                                                                                            >  2:176984/1‑148 (MQ=255)                                tCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGc                                                                                                           <  7:18808/149‑1 (MQ=255)                                   gCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCAtt                                                                                                        >  3:255781/1‑149 (MQ=255)                                    cgacgaCACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCAttt                                                                                                       <  5:2886/149‑1 (MQ=255)                                       cgacACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCa                                                                                                    >  3:396457/1‑149 (MQ=255)                                        gacACCTTCACCACCTTCACCTGGGAGCGCTCCGACCGCGCCGAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCTGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCa                                                                                                    <  6:277777/148‑1 (MQ=255)                                        gacACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCAt                                                                                                   >  1:231972/1‑149 (MQ=255)                                        gacACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGCGC‑CTGTTACGGCATTTTCAt                                                                                                   >  7:232420/1‑149 (MQ=255)                                          caccttcaccaccttcacCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCg                                                                                                 >  8:46935/1‑149 (MQ=255)                                           accttcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAGGGGc                                                                                                                                >  5:368539/1‑118 (MQ=255)                                           accttcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAGGGGc                                                                                                                                <  6:368539/118‑1 (MQ=255)                                             cttcaccaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGccc                                                                                              <  5:176724/149‑1 (MQ=255)                                                 accaccttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTc                                                                                           <  5:78895/148‑1 (MQ=255)                                                      cttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATc                                                                                      >  7:152765/1‑148 (MQ=255)                                                       ttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCa                                                                                                    >  6:143118/1‑133 (MQ=255)                                                       ttcaccTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCa                                                                                                    <  5:143118/133‑1 (MQ=255)                                                                  gcgcACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGa                                                                                    <  3:24425/138‑1 (MQ=255)                                                                  gcgcACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGa                                                                                    >  4:24425/1‑138 (MQ=255)                                                                     cACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGtt                                                                      >  3:121754/1‑149 (MQ=255)                                                                     cACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGtt                                                                      >  4:179717/1‑149 (MQ=255)                                                                      accgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGttt                                                                     <  6:259010/149‑1 (MQ=255)                                                                      accgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGttt                                                                     >  1:323282/1‑149 (MQ=255)                                                                      accgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGttt                                                                     >  5:259010/1‑149 (MQ=255)                                                                      accgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGttt                                                                     >  5:21314/1‑149 (MQ=255)                                                                            cgcgCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCt                                                               <  4:396457/149‑1 (MQ=255)                                                                            cgcgCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCt                                                               <  2:323282/149‑1 (MQ=255)                                                                              cgcCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCt                                                               <  4:4753/147‑1 (MQ=255)                                                                              cgcCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCt                                                               >  3:4753/1‑147 (MQ=255)                                                                                     tCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGACCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCgg                                                      >  2:338996/1‑149 (MQ=255)                                                                                          gCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCtt                                                  <  3:103154/148‑1 (MQ=255)                                                                                           cGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTgg                                                <  8:152765/149‑1 (MQ=255)                                                                                           cGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTg                                                 <  2:207567/148‑1 (MQ=255)                                                                                           cGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGCGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTg                                                 <  8:232420/148‑1 (MQ=255)                                                                                               aCGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgaga                                             <  3:179717/148‑1 (MQ=255)                                                                                                 gCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGa                                          <  7:46935/149‑1 (MQ=255)                                                                                                   tGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCtt                                                                  >  1:198853/1‑123 (MQ=255)                                                                                                   tGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCtt                                                                  <  2:198853/123‑1 (MQ=255)                                                                                                      cGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCAc                                      >  8:269400/1‑148 (MQ=255)                                                                                                       ggCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCa                                    <  1:176984/149‑1 (MQ=255)                                                                                                                    ccTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCa                         <  4:255781/147‑1 (MQ=255)                                                                                                                     cTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCAcc                       <  4:44817/148‑1 (MQ=255)                                                                                                                       gCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgaga                                             >  6:217184/1‑124 (MQ=255)                                                                                                                       gCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgaga                                             <  5:217184/124‑1 (MQ=255)                                                                                                                        cAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGa                                          >  4:223114/1‑126 (MQ=255)                                                                                                                        cAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGa                                          <  3:223114/126‑1 (MQ=255)                                                                                                                             gcacaaaaaaGCCCCTGCCGATCGCTCGGCAGGGGCTTTTTTGTGCGCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACt               >  1:294074/8‑149 (MQ=38)                                                                                                                               aGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgca            >  7:43940/1‑149 (MQ=255)                                                                                                                                 aaaaaaGCCCCTGCCGATCGCTCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgcacg          >  5:352211/4‑149 (MQ=255)                                                                                                                                   gAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGTCACGAGCACCAGCGAAATATCACCTGGCAACTgcacgcc        >  1:405468/1‑148 (MQ=255)                                                                                                                                         cccTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgag                                              <  5:418769/105‑1 (MQ=255)                                                                                                                                         cccTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgag                                              >  6:418769/1‑105 (MQ=255)                                                                                                                                         cccTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCAccgccc  >  1:26621/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                             CACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGAGCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGC‑CTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCACGCACCGCCC  >  NC_002947/5660037‑5660321

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