Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,743,335	T→C	7.9%	intergenic (‑168/+79)	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,335	0	T	C	7.9%	179.2 / 4.9	63	intergenic (‑168/+79)	fbp/glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (29/29);  new base C (3/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 = 3.57e-01

GGGGGATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCC  >  NC_002947/5743188‑5743479                                                                                                                                                     |                                                                                                                                                 gggggATGAAAAACCCACGGAGTTTACTGGGAGAGTGACGGGCTAGCGAGTTATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTTGTCTGG‑TTGTGTGTCCATTCGCGGGCGcgct                                                                                                                                                  <  6:20589/148‑1 (MQ=255)     gATGAAAAACCCAGGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCGTTCGCGGGCACgcccgctc                                                                                                                                             >  6:98722/1‑149 (MQ=255)      aTGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTcc                                                                                                                                            >  1:287847/1‑149 (MQ=255)                 aCGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGt                                                                                                                                     <  3:353798/145‑1 (MQ=255)                 aCGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGt                                                                                                                                     >  4:353798/1‑145 (MQ=255)                 aCGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAg                                                                                                                                 >  7:106295/1‑149 (MQ=255)                 aCGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAg                                                                                                                                 >  3:215613/1‑149 (MQ=255)                 aCGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAg                                                                                                                                 >  7:283556/1‑149 (MQ=255)                  cGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTTGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAgc                                                                                                                                >  4:79909/1‑149 (MQ=255)                           cAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCtt                                                                                                                       <  1:338485/149‑1 (MQ=255)                            aGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCtt                                                                                                                       >  3:235914/1‑148 (MQ=255)                            aGGGAGAGTGACGGGCTAGCGAGGGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTg                                                                                                                      <  1:177249/149‑1 (MQ=255)                               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‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGcc                                                                                                                         >  4:12966/1‑132 (MQ=255)                                          gCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGcc                                                                                                                         <  3:12966/132‑1 (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                                                                                       >  7:119539/1‑145 (MQ=255)                                                               ccACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGt                                                                                       <  8:119539/145‑1 (MQ=255)                                                                cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgtg                                                                                                  >  4:115611/1‑133 (MQ=255)                                                                     ctctCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCGTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGGGCCCGCGAAc                                                                             <  8:106295/149‑1 (MQ=255)                                                                        tcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCACGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGa                                                                                               >  7:43001/1‑128 (MQ=17)                                                                        tcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCACGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGa                                                                                               <  8:43001/128‑1 (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                              >  8:143852/1‑121 (MQ=255)                                                                                                                                                 gcgcTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcgc                              <  7:143852/121‑1 (MQ=255)                                                                                                                                                 gcgcTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGTTCCGACGAGAACCAGcc  >  2:60086/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GGGGGATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGG‑TTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCC  >  NC_002947/5743188‑5743479

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