Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	962,364	A→G	5.7%	G17G (GGT→GGC)	selD ←	selenide, water dikinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	962,364	0	A	G	5.7%	201.6 / 2.0	70	G17G (GGT→GGC)	selD	selenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base A (30/36);  new base G (2/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

GTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCTTTAGGCCCTTGCCTCTTATAATTCCCGGTGTCGGCAAATGTGTC  >  NC_002947/962221‑962507                                                                                                                                                |                                                                                                                                                   gtCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACcaca                                                                                                                                               <  3:155136/149‑1 (MQ=255)  tCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAg                                                                                                                                              >  6:84657/1‑149 (MQ=255)    gagaCCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGcc                                                                                                                                            >  4:142022/1‑149 (MQ=255)          aCGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAGCCACAGCCGgc‑‑gc                                                                                                                                     >  7:8926/1‑148 (MQ=255)          aCGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAGCCACAGCCGgc‑‑gc                                                                                                                                     >  7:14206/1‑148 (MQ=255)            gccgcGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCAt                                                                                                                                  <  4:237615/149‑1 (MQ=255)            gccgcGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCAt                                                                                                                                  >  2:338120/1‑149 (MQ=255)                 gcTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTg                                                                                                                             <  4:199992/149‑1 (MQ=255)                     gtcgtcCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACt                                                                                                                         >  5:267748/1‑149 (MQ=255)                        gtcCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTggg                                                                                                                      >  4:228233/1‑149 (MQ=255)                           cAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGt                                                                                                                     >  2:57577/1‑147 (MQ=255)                              cgcgTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑c                                                                                                               <  6:267748/149‑1 (MQ=255)                              cgcgTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑c                                                                                                               <  1:308796/149‑1 (MQ=255)                              cgcgTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑c                                                                                                               <  1:338120/149‑1 (MQ=255)                                 gTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGaa                                                                                                            >  5:274016/1‑149 (MQ=255)                                   acacTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAAc                                                                                                                                                   <  4:234157/111‑1 (MQ=255)                                   acacTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAAc                                                                                                                                                   >  3:234157/1‑111 (MQ=255)                                       tgctgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGgct                                                                                                      >  2:22603/1‑149 (MQ=255)                                       tgctgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGgct                                                                                                      >  8:62350/1‑149 (MQ=255)                                       tgctgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCAgct                                                                                                      >  3:62040/1‑149 (MQ=255)                                       tgctgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCCTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGgct                                                                                                      >  8:55448/1‑149 (MQ=255)                                         ctgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAGCCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGgctcg                                                                                                    <  8:8926/149‑1 (MQ=255)                                         ctgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAGCCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGgctcg                                                                                                    <  8:14206/149‑1 (MQ=255)                                          tgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACt                                                                                                                         >  1:95166/1‑128 (MQ=255)                                          tgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACt                                                                                                                         <  2:95166/128‑1 (MQ=255)                                          tgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGgctcgc                                                                                                   >  3:312899/1‑149 (MQ=255)                                          tgcATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGgctcgc                                                                                                   >  5:209108/1‑149 (MQ=255)                                          tgcATCGTCGAGCGATGCGTTGCCGACCAACAGTTTCTGGTCCATGTCCTGGGTGCCGCTTTCGTCAAGGATCACGTCCAGCAGCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCCGTACTGGGTCAGG‑CGAATCGgctcgc                                                                                                   <  7:312664/149‑1 (MQ=255)                                            cATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGgctcgc                                                                                                   <  2:224425/147‑1 (MQ=255)                                            cATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGgctcgc                                                                                                   <  7:51200/147‑1 (MQ=255)                                              tcgtcgCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCAt                                                                                               <  4:312899/149‑1 (MQ=255)                                              tcgtcgCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCAt                                                                                               <  6:209108/149‑1 (MQ=255)                                              tcgtcgCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCAGCTCGCTCAt                                                                                               <  4:62040/149‑1 (MQ=255)                                                gtcgCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCCTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATa                                                                                              <  7:55448/148‑1 (MQ=255)                                                    cgcgATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGgctc                                                                                                     >  1:85653/1‑137 (MQ=255)                                                    cgcgATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGgctc                                                                                                     <  2:85653/137‑1 (MQ=255)                                                             ttGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAgg                                                                                >  7:323224/1‑149 (MQ=255)                                                                   aCCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGGCCGAATCGGCTCGCTCATACGCACCTCTGCAGGTCaaa                                                                           <  1:57577/149‑1 (MQ=255)                                                                    cccACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAAtt                                                                         <  1:22603/149‑1 (MQ=255)                                                                     ccACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAAttg                                                                        >  1:56836/1‑149 (MQ=255)                                                                     ccACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAAttg                                                                        >  6:299786/1‑149 (MQ=255)                                                                       acaGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAAttgtt                                                                      <  7:62350/149‑1 (MQ=255)                                                                       acaGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAAttgtt                                                                      <  5:272458/149‑1 (MQ=255)                                                                       acaGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAAttgtt                                                                      <  5:299786/149‑1 (MQ=255)                                                                             tCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGcgcg                                                                 <  3:228233/148‑1 (MQ=255)                                                                             tCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGcgcg                                                                 <  3:142022/148‑1 (MQ=255)                                                                                 gTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTg                                                             >  5:206320/1‑148 (MQ=255)                                                                                  tCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTg                                                             >  4:50831/1‑147 (MQ=255)                                                                                    cAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAg                                                          <  6:274016/148‑1 (MQ=255)                                                                                    cAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGc                                                         <  2:56836/149‑1 (MQ=255)                                                                                     aGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGc                                                         >  5:233559/1‑148 (MQ=255)                                                                                       ggCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCaaa                                                      <  3:50831/149‑1 (MQ=255)                                                                                         ccTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCTGGGCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCa                                                        >  1:282586/1‑147 (MQ=255)                                                                                          cTGGGTGCAGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTCTACTGGGTCAGG‑CGAATCGGCTCGCTCAAACGCACCTATCCATGTCAAATAGTTCCGCGAATGTAGCAAAGCt                                                   >  2:85323/1‑149 (MQ=255)                                                                                                   gCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGTGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATTCGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCTTTAg                                          <  1:85323/149‑1 (MQ=255)                                                                                                         ggCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCAATTGTAGCAAAGCTGGCCTTTAGGCCCtt                                    >  5:351335/1‑149 (MQ=255)                                                                                                          gCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAAttgt                                                                       >  4:27337/1‑113 (MQ=255)                                                                                                          gCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAAttgt                                                                       <  3:27337/113‑1 (MQ=255)                                                                                                           cAAGGATCACCTCCAGGACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCTTTAGGCCCTTGc                                  <  8:170901/149‑1 (MQ=255)                                                                                                                 tCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCACATTGTTGCGCGATTGTAGCAAATCTGGCCTTTAGGCCCTGGCctcgt                             >  4:230245/1‑146 (MQ=255)                                                                                                                 tCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCTTTAGGccc                                      <  7:285115/139‑1 (MQ=255)                                                                                                                 tCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCTTTAGGccc                                      >  8:285115/1‑139 (MQ=255)                                                                                                                 tCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCTTTAGGCCCTTGCCTCTta                            <  2:282586/149‑1 (MQ=255)                                                                                                                   aCGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCTTTAGGCCCTTGCCTCTtat                           <  6:206320/148‑1 (MQ=255)                                                                                                                          ggaCCTTGGGGGAGATCGTGCAACCACCGCCGGC‑‑GCCATGGCGGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCTTTAGGCCCTTGCCTCTTATAATTCCCg                   <  3:230245/147‑1 (MQ=255)                                                                                                                                 gggggAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCCCCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCTTTAGGCCCTTGCCTCTTATAATTCCCGGTGTCgg            >  2:180311/1‑149 (MQ=255)                                                                                                                                      agaTCTTGCAACTACAGCCGGC‑‑GCCATGGCTGTAGTGGGTCAGG‑CGAATAGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCTTTAGGCCCTTGCCTCTTATAATTCCCGGTGTCGGCAAAt       <  8:323224/149‑1 (MQ=255)                                                                                                                                          cTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCttt                                            <  3:245907/108‑1 (MQ=255)                                                                                                                                          cTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCttt                                            >  4:245907/1‑108 (MQ=255)                                                                                                                                             gCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCTTTAGGCCCTTGCCTCTTATAATTCCCGGTGTCGGCAAATGTGTc  <  6:233559/147‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   GTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGC‑‑GCCATGGCTGTACTGGGTCAGG‑CGAATCGGCTCGCTCATACGCACCTCTGCAGGTCAAATTGTTGCGCGATTGTAGCAAAGCTGGCCTTTAGGCCCTTGCCTCTTATAATTCCCGGTGTCGGCAAATGTGTC  >  NC_002947/962221‑962507

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑