Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,885,006	G→C	5.3%	R130P (CGG→CCG)	PP_1692 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,885,006	0	G	C	5.3%	313.2 / 5.7	95	R130P (CGG→CCG)	PP_1692	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (48/42);  new base C (2/3);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

ATAGGCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCGGGTGG  >  NC_002947/1884861‑1885153                                                                                                                                                  |                                                                                                                                                     ataGGCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCgcggcg                                                                                                                                                   >  6:222732/1‑149 (MQ=255)  taGGCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGgcgc                                                                                                                                                  <  5:68004/149‑1 (MQ=255)   aGGCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCt                                                                                                                                                 >  7:37933/1‑149 (MQ=255)    ggCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCtt                                                                                                                                                >  5:176636/1‑149 (MQ=255)    ggCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCtt                                                                                                                                                >  5:98130/1‑149 (MQ=255)        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:258122/1‑149 (MQ=255)                    gTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGttt                                                                                                                                >  5:237110/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                                                                                                                           <  1:207870/149‑1 (MQ=255)                         gCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTtctgt                                                                                                                           >  8:204785/1‑149 (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                                                                      >  1:302295/1‑149 (MQ=255)                                                                              gCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTc                                                                      <  3:78459/149‑1 (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                                              >  7:54850/5‑122 (MQ=38)                                                                                                                                cactgTCGCAACAAGCGCCGCGCTTCCTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGtt                                              <  8:54850/118‑1 (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            >  8:20278/1‑149 (MQ=255)                                                                                                                                       gCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGaa            >  4:248766/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      <  7:253062/149‑1 (MQ=255)                                                                                                                                             aGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCgg      <  3:81823/149‑1 (MQ=255)                                                                                                                                             aGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCgg      >  8:11043/1‑149 (MQ=255)                                                                                                                                                 cggcggTTGTTGCGTCAGTGGTTTTCTGGCTGTGGTTTTTTGCCTTCACGGGGCCTTTCGCGAGCAATCACCGCTGTGTGGCCACCGGCGTGCGTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCGGGTg   <  1:107613/149‑1 (MQ=25)                                                                                                                                                 cggcgCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCGGGTgg  <  3:248766/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     ATAGGCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTG‑GTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGGCGGGTGG  >  NC_002947/1884861‑1885153

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