Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	160,310	C→T	5.9%	E27K (GAG→AAG)	PP_0149 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	160,310	0	C	T	5.9%	250.9 / 4.2	85	E27K (GAG→AAG)	PP_0149	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/41);  new base T (2/3);  total (40/45)
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.33e-01

GCAGGCATGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGG  >  NC_002947/160168‑160458                                                                                                                                               |                                                                                                                                                     gcaggcaTGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTg                                                                                                                                                <  6:263615/149‑1 (MQ=255)  caggcaTGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTg                                                                                                                                                >  2:317346/1‑148 (MQ=255)   aggcaTGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCa                                                                                                                                              <  7:299413/149‑1 (MQ=255)    ggcaTGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTTGAGCTGCAg                                                                                                                                             <  6:189264/149‑1 (MQ=255)      caTGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGt                                                                                                                                           >  2:267007/1‑149 (MQ=255)          ttgtCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGc                                                                                                                                       <  1:193145/149‑1 (MQ=255)            gtCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGcc                                                                                                                                      >  1:201517/1‑148 (MQ=255)            gtCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACACCTGCAGCTTGAGCTGCAggcgtccg                                                                                                                                      <  7:191219/148‑7 (MQ=255)               gCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAg                                                                                                                                   <  3:253888/148‑1 (MQ=255)                cTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAggc                                                                                                                                 >  3:65970/1‑149 (MQ=255)                 tCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAggcg                                                                                                                                >  3:147291/1‑149 (MQ=255)                    aggtGACCCGCACGCGGCGGTCCAGGTGTTGCTTGAAGCGCGGTGGCTGCTTGCCCATGGCCTCGTCGAGCTGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAggcgggc                                                                                                                             <  4:173868/146‑1 (MQ=255)                    aGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAggcgggc                                                                                                                             <  7:53757/149‑1 (MQ=255)                      cTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCa                                                                                                                                                     >  3:222869/1‑123 (MQ=255)                      cTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCa                                                                                                                                                     <  4:222869/123‑1 (MQ=255)                       tGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTg                                                                                                                          >  1:242653/1‑149 (MQ=255)                         aCCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGcc                                                                                                                        >  1:81383/1‑149 (MQ=255)                         aCCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGcc                                                                                                                        <  2:201517/149‑1 (MQ=255)                            cGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGag                                                                                                                     <  4:147291/149‑1 (MQ=255)                                   ggcggTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCAGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGt                                                                                                               >  4:144387/1‑148 (MQ=255)                                    gcggTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCa                                                                                                             <  2:18362/149‑1 (MQ=255)                                       gTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAggg                                                                                                          >  7:156557/1‑149 (MQ=255)                                        tCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTTGAGCTGCAggc                                                                                                                                           <  5:20855/115‑2 (MQ=255)                                        tCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTTGAGCTGCAggc                                                                                                                                           >  6:20855/1‑114 (MQ=255)                                         ccAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGgcg                                                                                                        <  7:161403/149‑1 (MQ=255)                                              tGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGcc                                                                                                   <  1:317346/149‑1 (MQ=255)                                                ttGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCaa                                                                                                 <  7:215240/149‑1 (MQ=255)                                                ttGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCaa                                                                                                 <  1:267007/149‑1 (MQ=255)                                                   cTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAag                                                                                                >  1:270323/1‑147 (MQ=255)                                                     tGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTg                                                                                                                                                >  2:138006/1‑97 (MQ=255)                                                     tGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTg                                                                                                                                                <  1:138006/97‑1 (MQ=255)                                                      gAAGCGCGGTGGCAGCTTGCCCATGGCCACGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAagcagc                                                                                            >  1:180664/1‑148 (MQ=255)                                                              gTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGgcacgc                                                                                   <  5:49728/149‑1 (MQ=255)                                                                 gCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGgcacgcac                                                                                 <  2:180664/148‑1 (MQ=255)                                                                 gCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGcacgcacg                                                                                >  6:114431/1‑149 (MQ=255)                                                                  cAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCacgcacgc                                                                               <  6:294679/149‑1 (MQ=255)                                                                  cAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTAGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCATCCAGgcacgaccgc                                                                               >  8:173654/1‑144 (MQ=255)                                                                    gCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCTGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACgcgc                                                                             <  3:144387/149‑1 (MQ=255)                                                                        gCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCtt                                                                           <  4:65970/147‑1 (MQ=255)                                                                          ccATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTcaca                                                                       <  7:173654/149‑1 (MQ=255)                                                                          ccATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTca                                                                         >  1:186684/1‑147 (MQ=255)                                                                          ccATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCACGCACGCACGCCCTTca                                                                         >  3:11853/1‑147 (MQ=255)                                                                           aaTGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGATGCAGGTCGGCCAAGGCGGGCGTGCCGCGAAGCGTCAGGGCGCAGCCAAGAAGCCAGGCACGCACGCGATTCACAg                                                                      >  2:254620/2‑149 (MQ=255)                                                                             tGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTg                                                                    <  2:270323/149‑1 (MQ=255)                                                                              ggCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGTCAAGCAGcc                                                                                           <  5:29930/125‑1 (MQ=255)                                                                              ggCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGTCAAGCAGcc                                                                                           >  6:29930/1‑125 (MQ=255)                                                                               gCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGGCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGc                                                                   <  3:171502/148‑1 (MQ=255)                                                                                    gtcgAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTAt                                                             >  8:25765/1‑149 (MQ=255)                                                                                     tcgAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAgg                                                                                        >  6:2333/1‑121 (MQ=255)                                                                                     tcgAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAgg                                                                                        <  5:2333/121‑1 (MQ=255)                                                                                     tcgAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATg                                                            <  8:156557/149‑1 (MQ=255)                                                                                       gAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGc                                                          >  2:121138/1‑149 (MQ=255)                                                                                        aGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCt                                                         >  1:247256/1‑149 (MQ=255)                                                                                           aGGGCTTGGGAGGCTTTAGTTTGTTGTGTGTCGAGGCCGGACGCCTGCAGCTAAAGCTGCAGTTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTTTGGCTTCa                                                      <  4:125628/149‑1 (MQ=255)                                                                                             ggCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGc                                                    >  2:29601/1‑149 (MQ=255)                                                                                             ggCTTGGGAGGCTTTAGTTTGTTGGGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGc                                                    <  7:219443/149‑1 (MQ=255)                                                                                              gCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAagcagc                                                                                            >  8:35677/1‑108 (MQ=255)                                                                                              gCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAagcagc                                                                                            <  7:35677/108‑1 (MQ=255)                                                                                              gCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGa                                                                 <  6:235823/135‑1 (MQ=255)                                                                                              gCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGa                                                                 >  5:235823/1‑135 (MQ=255)                                                                                              gCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGcc                                                   <  2:186684/149‑1 (MQ=255)                                                                                               cTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGcc                                                   >  2:227214/1‑148 (MQ=255)                                                                                                 ggggAGGCTTTAGTGTGTTGTGGTGCTAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCaa                                                 <  4:11853/147‑1 (MQ=255)                                                                                                  gggAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAg                                                                >  5:70677/1‑132 (MQ=255)                                                                                                  gggAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAg                                                                <  4:163372/132‑1 (MQ=255)                                                                                                  gggAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAg                                                                <  6:70677/132‑1 (MQ=255)                                                                                                  gggAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAg                                                                >  3:163372/1‑132 (MQ=255)                                                                                                   ggAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATAc                                              <  1:227214/149‑1 (MQ=255)                                                                                                    gAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGag                                                                                                                     >  4:247379/1‑77 (MQ=255)                                                                                                    gAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGag                                                                                                                     <  3:247379/77‑1 (MQ=255)                                                                                                    gAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATAcc                                             >  5:293779/1‑149 (MQ=255)                                                                                                              tttGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTgg                                   <  4:106469/149‑1 (MQ=255)                                                                                                                     tgGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCcgcg                            >  7:70063/1‑149 (MQ=255)                                                                                                                        gTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGCGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCcgcgcg                          >  4:18267/1‑148 (MQ=255)                                                                                                                           gAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCtt                       <  8:70063/148‑1 (MQ=255)                                                                                                                              gCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGagcagc                                                                                                                 <  6:131814/55‑1 (MQ=255)                                                                                                                              gCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGagcagc                                                                                                                 >  5:131814/1‑55 (MQ=255)                                                                                                                                  ccgaccTGCAGCTTGAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTAct               >  1:256316/5‑149 (MQ=255)                                                                                                                                   aCGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActc              >  6:63361/1‑149 (MQ=255)                                                                                                                                       cTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcgctc          >  6:79177/1‑149 (MQ=255)                                                                                                                                       cTGCAGCTCAAGCTGCAGGGCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCATCGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGATTCCGGTActctctc          <  1:254620/149‑1 (MQ=255)                                                                                                                                          cAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcgctcg         <  5:114431/147‑1 (MQ=255)                                                                                                                                             cTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGcgccg                              <  1:255944/123‑1 (MQ=255)                                                                                                                                             cTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGcgccg                              >  2:255944/1‑123 (MQ=255)                                                                                                                                             cTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAg    >  5:135844/1‑149 (MQ=255)                                                                                                                                               cAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAggg  >  6:64059/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     GCAGGCATGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGG  >  NC_002947/160168‑160458

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