Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,432,011	G→A	6.6%	V203M (GTG→ATG) ‡	PP_4771 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,432,011	0	G	A	6.6%	221.9 / 3.8	76	V203M (GTG→ATG) ‡	PP_4771	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (27/44);  new base A (2/3);  total (29/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-01

CATGCGCATCATGATCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTATTCAAGGAGTGCACCATGAACGAAAC  >  NC_002947/5431866‑5432157                                                                                                                                                  |                                                                                                                                                    cATGCGCATCATGATCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCgtg                                                                                                                                                   <  7:10333/148‑1 (MQ=255)  aTGCGCATCATGATCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCgtggt                                                                                                                                                 >  2:63772/1‑149 (MQ=255)     cgcATCATGATCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTTAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCt                                                                                                                                              <  2:100467/149‑1 (MQ=255)     cgcATCATGATCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGc                                                                                                                                               <  1:295627/148‑1 (MQ=255)     cgcATCATGATCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTACAGCACCACGGTGc                                                                                                                                               <  7:136242/148‑1 (MQ=255)     cgcATCATGATCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGACCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGc                                                                                                                                               <  1:363769/148‑1 (MQ=255)         tcatGATCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTAc                                                                                                                                          <  1:211563/149‑1 (MQ=255)           atGATCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCg                                                                                                                                        >  5:10744/1‑149 (MQ=255)           atGATCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCg                                                                                                                                        >  5:248283/1‑149 (MQ=255)           atGATCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCg                                                                                                                                        >  1:21925/1‑149 (MQ=255)              aTCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTACAGCACCACGGTGCTGAACAGctg                                                                                                                                     <  7:137218/149‑3 (MQ=255)               tCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCACGGTGCTGAAc                                                                                                                                          <  5:64050/143‑1 (MQ=255)               tCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCACGGTGCTGAAc                                                                                                                                          >  6:64050/1‑143 (MQ=255)                  gCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGtt                                                                                                                                  >  2:62386/1‑148 (MQ=255)                   cGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCTTGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTc                                                                                                                                 <  3:307544/148‑1 (MQ=255)                   cGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGt                                                                                                                                   <  8:101055/146‑1 (MQ=255)                   cGTCTTAGGCGGGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCa                                                                                                                                <  7:168510/149‑1 (MQ=255)                    gTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCAt                                                                                                                               <  1:107404/149‑1 (MQ=255)                    gTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCAt                                                                                                                               <  5:301160/149‑1 (MQ=255)                     tCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATg                                                                                                                              >  6:41004/1‑149 (MQ=255)                     tCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATg                                                                                                                              <  6:10744/149‑1 (MQ=255)                         aGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCAAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGc                                                                                                                          >  1:218761/1‑149 (MQ=255)                          ggCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGc                                                                                                                          <  5:378333/148‑1 (MQ=255)                          ggCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGc                                                                                                                          <  5:228803/148‑1 (MQ=255)                              ttcTGGCGCCGTTCTGGGTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCGTCAGGGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGaa                                                                                                                     <  2:377659/147‑1 (MQ=255)                              tGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGaa                                                                                                                     <  5:41004/149‑1 (MQ=255)                              tGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGaa                                                                                                                     <  5:229761/149‑1 (MQ=255)                                     gcCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTACATGTTGCTGGAAGGTGAAg                                                                                                              >  8:74449/1‑149 (MQ=255)                                      ccGATCTGGCTCATCGACGGGCTGCTGTTGATGACCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGc                                                                                                             <  5:154115/149‑1 (MQ=255)                                       cGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGc                                                                                                             >  1:68294/1‑148 (MQ=255)                                          tCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCAAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAg                                                                                                          <  2:218761/148‑1 (MQ=255)                                              gCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCa                                                                                                           <  1:63772/143‑1 (MQ=255)                                              gCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCAc                                                                                                       <  2:156412/147‑1 (MQ=255)                                                  aTCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATg                                                                                                                              <  8:98623/120‑1 (MQ=255)                                                  aTCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATg                                                                                                                              >  7:98623/1‑120 (MQ=255)                                                  aTCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGa                                                                                                   <  4:319186/147‑1 (MQ=255)                                                  aTCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGa                                                                                                   <  4:324547/147‑1 (MQ=255)                                                   tCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGc                                                                                                >  8:205388/1‑149 (MQ=255)                                                    cGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCgtg                                                                                                                                                   >  6:300788/1‑97 (MQ=255)                                                    cGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCgtg                                                                                                                                                   <  5:300788/97‑1 (MQ=255)                                                          gctgctGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCggg                                                                                         >  5:187664/1‑149 (MQ=255)                                                            tgctgTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCAGCAGCTGATCAGCACCGTGGTGCTGTACCGCAAGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCgggg                                                                                        >  2:269501/1‑148 (MQ=255)                                                              ctgTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAAt                                                                                     >  5:330684/1‑149 (MQ=255)                                                                    aTGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGtttt                                                                                <  2:68294/148‑1 (MQ=255)                                                                    aTGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTc                                                                               <  8:104784/149‑1 (MQ=255)                                                                    aTGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTc                                                                               >  3:299320/1‑149 (MQ=255)                                                                       gCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTGCAGCACCGTGGGGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGAc                                                                            <  8:170520/149‑1 (MQ=255)                                                                        ccTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTc                                                                               <  3:40278/145‑1 (MQ=255)                                                                        ccTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTc                                                                               >  4:40278/1‑145 (MQ=255)                                                                         cTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCt                                                                          <  2:21925/149‑1 (MQ=255)                                                                             cccGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCg                                                                      >  7:84531/1‑149 (MQ=255)                                                                                  gacATCACCTGGGATGGAACTGTCCTTTGAACGCCTGAGCGGCTTCCTGCAGCTGTTCAGCACCGTTGTTCTGTACCTCCTGTTCATGTTGCTTGAAGGTGTAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGc                                                                 <  8:251159/147‑1 (MQ=255)                                                                                   ccATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGcc                                                                >  2:104610/1‑149 (MQ=255)                                                                                     aTCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCa                                                              <  3:271947/149‑1 (MQ=255)                                                                                      tCACCTGGGATGGAACTGGACATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAg                                                             >  8:253898/1‑149 (MQ=255)                                                                                        aCCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGcc                                                           >  4:366515/1‑149 (MQ=255)                                                                                         ccTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGccc                                                          >  7:65866/1‑149 (MQ=255)                                                                                               aTGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGgcg‑c                                                   <  4:316865/149‑1 (MQ=255)                                                                                                  gAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGcc                                                           <  3:227683/139‑1 (MQ=255)                                                                                                  gAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGcc                                                           >  4:227683/1‑139 (MQ=255)                                                                                                   aaCTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑Ccaca                                               <  3:366515/149‑1 (MQ=255)                                                                                                   aaCTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTACCCAGCCCTGGCG‑Cca                                                 >  4:296999/1‑147 (MQ=255)                                                                                                     cTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGt                                             <  4:299320/149‑1 (MQ=255)                                                                                                               aCAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTACATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTc                                   <  7:74449/149‑1 (MQ=255)                                                                                                                    cTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGc                              >  6:96257/1‑149 (MQ=255)                                                                                                                       aGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTa                            <  3:304991/148‑1 (MQ=255)                                                                                                                          ggCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGTGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTATTCa                        <  7:253898/149‑1 (MQ=255)                                                                                                                            cTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTATTCAAg                      <  8:65866/149‑1 (MQ=255)                                                                                                                            cTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTATTCAAg                      <  8:84531/149‑1 (MQ=255)                                                                                                                             ttCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTATTCAAg                      >  5:338849/1‑148 (MQ=255)                                                                                                                                  gCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTATTCAAGGAGTg                 >  1:100782/1‑148 (MQ=255)                                                                                                                                  gCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTATTCAAGGAGTg                 >  2:79542/1‑148 (MQ=255)                                                                                                                                   cAGGTGTTCAGGACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCATCAACTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCGCCCACAGTGCACATCCTCAATGCTATTCAAGGAGTGc                <  1:62386/149‑1 (MQ=255)                                                                                                                                        gtacAGCACCACGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTATTCAAGGAGTGCACCAt           >  5:83919/4‑148 (MQ=255)                                                                                                                                         ttCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTAt                           <  6:167402/131‑1 (MQ=255)                                                                                                                                         ttCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTAt                           >  5:167402/1‑131 (MQ=255)                                                                                                                                               aCCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTATTCAAGGAGTGCACCATGAACGaaa   <  3:178276/149‑1 (MQ=255)                                                                                                                                               aCCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTATTCAAGGAGTGCACCATGAACGaaa   <  6:135806/149‑1 (MQ=255)                                                                                                                                               aCCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTACCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTATTCAAGGAGTGCACCATGAACGaaa   <  3:296999/149‑1 (MQ=255)                                                                                                                                                ccGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTATTCAAGGAGTGCACCATGAACGAAAc  <  6:330684/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CATGCGCATCATGATCTGCGTCTTAGGCGTGCTGGCGCCGATCTGGCTCATCGACGGGCTGCTGTTGATGGCCTTCCCCGAGCCATCACCTGGGATGGAACTGGCCATGGACAGCCTCAGCGGCTTCCTGCAGCTGTTCAGCACCGTGGTGCTGTACCGCCTGTTCATGTTGCTGGAAGGTGAAGCCAGCACCTGAGGCCACCGGGGAATGTTTTCGACCTGCCGGCTGCCCAGCCCTGGCG‑CCACAGTGCACATCCTCAATGCTATTCAAGGAGTGCACCATGAACGAAAC  >  NC_002947/5431866‑5432157

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