Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,743,365	A→G	12.5%	intergenic (‑198/+49)	fbp ← / ← glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,743,365	0	A	G	12.5%	141.0 / 10.0	57	intergenic (‑198/+49)	fbp/glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base G (3/4);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

GAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCA  >  NC_002947/5743218‑5743512                                                                                                                                                     |                                                                                                                                                    gagaGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGCTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGaa                                                                                                                                                     <  1:294323/149‑1 (MQ=255)  agagTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTCTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGc                                                                                                                                                   <  1:98537/149‑1 (MQ=255)        aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTTGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCg                                                                                                                                              <  3:79909/148‑1 (MQ=255)        aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                                               >  6:6817/1‑147 (MQ=255)         cGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCTTCAAGGCTTGCgct                                                                                                                                            <  8:88641/149‑3 (MQ=255)          gggCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCgaga                                                                                                                                           >  7:33616/1‑149 (MQ=255)           ggCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCgaga                                                                                                                                           >  5:240788/1‑148 (MQ=255)            gCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGaca                                                                                                                                         >  4:115733/1‑149 (MQ=255)                gCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTGCAGCCTTCGCCAGGCCGGTGGGCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCt                                                                                                                                      <  7:266321/148‑1 (MQ=255)                gCGAGTGATGAGACTCGCCACGGATCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGTTCTGG‑TTGTGTGTCCATTCGCGTGCGCGCTCGCTCCCACAGGTACCTCGCAAGCCATGAAGCTTGCGAGACACCtg                                                                                                                                     >  4:144831/1‑149 (MQ=255)                  gAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgtg                                                                                                                                   >  8:80801/1‑149 (MQ=255)                   aGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTgg                                                                                                                                  <  4:215613/149‑1 (MQ=255)                     tgatgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGa                                                                                                                                <  8:283556/149‑1 (MQ=255)                          agaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCgg                                                                                                                            <  8:33616/148‑1 (MQ=255)                          agaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCgctgtac                                                                                                                                        >  1:26497/1‑130 (MQ=21)                          agaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCgctgtac                                                                                                                                        <  2:26497/136‑7 (MQ=21)                          agaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCACCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGAGCga                                                                                                                            <  2:35124/148‑2 (MQ=18)                                gCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGcc                                                                                                                     >  7:6329/1‑149 (MQ=255)                                gCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑GTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGc                                                                                                                      <  4:235914/148‑1 (MQ=255)                                 ccACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGt                                                                                                                        <  8:119539/145‑1 (MQ=255)                                 ccACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGt                                                                                                                        >  7:119539/1‑145 (MQ=255)                                  cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgtg                                                                                                                                   >  4:115611/1‑133 (MQ=255)                                       ctctCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCGTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGGGCCCGCGAAc                                                                                                              <  8:106295/149‑1 (MQ=255)                                          tcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCACGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGa                                                                                                                                <  8:43001/128‑1 (MQ=17)                                          tcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCACGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGa                                                                                                                                >  7:43001/1‑128 (MQ=17)                                                ggCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctcacc                                                                                                     <  8:6329/149‑1 (MQ=255)                                                  cacaATAGGTCGGTAGTTCAGCCTTCCCCAACCCGGAGGTCTCG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCATGAAGCTTGCGAGACACCTGTCGGAGCGGGCGTGCCCGCGAACCGcctcaact                                                                                                    >  5:380155/1‑148 (MQ=255)                                                          ggcGGTAGTTCTGCGTTCGCGTGGGCGTTGGGCTGG‑TTGTGTGGCCATTCGCGGGGGCGCTGGTTCCCACAGTTAAAGCGCAAGCGTTTAAGCTTGCGATACACGTTTGGGAGGGGGGGGGCCCGCGAACGGCCTCACCTCAATCCaa                                                                                            <  8:46611/146‑1 (MQ=255)                                                          gTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAAt                                                                                           >  2:7236/1‑149 (MQ=255)                                                           tCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATg                                                                                          >  7:112855/1‑149 (MQ=255)                                                                aGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTc                                                                                     >  6:263771/1‑149 (MQ=255)                                                                         ttCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCg                                                                            >  1:192848/1‑149 (MQ=255)                                                                         ttCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTCTGGGAGCGGGCGTCCCAGCGAACGTCCTCACCGCAATCGAAAGCCATCCAGaacacg                                                                            >  4:274804/1‑144 (MQ=255)                                                                           cGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCgg                                                                           <  5:6817/148‑1 (MQ=255)                                                                             ccAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTgg                                                                        >  8:154997/1‑149 (MQ=255)                                                                              cAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTgg                                                                        >  5:363773/1‑148 (MQ=255)                                                                                  ccGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCatat                                                                   <  7:154997/149‑1 (MQ=255)                                                                                    ggtggtCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTc                                                                 <  1:7236/149‑1 (MQ=255)                                                                                    ggtggtCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTc                                                                 <  5:94148/149‑1 (MQ=255)                                                                                        gtCTGG‑TGGTGTGGCCATGCGCGGGCGCGGTCGCTCCCACCGGTACAGCGCAAGCCTTGACGCTTGCGAGACACCTGGGGGAGGGGGGGGGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGa                                                             <  3:144831/149‑1 (MQ=255)                                                                                           tGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAAtc                                                          >  2:254069/1‑148 (MQ=255)                                                                                             g‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATg                                                          <  2:287847/147‑1 (MQ=255)                                                                                             g‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                                                        <  7:80801/149‑1 (MQ=255)                                                                                                      tCCATGCGCGGGCGCGCTCGCTCCCCCTGGAACCGCGCAAGCCTTCACGCTTGCGAGACACCTGTGGGAGGGGGCGGGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcg                                                <  6:240788/149‑1 (MQ=255)                                                                                                        cATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                                              >  4:306449/1‑149 (MQ=255)                                                                                                        cATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                                              <  3:115733/149‑1 (MQ=255)                                                                                                               gggCACGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGAGACACCTGTGGGAGTGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAAc                                        >  6:85678/1‑148 (MQ=25)                                                                                                                  cgcgcTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGc                                    <  7:101175/149‑1 (MQ=255)                                                                                                                   gcgcTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcgc                                                               <  7:143852/121‑1 (MQ=255)                                                                                                                   gcgcTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcgc                                                               >  8:143852/1‑121 (MQ=255)                                                                                                                   gcgcTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGTTCCGACGAGAACCAGcc                                   >  2:60086/1‑149 (MQ=255)                                                                                                                          ctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgt                            <  6:363773/149‑1 (MQ=255)                                                                                                                             ccACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtc                                                                               >  2:245762/1‑95 (MQ=255)                                                                                                                                  ggTACACCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCTTCCGTGCCGTGtt                    <  6:380155/149‑1 (MQ=255)                                                                                                                                         cgcAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATAcc             >  6:114057/1‑149 (MQ=255)                                                                                                                                          gcAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCg            <  3:306449/149‑1 (MQ=255)                                                                                                                                                cTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgc      >  4:205040/1‑149 (MQ=255)                                                                                                                                                 ttGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgcg     >  1:217598/1‑149 (MQ=255)                                                                                                                                                    aaGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGcca  >  3:231453/1‑149 (MQ=255)                                                                                                                                                    aaGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGcca  <  4:231453/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCA  >  NC_002947/5743218‑5743512

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