Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,884,998	G→A	5.3%	Q127Q (CAG→CAA)	PP_1692 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,884,998	0	G	A	5.3%	292.3 / 3.1	94	Q127Q (CAG→CAA)	PP_1692	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (49/40);  new base A (2/3);  total (51/43)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

ACCGGCGTATAGGCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGG  >  NC_002947/1884853‑1885146                                                                                                                                                  |                                                                                                                                                     aCCGGCGTATAGGCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGaa                                                                                                                                                   <  1:392401/149‑1 (MQ=255)   cGGCGTATAGGCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAgc                                                                                                                                                 >  2:128848/1‑149 (MQ=255)         ataGGCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCgcggcg                                                                                                                                           >  6:222732/1‑149 (MQ=255)          taGGCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGgcgc                                                                                                                                          <  5:68004/149‑1 (MQ=255)           aGGCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCt                                                                                                                                         >  7:37933/1‑149 (MQ=255)            ggCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCtt                                                                                                                                        >  5:98130/1‑149 (MQ=255)            ggCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCtt                                                                                                                                        >  5:176636/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                                                                                 <  6:258122/149‑1 (MQ=255)                                                                   aGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGa                                                                                 <  4:66724/149‑1 (MQ=255)                                                                   aGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGa                                                                                 <  5:290759/149‑1 (MQ=255)                                                                   aGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGa                                                                                 <  5:247814/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                                                                                     >  8:38321/1‑131 (MQ=255)                                                                                 tgtCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTc                                                                                     <  7:38321/131‑1 (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)                                                                                                 cGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTgg                                                   >  3:4428/1‑149 (MQ=255)                                                                                                  ggAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGt                                                  <  8:37933/149‑1 (MQ=255)                                                                                                   gAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCgccgccaccgta                                                 >  6:400573/1‑143 (MQ=255)                                                                                                          cgcgcTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCgg                                          <  7:352314/149‑1 (MQ=255)                                                                                                           gcgcTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGa                                                                                 <  3:53452/109‑1 (MQ=255)                                                                                                           gcgcTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGa                                                                                 >  4:53452/1‑109 (MQ=255)                                                                                                            cgcTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGt                                        >  4:133053/1‑149 (MQ=255)                                                                                                             gcTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGtt                                       <  2:249915/149‑1 (MQ=255)                                                                                                                 ttGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTAcc                                   >  3:130521/1‑149 (MQ=255)                                                                                                                   ggCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTATTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCcg                                 >  1:111531/1‑149 (MQ=255)                                                                                                                    gCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCcgc                                <  1:311778/149‑1 (MQ=255)                                                                                                                     cTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCTCGAGCAATCACCGATGTGTGTCAACCGCCGTGGTAGGCGCGGGTTTACCcgc                                >  6:327440/1‑148 (MQ=255)                                                                                                                      ttGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCcgcg                               >  8:253062/1‑148 (MQ=255)                                                                                                                         aGCTTGCGGCAGGGCAGGGGGCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTAAACCcgcgcct                            <  6:38124/148‑4 (MQ=255)                                                                                                                           cTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACg                         >  3:226010/1‑149 (MQ=255)                                                                                                                           cTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACg                         >  5:200456/1‑149 (MQ=255)                                                                                                                                ggcagggcagGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTggc                    >  6:333329/1‑149 (MQ=255)                                                                                                                                   agggcagggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTggcggc                 <  6:237110/149‑1 (MQ=255)                                                                                                                                   agggcagggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGC‑TTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTggcggcg                >  1:29229/1‑149 (MQ=255)                                                                                                                                     ggccgggGGCGCTGCACGCGCGGCGCTTGTGCCCACTGTGGTTTTCTGTCTGAGCTTTTTTGCCCCCCCGGCCCCTTTCCCCACCAACCACCGATGTGTGCCCCCCGCCGTGGTCGGCGCGGGTTAACCCGCTAATACGGGggcgcc                 <  2:302295/147‑3 (MQ=255)                                                                                                                                     ggcagggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTggcggcgg               >  6:35776/1‑149 (MQ=255)                                                                                                                                     ggcagggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTggcggcgg               >  2:176278/1‑149 (MQ=255)                                                                                                                                     ggcagggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTggcggcgg               >  3:358183/1‑149 (MQ=255)                                                                                                                                      gcagggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTggcggcggc              <  4:130521/149‑1 (MQ=255)                                                                                                                                        cactgTCGCAACAAGCGCCGCGCTTCCTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGtt                                       >  7:54850/5‑122 (MQ=38)                                                                                                                                        cactgTCGCAACAAGCGCCGCGCTTCCTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGtt                                       <  8:54850/118‑1 (MQ=38)                                                                                                                                        agggGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCaa            >  1:350955/1‑149 (MQ=255)                                                                                                                                         ggggTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAAc           >  5:334204/1‑149 (MQ=255)                                                                                                                                            gTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGt        <  2:146713/149‑1 (MQ=255)                                                                                                                                               gCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGaa     >  8:20278/1‑149 (MQ=255)                                                                                                                                               gCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGaa     >  4:248766/1‑149 (MQ=255)                                                                                                                                                cAGTAAGCGCGGCGCTTGTTGCGACAGTTGTTGTCTGTCGGAGCTTTTTTGCCTGCCCGGGCCCGTGCGCGCGCAATCACCGATGTGTGTCCACCGCCGTTGTAGGCGCGGGTTTACACGCGAATACGGTGGCGGCGGCAACGGTGAAt    <  2:316312/149‑1 (MQ=25)                                                                                                                                                cAGGACGCGCGGCGCTTGTTGCGACAGTGGTGTTCTGTCTGAGCTGGTTTGCCTGCACGGGCCCTTTCGCGCGCAAGCCCCGATGTGGGGCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAAt    <  2:321088/149‑1 (MQ=25)                                                                                                                                                 aGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACggg                                                                                            >  1:45120/1‑60 (MQ=255)                                                                                                                                                 aGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACggg                                                                                            <  2:45120/60‑1 (MQ=255)                                                                                                                                                  ggAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATgg  <  7:72661/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     ACCGGCGTATAGGCCATGAGCTGCTGCGTCAGGCTGCCCGCGCAGGGGATCGGCGGGCGCAGCTGGAGTACGCCCGGCTGTGTCAGCACAACGAGCCGGAGCAGGCGCGCTATTGGCTTGAGCTTGCGGCAGGGCAGGGGTCGCAGGAAGCGCGGCGCTTGTTGCGACAGTGGTTTTCTGTCTGAGCTTTTTTGCCTGCACGGGCCCTTTCGCGAGCAATCACCGATGTGTGTCCACCGCCGTGGTAGGCGCGGGTTTACCCGCGAATACGGTGGCGGCGGCAACGGTGAATGG  >  NC_002947/1884853‑1885146

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