Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,761,928	A→T	7.5%	N359I (AAC→ATC) ‡	PP_1567 →	capsid protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,761,928	0	A	T	7.5%	236.8 / 7.7	80	N359I (AAC→ATC) ‡	PP_1567	capsid protein
Reads supporting (aligned to +/- strand):  ref base A (37/37);  new base T (4/2);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

GCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCACGGC  >  NC_002947/1761786‑1762076                                                                                                                                                |                                                                                                                                                      gCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑cgc                                                                                                                                                <  7:81734/149‑1 (MQ=255)   ggACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGTTCAT‑cgcgc                                                                                                                                              <  8:241231/149‑2 (MQ=255)   ggACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCgg                                                                                                                                              >  8:279504/1‑149 (MQ=255)       aTCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGcc                                                                                                                                            >  2:139764/1‑147 (MQ=255)        tCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAg                                                                                                                                         <  5:96585/149‑1 (MQ=255)               gACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGTTCAT‑cgcgc                                                                                                                                              >  6:23319/1‑136 (MQ=255)                aCCAAGGACGCCAACGTCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGTAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTc                                                                                                                                   <  4:388348/147‑1 (MQ=255)                aCCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTc                                                                                                                                   <  8:39665/147‑1 (MQ=255)                aCCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTc                                                                                                                                   <  7:74222/147‑1 (MQ=255)                  cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGTTCAT‑cgcgcctacca                                                                                                                                        <  6:19636/140‑8 (MQ=255)                  cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGTTCAT‑cgcgcctacca                                                                                                                                        >  5:19636/1‑133 (MQ=255)                  cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGAcc                                                                                                                               >  8:48428/1‑149 (MQ=255)                  cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGAcc                                                                                                                               >  2:64926/1‑149 (MQ=255)                  cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGAc                                                                                                                                >  5:72429/1‑148 (MQ=255)                  cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGAc                                                                                                                                <  6:72429/148‑1 (MQ=255)                    aGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACcgc                                                                                                                             >  1:407896/1‑149 (MQ=255)                           cAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAg                                                                                                                                         >  1:294345/1‑131 (MQ=255)                           cAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAg                                                                                                                                         <  2:294345/130‑1 (MQ=255)                            aaCGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGc                                                                                                                     <  2:231643/149‑1 (MQ=255)                            aaCGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGc                                                                                                                     >  6:285337/1‑149 (MQ=255)                             aCGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCg                                                                                                                    <  5:285337/149‑1 (MQ=255)                                  tCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGAGGCTGt                                                                                                               >  7:240700/1‑149 (MQ=255)                                       cATCTGGGTGAATGTCCAGGATGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAg                                                                                                           <  5:207502/148‑1 (MQ=255)                                              gTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCAc                                                                                                   <  8:265840/149‑1 (MQ=255)                                              gTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCAc                                                                                                   >  3:35407/1‑149 (MQ=255)                                                  aTGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGaa                                                                                               <  7:269846/149‑1 (MQ=255)                                                     tCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGAc                                                                                            >  8:35283/1‑149 (MQ=255)                                                     tCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGAc                                                                                            >  5:248985/1‑149 (MQ=255)                                                     tCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGAc                                                                                            >  6:256888/1‑149 (MQ=255)                                                       cAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACgg                                                                                          >  6:198093/1‑149 (MQ=255)                                                        aGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGt                                                                                         >  4:242117/1‑149 (MQ=255)                                                        aGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGt                                                                                         <  2:407896/149‑1 (MQ=255)                                                         ggAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTg                                                                                        <  6:18775/149‑1 (MQ=255)                                                          gAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGa                                                                                       >  5:223107/1‑149 (MQ=255)                                                           aaGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGAc                                                                                      <  5:256888/149‑1 (MQ=255)                                                                 gcgTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGGTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGAGCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTc                                                                                <  7:249396/149‑1 (MQ=255)                                                                             cTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCACGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGa                                                                                       >  6:32485/1‑131 (MQ=255)                                                                             cTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCACGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGa                                                                                       <  5:32485/131‑1 (MQ=255)                                                                                tgGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATg                                                                    <  5:114325/147‑1 (MQ=255)                                                                                tgGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGt                                                                  >  6:232098/1‑149 (MQ=255)                                                                                  gAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGa                                                                >  6:74649/1‑149 (MQ=255)                                                                                  gAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGa                                                                >  7:46931/1‑149 (MQ=255)                                                                                          cGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCcgc                                                        >  6:271963/1‑149 (MQ=255)                                                                                             tggtggACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGcc                                                     <  7:279504/149‑1 (MQ=255)                                                                                                tggACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgag                                                  >  7:70100/1‑149 (MQ=255)                                                                                                 ggACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgag                                                  >  8:97971/1‑148 (MQ=255)                                                                                                  gACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgagga                                                <  5:198093/149‑1 (MQ=255)                                                                                                  gACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgag                                                  >  2:152119/1‑147 (MQ=255)                                                                                                   aCAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgaggag                                               >  6:342113/1‑149 (MQ=255)                                                                                                      gCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAgcgc                                            <  6:223107/149‑1 (MQ=255)                                                                                                          cgcgGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctgg                                        <  1:64926/149‑1 (MQ=255)                                                                                                             ggTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctggcgc                                     <  5:271963/149‑1 (MQ=255)                                                                                                              gTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctggcgct                                    <  1:152119/149‑1 (MQ=255)                                                                                                                 ccAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATGCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctggcgctgg                                  >  1:60208/1‑149 (MQ=255)                                                                                                                  cAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctggcgctggcg                                <  5:256679/149‑1 (MQ=255)                                                                                                                   aGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCgg                               >  1:202968/1‑149 (MQ=255)                                                                                                                   aGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGAGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctgg                                        <  4:52804/140‑1 (MQ=255)                                                                                                                   aGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGAGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctgg                                        >  3:52804/1‑140 (MQ=255)                                                                                                                          gAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCg                        >  8:339999/1‑149 (MQ=255)                                                                                                                          gAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCg                        >  8:81449/1‑149 (MQ=255)                                                                                                                          gAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGCAGCGCCTGGCGCTGGCGGTCTACCg                        >  1:321495/1‑149 (MQ=255)                                                                                                                               ccTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCa                    <  8:240700/148‑1 (MQ=255)                                                                                                                               ccTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCa                    <  7:97971/148‑1 (MQ=255)                                                                                                                                cTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTAGACCGccagcg                  <  6:115740/149‑5 (MQ=255)                                                                                                                                 tGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCagag                 >  8:63108/1‑149 (MQ=255)                                                                                                                                  ctgggccgcgATGTTCAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGt                >  3:376822/7‑149 (MQ=255)                                                                                                                                   gTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCgcgc                                                      >  3:124878/1‑110 (MQ=255)                                                                                                                                   gTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCgcgc                                                      <  4:124878/110‑1 (MQ=255)                                                                                                                                      ggCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGtcgtt            <  8:70100/149‑1 (MQ=255)                                                                                                                                       ccgcgATGTTCAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGtcgttc           >  6:256953/2‑149 (MQ=255)                                                                                                                                        cgcgATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGtcgttcg          <  4:35407/149‑1 (MQ=255)                                                                                                                                         gcgATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTg                                                                 >  8:221018/1‑93 (MQ=255)                                                                                                                                         gcgATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTg                                                                 <  7:221018/93‑1 (MQ=255)                                                                                                                                         gcgATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGAGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCagag                 >  5:202951/1‑141 (MQ=255)                                                                                                                                         gcgATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGAGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCagag                 <  6:202951/141‑1 (MQ=255)                                                                                                                                            aTGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCa      <  7:63108/149‑1 (MQ=255)                                                                                                                                             tGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCa                    <  8:27712/134‑1 (MQ=255)                                                                                                                                             tGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCa                    >  7:27712/1‑134 (MQ=255)                                                                                                                                             tGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCAc     <  7:339999/149‑1 (MQ=255)                                                                                                                                              gAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCACg    >  1:391991/1‑149 (MQ=255)                                                                                                                                                aCAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCACgg   <  2:202968/148‑1 (MQ=255)                                                                                                                                                aCAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCACGGc  <  7:35283/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                      GCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACA‑GCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAT‑CGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCACGGC  >  NC_002947/1761786‑1762076

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