Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,117,004	T→A	26.2%	D73V (GAT→GTT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,004	0	T	A	26.2%	147.2 / 53.5	93	D73V (GAT→GTT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (20/47);  new base A (13/11);  total (33/60)
Fisher's exact test for biased strand distribution p-value = 4.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

TTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCAT  >  NC_000913/4116856‑4117148                                                                                                                                                     |                                                                                                                                                 ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                                                                                  <  2:907148/149‑1 (MQ=255)     aaaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGa                                                                                                                                              >  2:1100712/1‑149 (MQ=255)      aaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGaa                                                                                                                                             <  1:1032969/149‑1 (MQ=255)      aaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCCTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGaa                                                                                                                                             <  1:1127103/149‑1 (MQ=255)        ccGTCCCGTTTTAATGGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACt                                                                                                                                           <  2:190058/149‑1 (MQ=255)        ccGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACt                                                                                                                                           <  1:17025/149‑1 (MQ=255)         cctcccGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTg                                                                                                                                          <  1:326514/144‑1 (MQ=255)         cGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTg                                                                                                                                          >  1:1109670/1‑149 (MQ=255)          gTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGa                                                                                                                                         <  1:785181/149‑1 (MQ=255)          gTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGGGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGa                                                                                                                                         <  2:1325824/149‑1 (MQ=255)           tCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCCGAACTGat                                                                                                                                        >  1:1046406/1‑149 (MQ=255)               gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCg                                                                                                                                    <  2:275557/149‑1 (MQ=255)                ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCgg                                                                                                                                   >  2:1100494/1‑149 (MQ=255)                ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCgg                                                                                                                                   >  1:968563/1‑149 (MQ=255)                   tAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCtt                                                                                                                                >  2:561388/1‑149 (MQ=255)                    aaaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCttt                                                                                                                               <  1:1402023/149‑1 (MQ=255)                       tGCTCGGAGATTTTGGCGGTACGCCGGCCCTGCCAGACAATGGCGTTAACGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGTTATCGGCTTTCGc                                                                                                                            <  2:887275/149‑1 (MQ=255)                       tGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGCTAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTTGTATCGGCTTTCGc                                                                                                                            <  2:1258439/149‑1 (MQ=255)                            gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCa                                                                                                                       >  2:400494/1‑149 (MQ=255)                                 tttCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCt                                                                                                                  <  2:688282/149‑1 (MQ=255)                                   tCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAc                                                                                                                <  2:473013/149‑1 (MQ=255)                                   tCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAc                                                                                                                <  1:561388/149‑1 (MQ=255)                                       cGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACACTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACCTCTACCAGc                                                                                                            >  2:278533/1‑149 (MQ=255)                                           acgacgGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTgg                                                                                                        <  2:745329/149‑1 (MQ=255)                                                  cACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCATACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTg                                                                                                 <  2:439132/149‑1 (MQ=255)                                                  cACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGGCAATAGTGGTTTCACGCTGGTTGGTAATACAGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCGCTTCTACCAGCGTGGCGCTTTg                                                                                                 >  2:1093634/1‑149 (MQ=255)                                                   aCTGCCAGACAATGGCGTTATATATACGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTgg                                                                                                <  1:845324/149‑1 (MQ=255)                                                   aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgg                                                                                                <  1:957301/149‑1 (MQ=255)                                                   aCTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTGCTTTTCCCCGGACAAGAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTgg                                                                                                <  2:152718/149‑1 (MQ=255)                                                     tGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTACCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGCTAGCTGCAAATTGAACGAAACTGATATCGGATTTCGCCAGCACTTCTACCCGAGTGGAGCTTTagga                                                                                              >  2:692456/1‑148 (MQ=255)                                                     tGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTGTTTCCCAGATAATTGTGGTTTCGCGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGt                                                                                              <  2:1165193/149‑1 (MQ=255)                                                       ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTgg                                                                                            >  2:835386/1‑149 (MQ=255)                                                       ccAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTgggggg                                                                                            >  2:1374530/1‑149 (MQ=255)                                                        cAGACAATGGCGTTATGGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGc                                                                                           <  2:1156612/149‑1 (MQ=255)                                                        cAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACCATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACATCTAACCCCGTGGAGCATTgggccgg                                                                                           >  2:1084284/1‑145 (MQ=255)                                                         aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCCGCGTGGAGCGTTGGGTGGcc                                                                                          >  2:470698/1‑149 (MQ=255)                                                          gaacattgCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGccc                                                                                         <  1:1417319/142‑1 (MQ=255)                                                          gACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGccc                                                                                         <  2:1166693/149‑1 (MQ=255)                                                            cAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTCCCAGCGTGGAGCTTGGGGTGGGCCAg                                                                                       >  1:34120/1‑149 (MQ=255)                                                             aaTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCAAgc                                                                                      <  2:1109670/149‑2 (MQ=255)                                                                  cGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTcc                                                                                 <  1:1220030/149‑1 (MQ=255)                                                                   gTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCa                                                                                >  1:1408276/1‑149 (MQ=255)                                                                    ttATAGATATGCCTGCTCGTTTTTTTTTCCCAGTCACTAGTGGTGTCACGCTGGTTCGTAATACCGCTAGCTGCAATTTGTTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAt                                                                               <  2:765452/149‑1 (MQ=255)                                                                     tataGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTTGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCAtt                                                                              <  2:968563/149‑1 (MQ=255)                                                                      atagatagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCGATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGACTTCCATTg                                                                             >  1:851994/1‑149 (MQ=255)                                                                       tagttagGCCTGCCGGTTGCTTTTTCCCAGACAATAGTGGTTTCACGCTCGTTCGTAATACCGATAGCTGCAATTTGACCGGAACCGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTgg                                                                            <  2:47731/149‑1 (MQ=255)                                                                                ttGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCCCGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTc                                                                   <  1:190352/149‑1 (MQ=255)                                                                                ttGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTc                                                                   <  1:1100494/149‑1 (MQ=255)                                                                                 tGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGCGGCCCAGATTTCCATTGGGTCGTGTTCt                                                                  <  1:1088572/149‑1 (MQ=255)                                                                                   ccGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTTTTCTAc                                                                <  2:1451222/149‑1 (MQ=255)                                                                                   ccGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAc                                                                <  1:49539/149‑1 (MQ=255)                                                                                     ggTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAccc                                                              <  1:1022191/149‑1 (MQ=255)                                                                                     ggTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAccc                                                              <  1:1100712/149‑1 (MQ=255)                                                                                        ttCTTTTTCCCCGACTATAGTGGTTTCACGCTGGTTCCTAATACCGATAGCGTCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCTGATTTCCATTGGGTCGTGTTCTACCCAAc                                                           <  2:478078/149‑1 (MQ=255)                                                                                        ttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAc                                                           <  2:548636/149‑1 (MQ=255)                                                                                           tttttCTCAGACAATAGTGGTTTCACGCCGGAAAGTAATACCGCTAGCTGAAAATTGAACGGAAATGATATCTGCTTTAGACAGCAATCCTGACAGCGCGGAGCTTAGGGTGGCCCAGAGTTCCATTGGGGCGTGTGCTAGCCAGCCTg                                                        >  2:622455/1‑149 (MQ=255)                                                                                           tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATGTCGGCTTTCGCCAGCACTTGTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTGGTGTTCTAGCCCACCTg                                                        >  2:629386/1‑149 (MQ=255)                                                                                           tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                        <  1:835386/149‑1 (MQ=255)                                                                                              ttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtt                                                     <  1:200755/149‑1 (MQ=255)                                                                                                cccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGAAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                   <  2:782875/149‑1 (MQ=255)                                                                                                  cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTCGGGTGGCCCAGATTTCCATTGGGTCGTGTTCCACCCAACCCGGTTTTgg                                                 >  1:316771/1‑149 (MQ=255)                                                                                                  cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGAGAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTgg                                                 <  2:654697/149‑1 (MQ=255)                                                                                                   aGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTggg                                                >  2:794559/1‑149 (MQ=255)                                                                                                    gACAATAGTGGTTTCACGCTGGTTCGTAATGCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGt                                               <  2:1014057/149‑1 (MQ=255)                                                                                                     aCAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTa                                              <  1:1216789/149‑1 (MQ=255)                                                                                                       aaTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAgg                                            >  1:1142321/1‑148 (MQ=255)                                                                                                       aaTAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCCACCTGGTCTTGGGTAGa                                            >  1:1210157/1‑149 (MQ=255)                                                                                                            tGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGc                                       >  1:592995/1‑149 (MQ=255)                                                                                                             ggTTTCCCGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCt                                      <  1:1432395/149‑1 (MQ=255)                                                                                                             ggTTTCACGCTGGTTCGTAGTACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCt                                      <  2:823385/149‑1 (MQ=255)                                                                                                             ggTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCt                                      <  2:143494/149‑1 (MQ=255)                                                                                                              gTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCGCTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTCCCCAACCTGGTTTTGGGTAGATTTGCTc                                     >  1:846691/1‑149 (MQ=255)                                                                                                               tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCGAGCTGGTTCTGGGATGATTTGCTCa                                    >  1:1246616/1‑149 (MQ=255)                                                                                                                ttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCGCTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGGCGCTTTCTTCCCAAGCTGGTTTTGGGTAGATTTGCTCaa                                   >  1:1322335/1‑149 (MQ=255)                                                                                                                 tCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaaa                                  <  1:1190843/149‑1 (MQ=255)                                                                                                                 tCACGCGGGGTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaaa                                  <  2:383906/149‑1 (MQ=255)                                                                                                                          ttCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTg                         >  1:276368/1‑149 (MQ=255)                                                                                                                          gtCGTAATACCGATAGCTGCAATTGGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTg                         <  2:348834/148‑1 (MQ=255)                                                                                                                            cGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCGAATTCGCGCTGCg                       <  1:993533/149‑1 (MQ=255)                                                                                                                              tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGATTTGGGTAGATTTTCTCAAATTCGCGATGCGac                     >  2:841206/1‑149 (MQ=255)                                                                                                                              tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTTCCGCTGCCac                     >  1:250002/1‑149 (MQ=255)                                                                                                                              tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGac                     <  1:67073/149‑1 (MQ=255)                                                                                                                               aagaCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGaca                    <  2:521541/146‑1 (MQ=255)                                                                                                                                  aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGc                 <  1:61784/149‑1 (MQ=255)                                                                                                                                    cGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGGTTTGCTCAAATTCGCGCTGCGACCCCCTa               >  2:557488/1‑149 (MQ=255)                                                                                                                                     gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCGGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTGCGAGCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTaa              <  1:1001069/149‑1 (MQ=255)                                                                                                                                      aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAAt             <  1:165230/149‑1 (MQ=255)                                                                                                                                        aGCTGCAATTTGATCGGAACTGGTATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGa           <  1:177798/149‑1 (MQ=255)                                                                                                                                          ctgtcagttGACCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGata         <  1:622455/142‑1 (MQ=255)                                                                                                                                           tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat        <  1:1095221/149‑1 (MQ=255)                                                                                                                                           tGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGatat        >  1:214899/1‑149 (MQ=255)                                                                                                                                            gCAATTTGATCGGAATTCATATCGGCTCTCGCCAGCACTTCTCCCAGCGCGGAGCGTCGGGGGGCCCACATTGCCATTGGGTCGTGTTCTACCCACCCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGCTAtt       <  2:698493/149‑1 (MQ=255)                                                                                                                                            gCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAtt       <  1:1333726/149‑1 (MQ=255)                                                                                                                                                 ttGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAAGGATATTGGCAt  >  1:841730/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCAT  >  NC_000913/4116856‑4117148

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