Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,388,044	T→A	6.1%	D149V (GAC→GTC)	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,044	0	T	A	6.1%	198.5 / 2.9	66	D149V (GAC→GTC)	PP_1209	cold‑shock protein
Reads supporting (aligned to +/- strand):  ref base T (26/36);  new base A (2/2);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

TCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCAC  >  NC_002947/1387896‑1388189                                                                                                                                                      |                                                                                                                                                  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                                                                                                                                         <  4:70011/149‑1 (MQ=255)           tttACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGAt                                                                                                                                         <  6:3319/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                                                                                                                          <  5:221594/149‑1 (MQ=255)                          tCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                          <  5:25552/149‑1 (MQ=255)                          tCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                          <  3:308346/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)                                                                                                                                               c‑ccGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacggc     >  3:87868/1‑149 (MQ=255)                                                                                                                                                 ccGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGGCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacggcc    <  6:273131/149‑1 (MQ=255)                                                                                                                                                 ccGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacggc     >  3:317774/1‑148 (MQ=255)                                                                                                                                                   gAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacggccac  >  6:113499/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                  TCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGC‑CCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCAC  >  NC_002947/1387896‑1388189

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