Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	501,338	G→T	5.3%	S405I (AGC→ATC)	gsk →	inosine/guanosine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,338	0	G	T	5.3%	302.4 / 2.2	89	S405I (AGC→ATC)	gsk	inosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base G (44/40);  new base T (2/3);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTG  >  NC_000913/501193‑501475                                                                                                                                                  |                                                                                                                                          ggggATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTa                                                                                                                                        >  6:13042/1‑149 (MQ=255)    gATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTAt                                                                                                                                       >  4:256754/1‑147 (MQ=255)        gcgcATTGGCAGCGTTGCTGCATGACATTACCTCCATCAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGTTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTa                                                                                                                                 <  3:394120/149‑1 (MQ=255)        gcgcATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCGGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTa                                                                                                                                 <  2:35297/149‑1 (MQ=255)                gCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTTACCGTGTGAGCTATCAGGTACTGAAc                                                                                                                           >  3:334506/1‑147 (MQ=255)                gCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAAc                                                                                                                           >  7:262987/1‑147 (MQ=255)                gCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCa                                                                                                                         <  5:415692/149‑1 (MQ=255)                     gTTGCTTCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCAtt                                                                                                                    <  4:375911/149‑1 (MQ=255)                         ctgcATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTAt                                                                                                                                       <  6:235898/126‑1 (MQ=255)                         ctgcATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTAt                                                                                                                                       >  5:235898/1‑126 (MQ=255)                                aCATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTCAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGttt                                                                                                          >  8:84744/1‑148 (MQ=255)                                        gCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTTACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCgg                                                                                                 <  4:334506/149‑1 (MQ=255)                                        gCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGt                                                                                                            >  7:128536/1‑138 (MQ=255)                                        gCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGt                                                                                                            <  8:128536/138‑1 (MQ=255)                                         ccAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGc                                                                                                >  1:55407/1‑149 (MQ=255)                                          cAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCt                                                                                               <  3:198887/149‑1 (MQ=255)                                            aCAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTg                                                                                             >  3:347814/1‑149 (MQ=255)                                             cAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGc                                                                                            <  7:366085/149‑1 (MQ=255)                                              aGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGcc                                                                                           <  2:382323/149‑1 (MQ=255)                                               gatACCATCGTAGCAACGTACCAAACTCCAGCAAAAATAAATTCACCTGTTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCATGTACTGAACAAGCATTCACCTCTTTTAACGCGCGGCTTGCCg                                                                                          >  8:8467/3‑149 (MQ=255)                                               gCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCg                                                                                          >  8:14837/1‑149 (MQ=255)                                               gCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCAGCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCg                                                                                          <  6:250373/149‑1 (MQ=255)                                                cTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCg                                                                                          >  8:347927/1‑148 (MQ=255)                                                 tACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCg                                                                                          >  1:86559/1‑147 (MQ=255)                                                  aCCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAg                                                                                       >  8:377812/1‑149 (MQ=255)                                                  aCCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAg                                                                                       >  7:334461/1‑149 (MQ=255)                                                   ccATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGa                                                                                        <  5:13042/147‑1 (MQ=255)                                                   ccATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTACCCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGc                                                                                      >  6:334106/1‑149 (MQ=255)                                                    cATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCg                                                                                     <  3:256754/149‑1 (MQ=255)                                                     aTCGTAGCAACGTACCAACCTCCAGCATACATAAATTCACCTGGTTAACTTATACATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTGTCCTCGTTTAACGCGCGGCTTCCCGGAGCGt                                                                                    >  6:329557/1‑149 (MQ=255)                                                     aTCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTCAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGt                                                                                    <  7:84744/149‑1 (MQ=255)                                                         tAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAg                                                                                >  4:274818/1‑149 (MQ=255)                                                         tAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAg                                                                                >  2:253042/1‑149 (MQ=255)                                                            cAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATAAGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACa                                                                             >  5:172665/1‑149 (MQ=255)                                                                 tACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATAGTTAGCGAAGGTGTGTAAATATGCTAACAGTGTGAGCTATCAGGTACAGAACCAGCATTCACCTAGTTTAAcgcgc                                                                                                   >  3:80932/1‑122 (MQ=255)                                                                 tACCAAACTCAATCAAACATAATTTCACCTGGTTAAGTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAAcgcgc                                                                                                   <  4:80932/122‑1 (MQ=255)                                                                  aCCAAACTCCAGCAAACTTAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTg                                                                        >  4:187140/1‑148 (MQ=255)                                                                  aCCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTg                                                                        >  5:364274/1‑148 (MQ=255)                                                                    cAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAAc                                                                                                                           >  7:363502/1‑95 (MQ=255)                                                                    cAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAAc                                                                                                                           <  8:363502/95‑1 (MQ=255)                                                                    cAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGaa                                                                     <  8:262987/149‑1 (MQ=255)                                                                      aaCTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGaa                                                                     >  6:345604/1‑147 (MQ=255)                                                                         tCCAGCAAACATAAATTAACCAGGTTAACTTATTCATCGTTTGCGGAGGTGTGTTAATTTGCTAACCGTGTGAGCGATCAGGTACGGCACCAGGATTAACCGAGTTTTACGCGCGGCGTGCCGGAGCGTGAAGACAGCCGGGATGAGTc                                                                >  4:233945/1‑149 (MQ=255)                                                                           cAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCAtt                                                                                                                    <  4:245195/95‑1 (MQ=255)                                                                           cAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCAtt                                                                                                                    >  3:245195/1‑95 (MQ=255)                                                                           cAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTc                                                                <  5:189912/147‑1 (MQ=255)                                                                           cAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTc                                                                <  2:86559/147‑1 (MQ=255)                                                                            aGCAAACATAAATTCACCTGGTTAACTTATTCATCGGTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTCCTGACCCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGGCTTa                                                             <  6:342095/149‑1 (MQ=255)                                                                             gCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTAc                                                            <  2:55407/149‑1 (MQ=255)                                                                              cAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACt                                                           >  2:6277/1‑149 (MQ=255)                                                                               aaaCATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTg                                                          <  7:347927/149‑1 (MQ=255)                                                                                aaCATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTgg                                                         >  8:190702/1‑149 (MQ=255)                                                                                   aTAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAATAGTCTTACTGGGAt                                                      >  6:280841/1‑149 (MQ=255)                                                                                   aTAAATTCACCTGGTTAAATTATTCATCGTTAGCGCAGTTGTGTAAATATGCTAAACGTGTGAGCTATCAGGTAATGCACCCGCATTCACCTCGTTTAACGCGCGGCTTGCAGGAGCGGGAAGACAGCAGGGAAGAGTCTTACTGGGAt                                                      >  6:280840/1‑149 (MQ=255)                                                                                          caacTGGTTAACTTATTCATCGTTAGCGAAGGTGTTTAACTTTTCTAACCGTGTTATCTATCAGTGACTGAACCATCCTTCACCGCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGAGCGTTac                                                <  5:280840/145‑2 (MQ=255)                                                                                          cACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGGACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTaa                                                <  5:280841/148‑1 (MQ=255)                                                                                           aCCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGt                                              >  3:38965/1‑149 (MQ=255)                                                                                           aCCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCATGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGt                                              >  3:145342/1‑149 (MQ=255)                                                                                            ccTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGtt                                             <  6:172665/149‑1 (MQ=255)                                                                                             cTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGTGGCTTGCCGGAGCGTGAAGACAGCCTg                                                                        <  3:111475/121‑1 (MQ=255)                                                                                             cTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGGGGCTTGCCGGAGCGTGAAGACAGCCTg                                                                        >  4:111475/1‑121 (MQ=255)                                                                                                gTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTg                                                                        >  1:261312/1‑118 (MQ=255)                                                                                                gTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTg                                                                        <  2:261312/118‑1 (MQ=255)                                                                                                 ttAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTAtcgt                                        >  3:125106/1‑149 (MQ=255)                                                                                                 ttAACTTATTCATCGTTAGCGCAGGTGGGTAAATATGCTAACCGTGTGAGCTATCAGGTACTTAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTAtcgg                                        <  5:329557/149‑2 (MQ=255)                                                                                                    aCTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCgg                                     >  5:3049/1‑149 (MQ=255)                                                                                                       tATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTc                                  >  6:229149/1‑149 (MQ=255)                                                                                                         ttCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCg                                 <  1:6277/148‑1 (MQ=255)                                                                                                         ttCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCg                                 >  6:341058/1‑148 (MQ=255)                                                                                                          tCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGATCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGc                                                                                      >  7:255128/1‑94 (MQ=255)                                                                                                          tCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGATCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGc                                                                                      >  7:261581/1‑94 (MQ=255)                                                                                                          tCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGATCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGc                                                                                      <  8:255128/94‑1 (MQ=255)                                                                                                          tCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGATCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGc                                                                                      <  8:261581/94‑1 (MQ=255)                                                                                                          tCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTAtc                                          >  2:404220/1‑138 (MQ=255)                                                                                                          tCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTAtc                                          <  1:404220/138‑1 (MQ=255)                                                                                                           cATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTa                               <  7:14837/148‑1 (MQ=255)                                                                                                           cATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAg                              >  4:151299/1‑149 (MQ=255)                                                                                                           cATCGTTAGCGCAGGTGTGTAAATATGATAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCTCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTa                               <  7:8467/148‑1 (MQ=255)                                                                                                            aTCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAgg                             >  3:290577/1‑149 (MQ=255)                                                                                                             tCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGc                            <  7:377812/149‑1 (MQ=255)                                                                                                               gTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCa                          <  6:364274/149‑1 (MQ=255)                                                                                                                ttAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAg                         <  3:187140/149‑1 (MQ=255)                                                                                                                   gcgcAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATa                      <  8:334461/149‑1 (MQ=255)                                                                                                                    cgcAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATaa                     <  7:190702/149‑1 (MQ=255)                                                                                                                     gcAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAg                    <  3:151299/149‑1 (MQ=255)                                                                                                                               aaaTATGCTAACCGTGTGAGCTATCAGGTACTGCACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTCGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCAcgc          >  1:291767/1‑149 (MQ=255)                                                                                                                                aaTATTCTAACCGTGTGAGCTATCAGGTACTGAAACAGAATTCACCTCGTTTAACGCGCGGCTTGACGGAGCGTGAAGACAGCCTGTAATAGTCTTACTTGGAGCGTTAAGTTATCGTCGGTTATTAGGCCAGATAAGTCGTTCActcc         <  4:347814/149‑1 (MQ=255)                                                                                                                                     gCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATc    >  4:175645/1‑149 (MQ=255)                                                                                                                                         aCCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTg  <  1:253042/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                          GGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTG  >  NC_000913/501193‑501475

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