Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,660,185	A→C	7.7%	intergenic (+38/‑191)	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,185	0	A	C	7.7%	199.1 / 6.4	65	intergenic (+38/‑191)	metK/ligB	methionine adenosyltransferase/DNA ligase B
Reads supporting (aligned to +/- strand):  ref base A (21/39);  new base C (2/3);  total (23/42)
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.38e-01

GCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCA  >  NC_002947/5660044‑5660311                                                                                                                                                |                                                                                                                               gCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAg                                                                                                                         <  7:205217/149‑1 (MQ=255)  ccGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAgg                                                                                                                        <  6:134665/149‑1 (MQ=255)  ccGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAgg                                                                                                                        <  6:312691/149‑1 (MQ=255)  ccGTACTCCGCAGCAGAAGACTGTCGACGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAgg                                                                                                                        >  1:168823/1‑149 (MQ=255)   cGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCGCCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAggg                                                                                                                       <  1:296869/149‑1 (MQ=255)   cGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAggg                                                                                                                       <  4:4119/149‑1 (MQ=255)    atcgggtgcagcaGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAgggg                                                                                                                      <  8:166569/142‑1 (MQ=255)          cGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCttttt                                                                                                                >  3:101347/1‑149 (MQ=255)             agcagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtgt                                                                                                             <  2:133259/149‑1 (MQ=255)             agcagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtgt                                                                                                             <  7:95544/149‑1 (MQ=255)               cagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGc                                                                                                           >  6:309194/1‑149 (MQ=255)               cagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGc                                                                                                           <  1:48306/149‑1 (MQ=255)               cagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGc                                                                                                           <  3:241335/149‑1 (MQ=255)               cagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACAAAAGCCCCTGCCGAGCGATCGGCAGGGGCttttgtctgg                                                                                                           <  5:325206/149‑7 (MQ=255)               cagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACAAAAGCCCCTGCCGAGCGATCGGCAGGGGCttttgtctgg                                                                                                           <  3:7776/149‑7 (MQ=255)                 gaagaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCTTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCt                                                                                                         <  2:81524/149‑1 (MQ=255)                   agaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGTA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGt                                                                                                       <  2:205981/149‑1 (MQ=255)                   agaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTg                                                                                                        <  2:321744/148‑1 (MQ=255)                    gaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCtttttt                                                                                                               <  7:269709/140‑1 (MQ=255)                    gaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCtttttt                                                                                                               >  8:269709/1‑140 (MQ=255)                    gaCTGGCGGCGACGTCACCTTCACCACCTTCACCTGGGATCACGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTg                                                                                                        <  1:272154/149‑1 (MQ=255)                       tGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAGGGGCtttcgtctg                                                                                                            >  4:8250/1‑134 (MQ=255)                       tGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAGGGGCtttcgtctg                                                                                                            <  3:8250/140‑7 (MQ=255)                         tCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtgt                                                                                                             >  7:86546/1‑137 (MQ=255)                         tCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtgt                                                                                                             <  8:86546/137‑1 (MQ=255)                            gCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCAtt                                                                                              <  8:65050/149‑1 (MQ=255)                                cgacACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAgggg                                                                                                                      >  7:97421/1‑121 (MQ=255)                                cgacACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAgggg                                                                                                                      <  8:97421/121‑1 (MQ=255)                                 gacACCTTCACCACGTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCAt                                                                                         >  7:273500/1‑149 (MQ=255)                                   caccttcaccaccttcacCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCg                                                                                       <  2:97211/149‑1 (MQ=255)                                   caccttcaccaccttcacCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCg                                                                                       <  5:82173/149‑1 (MQ=255)                                           ccaccttcaccTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGACCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGt                                                                                  >  1:248571/1‑146 (MQ=255)                                            caccttcaccTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGt                                                                                  >  8:274478/1‑145 (MQ=255)                                             accttcaccTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGAt                                                                             <  5:45536/149‑1 (MQ=255)                                              ccttcaccTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCCTCGCCCGTCAGATc                                                                            <  7:274478/149‑1 (MQ=255)                                                  caccTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAgc                                                                        >  8:303958/1‑149 (MQ=255)                                                                cACCGACCGCGCCCGGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGtt                                                            >  8:193408/1‑149 (MQ=255)                                                                cACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGtt                                                            >  5:57215/1‑149 (MQ=255)                                                                   cgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCt                                                         <  6:282101/149‑1 (MQ=255)                                                                      ccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGc                                                      <  2:168823/149‑1 (MQ=255)                                                                      ccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGc                                                      <  3:203146/149‑1 (MQ=255)                                                                        gcgcCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTc                                                    >  2:276506/1‑149 (MQ=255)                                                                           cccAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTcgtc                                                 >  3:26336/1‑146 (MQ=255)                                                                                  gTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCggg                                           <  3:226960/148‑1 (MQ=255)                                                                                   ttGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCt                                         <  8:273500/149‑1 (MQ=255)                                                                                      cGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTg                                       <  7:193408/148‑1 (MQ=255)                                                                                      cGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTg                                       <  4:108139/148‑1 (MQ=255)                                                                                      cGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTg                                       >  2:269817/1‑148 (MQ=255)                                                                                        tGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgag                                    <  1:276506/149‑1 (MQ=255)                                                                                         gACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgaga                                   <  7:303958/149‑1 (MQ=255)                                                                                           cGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACg                                 <  6:57215/149‑1 (MQ=255)                                                                                             cTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAg                               <  2:38314/149‑1 (MQ=255)                                                                                                ccGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCAc                            >  7:298592/1‑149 (MQ=255)                                                                                                   gCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAg                         >  7:204621/1‑149 (MQ=255)                                                                                                       gTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGa                      >  5:83616/1‑148 (MQ=255)                                                                                                        tAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGaa                     <  8:298592/148‑1 (MQ=255)                                                                                                               ccTGCAGCTCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCAc              <  5:309194/148‑1 (MQ=255)                                                                                                               ccTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCa               <  8:204621/147‑1 (MQ=255)                                                                                                                cTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCt            <  1:269817/149‑1 (MQ=255)                                                                                                                cTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCt            >  6:190432/1‑149 (MQ=255)                                                                                                                  gCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCg                                             >  6:320526/1‑114 (MQ=255)                                                                                                                   cAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGa          <  4:101347/148‑1 (MQ=255)                                                                                                                   cAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGa          >  2:292464/1‑148 (MQ=255)                                                                                                                         cAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgc   <  2:248571/149‑1 (MQ=255)                                                                                                                          aGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgca  >  1:261539/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                               GCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCA  >  NC_002947/5660044‑5660311

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