Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,885,013	G→C	5.4%	L132F (TTG→TTC)	PP_1692 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,885,013	0	G	C	5.4%	303.4 / 5.1	93	L132F (TTG→TTC)	PP_1692	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (44/44);  new base C (2/3);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

ATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCGGGTGGAGATTG  >  NC_002947/1884868‑1885159                                                                                                                                                  |                                                                                                                                                    aTGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCttgttg                                                                                                                                                  >  5:375429/1‑149 (MQ=255)  tGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGc                                                                                                                                                 <  6:176636/149‑1 (MQ=255)      ctgctgCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCGTGTTGCGa                                                                                                                                               >  8:312278/1‑147 (MQ=255)       tgctgcGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTGAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAACAAGCGCCGCGCTTCCTGCGAccc                                                                                                                                            <  5:36524/149‑3 (MQ=255)         ctgcGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTg                                                                                                                                          >  5:23181/1‑149 (MQ=255)         ctgcGTCAGGCGTCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCGGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGTCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGt                                                                                                                                           <  8:96854/148‑1 (MQ=255)            cGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGt                                                                                                                                        <  5:373423/148‑1 (MQ=255)             gTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGttt                                                                                                                                      >  5:237110/1‑149 (MQ=255)             gTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGttt                                                                                                                                      >  5:258122/1‑149 (MQ=255)             gTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGGCGCAACAAGCGCCGCGCTTCCTGCGAcc                                                                                                                                             <  7:87956/142‑2 (MQ=255)             gTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGCCGCAACAAGCGCCGCGCTTCCTGCGAcc                                                                                                                                             >  8:87956/1‑141 (MQ=255)                aGGCTGACCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTtct                                                                                                                                   <  2:75348/149‑1 (MQ=255)                  gCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTtctgt                                                                                                                                 >  8:204785/1‑149 (MQ=255)                  gCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTtctgt                                                                                                                                 <  1:207870/149‑1 (MQ=255)                   cTGGCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGGCAGCACAACGTGGCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTtctgtc                                                                                                                                <  5:227420/149‑1 (MQ=255)                   cTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTtctgtc                                                                                                                                <  6:372140/149‑1 (MQ=255)                   agcctgacgcAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTtctgtc                                                                                                                                <  5:296551/142‑1 (MQ=255)                     gCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTtctgtct                                                                                                                               >  6:247814/1‑148 (MQ=255)                        cgcgCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGc                                                                                                                           >  2:311778/1‑149 (MQ=255)                                ggATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGc                                                                                                                   <  5:222732/149‑1 (MQ=255)                                  aTCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCt                                                                                                                 <  6:375429/149‑1 (MQ=255)                                   tCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTg                                                                                                                >  1:387640/1‑149 (MQ=255)                                     ggcgggcgCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGc                                                                                                               <  8:89660/148‑1 (MQ=255)                                      gcgggcgCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCAc                                                                                                             >  4:182569/1‑149 (MQ=255)                                      gcgggcgCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGCCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCAc                                                                                                             >  2:107613/1‑149 (MQ=255)                                       cgggcgCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACg                                                                                                            >  8:328013/1‑149 (MQ=255)                                        gggcgCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACgg                                                                                                           <  1:128848/149‑1 (MQ=255)                                         ggcgCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACggg                                                                                                          >  5:319311/1‑149 (MQ=255)                                            gcAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGtttt                                                                                                                                     <  6:226871/119‑1 (MQ=255)                                            gcAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGtttt                                                                                                                                     >  5:226871/1‑119 (MQ=255)                                                   gAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTcgc                                                                                                 <  6:98130/148‑1 (MQ=255)                                                    aGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGa                                                                                               <  5:247814/149‑1 (MQ=255)                                                    aGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGa                                                                                               <  6:258122/149‑1 (MQ=255)                                                    aGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGa                                                                                               <  5:290759/149‑1 (MQ=255)                                                    aGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGa                                                                                               <  4:66724/149‑1 (MQ=255)                                                      tACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAg                                                                                              >  4:81823/1‑148 (MQ=255)                                                         gCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAAt                                                                                          >  8:72661/1‑149 (MQ=255)                                                         gCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTAATGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCaa                                                                                           >  5:38124/1‑148 (MQ=255)                                                            cGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCAc                                                                                       <  6:23181/149‑1 (MQ=255)                                                              gCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCg                                                                                     <  2:71707/149‑1 (MQ=255)                                                                  tgtCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTc                                                                                                   <  7:38321/131‑1 (MQ=255)                                                                  tgtCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTc                                                                                                   >  8:38321/1‑131 (MQ=255)                                                                  tgtCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGAtgt                                                                                 <  6:319311/149‑1 (MQ=255)                                                                    tCAGCACAACGTGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGTCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGAtgtg                                                                                <  3:182569/148‑1 (MQ=255)                                                                      aGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATgtgtgt                                                                             >  1:270737/1‑149 (MQ=255)                                                                       gCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATgtgtgt                                                                             <  8:88585/148‑1 (MQ=255)                                                                       gCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATgtgtgt                                                                             >  7:88585/1‑148 (MQ=255)                                                                       gCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTc                                                                            <  3:78459/149‑1 (MQ=255)                                                                       gCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTc                                                                            >  1:302295/1‑149 (MQ=255)                                                                          caACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCAc                                                                         >  4:317671/1‑149 (MQ=255)                                                                           aaCGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCAcc                                                                        >  8:352314/1‑149 (MQ=255)                                                                           aaCGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCAcc                                                                        >  1:26535/1‑149 (MQ=255)                                                                                  cGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTg‑g                                                                >  3:4428/1‑149 (MQ=255)                                                                                   ggAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑Gt                                                               <  8:37933/149‑1 (MQ=255)                                                                                    gAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCgccgccaccgta                                                               >  6:400573/1‑143 (MQ=255)                                                                                           cgcgcTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCgg                                                       <  7:352314/149‑1 (MQ=255)                                                                                            gcgcTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGa                                                                                               <  3:53452/109‑1 (MQ=255)                                                                                            gcgcTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGa                                                                                               >  4:53452/1‑109 (MQ=255)                                                                                             cgcTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGt                                                     >  4:133053/1‑149 (MQ=255)                                                                                              gcTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGtt                                                    <  2:249915/149‑1 (MQ=255)                                                                                                  ttGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTAcc                                                >  3:130521/1‑149 (MQ=255)                                                                                                    ggCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTATTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCcg                                              >  1:111531/1‑149 (MQ=255)                                                                                                     gCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCcgc                                             <  1:311778/149‑1 (MQ=255)                                                                                                      cTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCTCGAGCAATCACCGATGTGTGTCAACCGCCGTG‑GTAGGCGCGGGTTTACCcgc                                             >  6:327440/1‑148 (MQ=255)                                                                                                       ttGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCcgcg                                            >  8:253062/1‑148 (MQ=255)                                                                                                          aGCTTGCGGCAGGGCAGGGGGCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTAAACCcgcgcct                                         <  6:38124/148‑4 (MQ=255)                                                                                                            cTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACg                                      >  5:200456/1‑149 (MQ=255)                                                                                                            cTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACg                                      >  3:226010/1‑149 (MQ=255)                                                                                                                 ggcagggcagGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTggc                                 >  6:333329/1‑149 (MQ=255)                                                                                                                    agggcagggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTggcggc                              <  6:237110/149‑1 (MQ=255)                                                                                                                    agggcagggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGC‑TTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTggcggcg                             >  1:29229/1‑149 (MQ=255)                                                                                                                      ggccgggGGCGCTGCACGCGCGGCGCTTGTGCCCACTGTGGTTTTCTGTCTGAGCTTTTTTGCCCCCCCGGCCCCTTTCCCCACCAACCACCGATGTGTGCCCCCCGCCGTG‑GTCGGCGCGGGTTAACCCGCTAATACGGGggcgcc                              <  2:302295/147‑3 (MQ=255)                                                                                                                      ggcagggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTggcggcgg                            >  2:176278/1‑149 (MQ=255)                                                                                                                      ggcagggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTggcggcgg                            >  6:35776/1‑149 (MQ=255)                                                                                                                      ggcagggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTggcggcgg                            >  3:358183/1‑149 (MQ=255)                                                                                                                       gcagggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTggcggcggc                           <  4:130521/149‑1 (MQ=255)                                                                                                                         cactgTCGCAACAAGCGCCGCGCTTCCTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGtt                                                    <  8:54850/118‑1 (MQ=38)                                                                                                                         cactgTCGCAACAAGCGCCGCGCTTCCTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGtt                                                    >  7:54850/5‑122 (MQ=38)                                                                                                                         agggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCaa                         >  1:350955/1‑149 (MQ=255)                                                                                                                          ggggTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAAc                        >  5:334204/1‑149 (MQ=255)                                                                                                                             gTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGt                     <  2:146713/149‑1 (MQ=255)                                                                                                                                gCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGaa                  >  4:248766/1‑149 (MQ=255)                                                                                                                                gCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGaa                  >  8:20278/1‑149 (MQ=255)                                                                                                                                 cAGTAAGCGCGGCGCTTGTTGCGACAGTTGTTGTCTGTCGGAGCTTTTTTGCCTGCCCGGGCCCGTGCGCGCGCAATCACCGATGTGTGTCCACCGCCGTT‑GTAGGCGCGGGTTTACACGCGAATACGGTGGCGGCGGCAACGGTGAAt                 <  2:316312/149‑1 (MQ=25)                                                                                                                                 cAGGACGCGCGGCGCTTGTTGCGACAGTGGTGTTCTGTCTGAGCTGGTTTGCCTGCACGGGCCCTTTCGCGCGCAAGCCCCGATGTGGGGCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAAt                 <  2:321088/149‑1 (MQ=25)                                                                                                                                  aGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACggg                                                                                                          <  2:45120/60‑1 (MQ=255)                                                                                                                                  aGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACggg                                                                                                          >  1:45120/1‑60 (MQ=255)                                                                                                                                   ggAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATgg               <  7:72661/149‑1 (MQ=255)                                                                                                                                      aGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCgg            >  8:11043/1‑149 (MQ=255)                                                                                                                                      aGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCgg            <  3:81823/149‑1 (MQ=255)                                                                                                                                      aGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCgg            <  7:253062/149‑1 (MQ=255)                                                                                                                                          cggcggTTGTTGCGTCAGTGGTTTTCTGGCTGTGGTTTTTTGCCTTCACGGGGCCTTTCGCGAGCAATCACCGCTGTGTGGCCACCGGCGTGCGTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCGGGTg         <  1:107613/149‑1 (MQ=25)                                                                                                                                          cggcgCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCGGGTgg        <  3:248766/149‑1 (MQ=255)                                                                                                                                              gcTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCGGGTGGAGAt    <  5:35776/149‑1 (MQ=255)                                                                                                                                                  gttgCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCGGGTGGAGATTg  <  2:350955/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    ATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCGGGTGGAGATTG  >  NC_002947/1884868‑1885159

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