Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,761,039	C→A	6.5%	L63I (CTT→ATT)	PP_1567 →	capsid protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,761,039	0	C	A	6.5%	288.7 / 6.2	89	L63I (CTT→ATT)	PP_1567	capsid protein
Reads supporting (aligned to +/- strand):  ref base C (40/43);  new base A (4/2);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

CGCAAGGAGCACGCCGACGACAAGATCCAGCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCT  >  NC_002947/1760892‑1761187                                                                                                                                                     |                                                                                                                                                     cgcAAGGAGCACGCCGACGACAAGATCCAGCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCt                                                                                                                                                     <  2:48973/149‑1 (MQ=255)      ggAGCACGCCGACGACAAGATCCAGCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATg                                                                                                                                                >  5:37369/1‑149 (MQ=255)       gAGGACGCCGTCTACTAGATCCATTTGAATGCGGTCAATGTCGTGCTGGAGCAGCGTGACTAGGTGATGTAGGCATTCGCTAAGTAGGCGGCTGACGAGATCAAGGAGCACTGGAACATGCTGTCC‑GAAACCAAGACCATCATTGTTgg                                                                                                                                               <  4:302430/149‑1 (MQ=255)       gAGCACGCCGACGACAAGATCCAGCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATgg                                                                                                                                               >  7:292455/1‑149 (MQ=255)         gCACGCCGACGACAAGATCCAGCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGcc                                                                                                                                             <  3:68743/149‑1 (MQ=255)         gCACGCCGACGACAAGATCCAGCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGcc                                                                                                                                             <  2:32962/149‑1 (MQ=255)           aCGCCGACGACAAGATCCAGCTGAAAGCGGTCATTGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCtt                                                                                                                                           <  2:341438/149‑1 (MQ=255)                 acgacAAGATCCAGCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAg                                                                                                                                     <  4:73688/149‑1 (MQ=255)                     cAAGATCCAGCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGAt                                                                                                                                                 <  4:190983/133‑1 (MQ=255)                     cAAGATCCAGCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGAt                                                                                                                                                 >  3:190983/1‑133 (MQ=255)                     cAAGATCCAGCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATg                                                                                                                                 <  7:222734/149‑1 (MQ=255)                         aTCCAGCTGAAAGCGGTCATTGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGTCCTTGAGAAGGAAGGCCt                                                                                                                             >  3:131046/1‑149 (MQ=255)                         aTCCAGCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGcca                                                                                                                             >  5:105345/1‑148 (MQ=255)                           ccAGCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTgg                                                                                                                           <  5:127866/149‑1 (MQ=255)                            cAGCTGAACGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTgg                                                                                                                           <  8:6296/148‑1 (MQ=255)                              gCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTc                                                                                                                        <  6:37369/149‑1 (MQ=255)                              gCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTc                                                                                                                        >  6:64082/1‑149 (MQ=255)                              gCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTc                                                                                                                        >  4:21070/1‑149 (MQ=255)                                  aaaGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGc                                                                                                                     <  7:187285/148‑1 (MQ=255)                                  aaaGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGc                                                                                                                     >  4:3677/1‑148 (MQ=255)                                       ggTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGc                                                                                                                     >  4:336062/1‑143 (MQ=255)                                        gTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCg                                                                                                              >  4:255497/1‑149 (MQ=255)                                           aaTGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCt                                                                                                           <  5:64082/149‑1 (MQ=255)                                               aCCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCaa                                                                                                       <  3:350336/149‑1 (MQ=255)                                                ccTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCaa                                                                                                       <  8:292455/148‑1 (MQ=255)                                                 cTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTtg                                                                                                                                          <  7:224103/112‑1 (MQ=255)                                                 cTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTtg                                                                                                                                          >  8:224103/1‑112 (MQ=255)                                                   gCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTAGTGAAGGATGGCCTCGGTCTGCAAGATCGCCTGCaa                                                                                                       >  4:67181/1‑145 (MQ=255)                                                       gACCTGCGTGACGAGGAGATCAAGGGATTCGCATAGAAGTCGGCTGACGAGATCAAGGAGGACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCGTTGTGAAGGATGGCCTGGGTCTTCAAGATCGCCTGCAAGAGATc                                                                                                 <  4:131046/147‑1 (MQ=255)                                                       gACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATc                                                                                                 >  3:13014/1‑147 (MQ=255)                                                       gACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGa                                                                                               <  3:327280/149‑1 (MQ=255)                                                        aCCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAg                                                                                              <  3:3677/149‑1 (MQ=255)                                                        aCCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAg                                                                                              <  5:155977/149‑1 (MQ=255)                                                               tGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCagaaga                                                                                       <  3:21070/149‑1 (MQ=255)                                                                 aCGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATg                                                                                     <  5:170561/149‑1 (MQ=255)                                                                   gaggagATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGc                                                                                   >  3:333914/1‑149 (MQ=255)                                                                         aTCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCg                                                                             <  3:90470/149‑1 (MQ=255)                                                                              ggCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAgaga                                                                                                   >  6:153162/1‑122 (MQ=255)                                                                              ggCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAgaga                                                                                                   <  5:153162/122‑1 (MQ=255)                                                                               gCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGAt                                                                                      <  5:307721/134‑1 (MQ=255)                                                                               gCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGAt                                                                                      >  6:307721/1‑134 (MQ=255)                                                                               gCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTtctc                                                                       >  2:16369/1‑149 (MQ=255)                                                                                 aTTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCCGGAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTcc                                                                      >  1:101648/1‑149 (MQ=255)                                                                                 aTTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCg                                                                     >  6:64080/1‑149 (MQ=255)                                                                                  ttCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGc                                                                    <  3:62382/149‑1 (MQ=255)                                                                                     gCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCaa                                                                 >  5:181241/1‑149 (MQ=255)                                                                                     gCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCa                                                                  >  8:161410/1‑148 (MQ=255)                                                                                     gCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCACTTCTCCGCCaa                                                                 >  8:212135/1‑149 (MQ=255)                                                                                      ccAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAAt                                                                >  6:231977/1‑149 (MQ=255)                                                                                       cAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAAt                                                                >  7:18037/1‑148 (MQ=255)                                                                                        aagaagGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCCCTTCTCCGCCAATg                                                               >  3:112981/1‑148 (MQ=255)                                                                                           aagGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATc                                                            >  2:227572/1‑148 (MQ=255)                                                                                              gcggcTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGt                                                        >  4:233743/1‑149 (MQ=255)                                                                                               cggcTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGtt                                                       <  5:64080/149‑1 (MQ=255)                                                                                                ggcTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGtt                                                       >  3:271162/1‑148 (MQ=255)                                                                                                 gcTGACGAGATCAAGGAGCACGGCACCTTCCTGGCC‑GAAACCAAGACCATCCTTTATTGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCGAGTTGa                                                     <  3:262352/149‑1 (MQ=255)                                                                                                  cTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGAc                                                    >  4:24568/1‑149 (MQ=255)                                                                                                     aCGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACtt                                                  >  7:29631/1‑148 (MQ=255)                                                                                                            cAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCgg                                          >  5:105787/1‑149 (MQ=255)                                                                                                            cAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATAGCCTGc                                                                                                         >  4:76716/1‑86 (MQ=255)                                                                                                             aaGGAGCACGGCACCATCCTGGCC‑GAAACCAAGAGCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCGAATGATCCAGTTGACTTCAAGTCgg                                          <  8:154253/148‑1 (MQ=255)                                                                                                              aGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTTGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGGCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCt                                        <  7:161410/149‑1 (MQ=255)                                                                                                               ggAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTg                                       <  6:181241/149‑1 (MQ=255)                                                                                                               ggAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGACCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTg                                       <  4:13014/149‑1 (MQ=255)                                                                                                                    aCGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCg                                    >  6:277922/1‑147 (MQ=255)                                                                                                                     cGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCgaag                                 <  3:24568/149‑1 (MQ=255)                                                                                                                      ggCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCCCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCgaaga                                <  7:212135/149‑1 (MQ=255)                                                                                                                           gaTCCTGTCC‑GAAACCAAGACGATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGa                           <  4:112981/148‑1 (MQ=255)                                                                                                                           cATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGa                           >  4:202487/1‑149 (MQ=255)                                                                                                                            aTCCTGGGG‑GAAACCAAGACGATCCTTGATGGGCTGGTGAGGGATGGCCTGGGTCTGCAGGATCGCCTGCAAGAGATCGCGCAGAAGGTGGCCCGTCGCTTCTCGGGCAATGATCGGGTTGACTGCAGGTCGGCTGGCGAGGACCTggc                          >  4:343136/1‑147 (MQ=255)                                                                                                                            aTCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGAc                          <  8:29631/149‑1 (MQ=255)                                                                                                                            aTCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCAGGCGAAGAACTGAc                          <  6:105787/149‑1 (MQ=255)                                                                                                                                  gCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAgt                     <  4:333914/148‑1 (MQ=255)                                                                                                                                  gCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAgt                     <  3:255497/148‑1 (MQ=255)                                                                                                                                  gCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAgt                     <  2:101648/148‑1 (MQ=255)                                                                                                                                   cc‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAgtgt                   <  3:233743/149‑1 (MQ=255)                                                                                                                                    tccttcacAAGGCCATCAAGGATGGTCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAgtgt                   >  4:328925/8‑149 (MQ=255)                                                                                                                                    c‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGtgtg                  <  3:67181/149‑1 (MQ=255)                                                                                                                                      gAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTga                 >  4:365835/1‑149 (MQ=255)                                                                                                                                       ttcacAAGGCCATCAAGGATGGTCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTgac                >  7:248221/5‑149 (MQ=255)                                                                                                                                       ttcacAAGGCCATCAAGGATGGTCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTgac                >  2:11905/5‑149 (MQ=255)                                                                                                                                       ttcacAAGGCCATCAAGGATGGTCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTACGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGtgtg                  >  5:44585/5‑147 (MQ=255)                                                                                                                                       aaaCCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTgac                >  3:92821/1‑149 (MQ=255)                                                                                                                                         aCCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTgacga              >  5:99838/1‑149 (MQ=255)                                                                                                                                           cAAGGCCATCAAGGATGGTCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACt            <  1:11905/149‑1 (MQ=255)                                                                                                                                               aCCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCaa        <  3:336062/149‑1 (MQ=255)                                                                                                                                                 cATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTgacgac             >  5:316361/1‑142 (MQ=255)                                                                                                                                                 cATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTgacgac             <  6:316361/142‑1 (MQ=255)                                                                                                                                                   tCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCg    >  3:146545/1‑149 (MQ=255)                                                                                                                                                     cTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCt  <  1:389945/149‑1 (MQ=255)                                                                                                                                                     cTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCt  >  1:156124/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CGCAAGGAGCACGCCGACGACAAGATCCAGCTGAAAGCGGTCAATGACCTGCTCGACCAGCGTGACGAGGAGATCAAGGCATTCGCCAAGAAGGCGGCTGACGAGATCAAGGAGCACGGCACCATCCTGGCC‑GAAACCAAGACCATCCTTGATGGCCTTGTGAAGGATGGCCTGGGTCTGCAAGATCGCCTGCAAGAGATCGAGCAGAAGATGGCCCGTCGCTTCTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCT  >  NC_002947/1760892‑1761187

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