Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,388,037	G→T	6.1%	G151G (GGC→GGA)	PP_1209 ←	cold‑shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,388,037	0	G	T	6.1%	213.1 / 3.7	66	G151G (GGC→GGA)	PP_1209	cold‑shock protein
Reads supporting (aligned to +/- strand):  ref base G (27/35);  new base T (2/2);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

GTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCA  >  NC_002947/1387891‑1388182                                                                                                                                                   |                                                                                                                                                   gTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGc‑c                                                                                                                                                 >  2:251252/1‑149 (MQ=255)     gTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAg                                                                                                                                             >  6:221594/1‑149 (MQ=255)      tCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGt                                                                                                                                            >  2:3635/1‑149 (MQ=255)       cAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGtc                                                                                                                                           >  5:30470/1‑149 (MQ=255)           ggcggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCgg                                                                                                                                       <  3:221876/149‑1 (MQ=255)              ggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGgag                                                                                                                                    <  4:23157/149‑1 (MQ=255)               gTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTTCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGgaga                                                                                                                                   <  6:10605/149‑1 (MQ=255)               gTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGgaga                                                                                                                                   <  2:34283/149‑1 (MQ=255)                tttACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGAt                                                                                                                                  <  6:3319/149‑1 (MQ=255)                tttACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGAt                                                                                                                                  <  4:70011/149‑1 (MQ=255)                tttACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGAt                                                                                                                                  <  8:161051/149‑1 (MQ=255)                       cAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAAcc                                                                                                                            >  3:229338/1‑148 (MQ=255)                         aCCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTC‑CCGAGACTCGGGCGATGacatcttc                                                                                                                         <  3:247537/149‑8 (MQ=255)                         aCCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCaaa                                                                                                                         <  5:310316/149‑1 (MQ=255)                          ccACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAg                                                                                                                        >  2:70312/1‑149 (MQ=255)                          ccACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACACCACGCGATGCCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAg                                                                                                                        >  8:187832/1‑149 (MQ=255)                           cACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGc                                                                                                                       >  5:258197/1‑149 (MQ=255)                           cACGGCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTC‑CCGAGACTCGGGCGAt                                                                                                                                  >  4:219197/1‑138 (MQ=255)                           cACGCCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTC‑CCGAGACTCGGGCGAt                                                                                                                                  <  3:219197/138‑1 (MQ=255)                            aCGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGcc                                                                                                                      <  1:268346/149‑1 (MQ=255)                               tCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                   <  3:308346/149‑1 (MQ=255)                               tCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                   <  5:25552/149‑1 (MQ=255)                               tCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                   <  5:221594/149‑1 (MQ=255)                               tCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCtt                                                                                                                    >  4:231626/1‑148 (MQ=255)                                  tCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGa                                                                                                                >  4:21184/1‑149 (MQ=255)                                   cGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTc                                                                                                                  <  7:152287/146‑1 (MQ=255)                                          cAGGCCCTTGTCGCGGTGCAGCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGtt                                                                                                          <  2:257109/147‑1 (MQ=255)                                               ccTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAt                                                                                                   >  3:154921/1‑149 (MQ=255)                                                    tCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGa                                                                                               <  6:258197/148‑1 (MQ=255)                                                      gcgGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACg                                                                                             <  7:249612/148‑1 (MQ=255)                                                                 ccGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCa                                                                                 >  5:310705/1‑149 (MQ=255)                                                                 ccGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCa                                                                                 >  5:255960/1‑149 (MQ=255)                                                                  cGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCAc                                                                                >  7:80703/1‑149 (MQ=255)                                                                   gAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACg                                                                               <  3:107887/149‑1 (MQ=255)                                                                     aaaCTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGt                                                                             <  3:21184/149‑1 (MQ=255)                                                                        cTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGc                                                                          >  8:201343/1‑149 (MQ=255)                                                                         tCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGcc                                                                         >  5:273131/1‑149 (MQ=255)                                                                            aCGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGACCCATTTGACGGTGCCGGTATCACGGTTGCCgg                                                                       >  5:251000/1‑148 (MQ=255)                                                                                 cTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCtt                                                                 <  2:293622/149‑1 (MQ=255)                                                                                    gCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgtg                                                              >  6:169887/1‑149 (MQ=255)                                                                                    gCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                               <  6:255960/148‑1 (MQ=255)                                                                                         ctaCCAGATTGCGTTGGGGCTCACAGCGGATGGCCCGAAAGTGCAAGAAGATGTCTTCGG‑CCGAGTCTCGGGAGTTGATACGAAAGCGTTTCGACGTGTTGAACCATTTGACGGTGCAGGTATCAGGGTTGCCGGCATCTTGTGCCAc                                                          <  4:229338/146‑1 (MQ=255)                                                                                         cGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCAc                                                          <  3:231626/148‑1 (MQ=255)                                                                                          gACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACgt                                                        >  5:316025/1‑149 (MQ=255)                                                                                              aGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTg                                                    <  4:49286/149‑1 (MQ=255)                                                                                               gAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGa                                                   <  6:30470/149‑1 (MQ=255)                                                                                               gAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGa                                                   <  4:101997/149‑1 (MQ=255)                                                                                                 aTGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAAGCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCATGTGGCACGTTTTGAGt                                                 >  2:245390/1‑149 (MQ=255)                                                                                                      gTGGCCCTCACCGCGTATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCcgc                                            <  7:187832/149‑1 (MQ=255)                                                                                                      gTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCAGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCcgc                                            <  1:70312/149‑1 (MQ=255)                                                                                                         gCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGt                                          <  5:132087/148‑1 (MQ=255)                                                                                                          cccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGc                                        <  5:169887/149‑1 (MQ=255)                                                                                                           ccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CAGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGcc                                       >  2:189015/1‑149 (MQ=255)                                                                                                             tCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCtt                                     >  4:12379/1‑149 (MQ=255)                                                                                                               aCCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTgc                                   <  3:194702/149‑1 (MQ=255)                                                                                                                ccGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTgcg                                  <  3:12379/149‑1 (MQ=255)                                                                                                                   cgGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGt                               <  6:316025/149‑1 (MQ=255)                                                                                                                        ggCCCGAAAGTGCACGAAGATGTCATCGC‑GCGAGTCTCGGGAGATGAAACGAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGGGTGCCTTGCGCGtgctg                           <  6:251000/148‑1 (MQ=255)                                                                                                                          cccGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCATAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCgg                        <  8:80703/149‑1 (MQ=255)                                                                                                                           ccGAAAGTGCACGAAGATGTCATCGCGCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCgg                        >  1:177767/1‑149 (MQ=255)                                                                                                                              aaaGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCaa                    <  6:128011/149‑1 (MQ=255)                                                                                                                                     aCGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCAc             >  8:58695/1‑149 (MQ=255)                                                                                                                                       gAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACgg           >  3:8477/1‑149 (MQ=255)                                                                                                                                             gTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGAt     >  5:220777/1‑149 (MQ=255)                                                                                                                                              tCATCTC‑CCGAGACTCGGGCGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATg    >  7:35554/1‑149 (MQ=255)                                                                                                                                                aTCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGca  >  3:206863/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCA  >  NC_002947/1387891‑1388182

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