Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,117,050	T→G	44.9%	S58R (AGC→CGC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,050	0	T	G	44.9%	‑0.3 / 117.0	96	S58R (AGC→CGC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (28/24);  new base G (11/32);  total (40/56)
Fisher's exact test for biased strand distribution p-value = 6.63e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-02

CACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCG  >  NC_000913/4116905‑4117196                                                                                                                                                  |                                                                                                                                                    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)                                                               ttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCGCTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGGCGCTTTCTTCCCAAGCTGGTTTTGGGTAGATTTGCTCa‑a                                                                                   >  1:1322335/1‑149 (MQ=255)                                                                tCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCa‑aa                                                                                  <  1:1190843/149‑1 (MQ=255)                                                                tCACGCGGGGTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCa‑aa                                                                                  <  2:383906/149‑1 (MQ=255)                                                                         ttCGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTg                                                                         >  1:276368/1‑149 (MQ=255)                                                                         gtCGTAATACCGATAGCTGCAATTGGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTg                                                                         <  2:348834/148‑1 (MQ=255)                                                                           cGTAATACCGATAGCTGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCG‑AATTCGCGCTGCg                                                                       <  1:993533/149‑1 (MQ=255)                                                                             tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGATTTGGGTAGATTTTCTCA‑AATTCGCGATGCGac                                                                     >  2:841206/1‑149 (MQ=255)                                                                             tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTTCCGCTGCCac                                                                     >  1:250002/1‑149 (MQ=255)                                                                             tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGac                                                                     <  1:67073/149‑1 (MQ=255)                                                                              aagaCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGaca                                                                    <  2:521541/146‑1 (MQ=255)                                                                                 aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGc                                                                 <  1:61784/149‑1 (MQ=255)                                                                                   cGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGGTTTGCTCA‑AATTCGCGCTGCGACCCCCTa                                                               >  2:557488/1‑149 (MQ=255)                                                                                    gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCGGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTGCGAGCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTaa                                                              <  1:1001069/149‑1 (MQ=255)                                                                                     aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAAt                                                             <  1:165230/149‑1 (MQ=255)                                                                                       aGCTGCAATTTGATCGGAACTGGTATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGa                                                           <  1:177798/149‑1 (MQ=255)                                                                                         ctgtcagttGACCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGata                                                         <  1:622455/142‑1 (MQ=255)                                                                                          tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGatat                                                        <  1:1095221/149‑1 (MQ=255)                                                                                          tGCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGatat                                                        >  1:214899/1‑149 (MQ=255)                                                                                           gCAATTTGATCGGAATTCATATCGGCTCTCGCCAGCACTTCTCCCAGCGCGGAGCGTCGGGGGGCCCACATTGCCATTGGGTCGTGTTCTACCCACCCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGCTAtt                                                       <  2:698493/149‑1 (MQ=255)                                                                                           gCAATTTGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATAtt                                                       <  1:1333726/149‑1 (MQ=255)                                                                                                ttGAACGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAAGGATATTGGCAt                                                  >  1:841730/1‑149 (MQ=255)                                                                                                     cGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGACCAGATTTCCATTGGGTCGTGTTCTACCCCACATGCAGCGGCCGAGGACGGCGTA‑AATTCGCGCTGCGACACGCTAATGATCTTGGCATCGTGa                                             >  2:309264/1‑149 (MQ=255)                                                                                                        accTGATATCGGCTTTCGCCCGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATcc                                          <  2:886355/147‑1 (MQ=255)                                                                                                        aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATcc                                          >  1:330115/1‑149 (MQ=255)                                                                                                        aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCACCGTGATcc                                          >  2:318913/1‑149 (MQ=255)                                                                                                         aCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCa                                         <  2:1408276/149‑1 (MQ=255)                                                                                                          cTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCAt                                        >  2:681664/1‑149 (MQ=255)                                                                                                          cTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTATGCCAAACCTGGTTTTCGGT‑GCAATGATGAGAATATGCGCCGAGACCCGCTCACTATATTGGCATCGTGATCCAt                                        >  1:1053474/1‑149 (MQ=255)                                                                                                           tGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCAtt                                       >  1:602768/1‑149 (MQ=255)                                                                                                               atCTGCTTTCGCCAGCACTTCTACCCGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGTTAGATTTGCTCA‑AATTCGCGCTGTGACACGCTAATGATATTCCCAGCGCGATCCATTACGa                                   <  2:863452/149‑1 (MQ=255)                                                                                                               atCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTGCCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTCCgc                                   >  2:45142/1‑148 (MQ=255)                                                                                                               accGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTCCAATTGGGAAGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGa                                   <  1:341184/147‑1 (MQ=255)                                                                                                                   gCTTTCGCCAGCACTTCTACCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcgc                               >  2:1379454/1‑149 (MQ=255)                                                                                                                       tCGCCAGCACTTCTCCCAGCGTGGAGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCgg                           <  1:599624/149‑1 (MQ=255)                                                                                                                       tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCgg                           >  1:740322/1‑149 (MQ=255)                                                                                                                        ggCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCggg                          <  1:1084284/148‑1 (MQ=255)                                                                                                                             gCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTTCTCAACCTGGTTTTGGGTAGCTTTGCTCA‑AATTCGCGCGGCGACACGCTGATGATATTGGCATCGTGATCCATTACGACCGCGCGGGCGCtg                     >  1:422531/1‑149 (MQ=255)                                                                                                                             gCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtg                     <  1:1374530/149‑1 (MQ=255)                                                                                                                              cACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtgg                    <  2:214899/149‑1 (MQ=255)                                                                                                                                   cgaCCAGCGTGCAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGgtggtg               <  1:91645/147‑1 (MQ=255)                                                                                                                                    tACCAGCGTGGAGCTGTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGc              <  2:1046406/149‑1 (MQ=255)                                                                                                                                          cGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGGGCTGGTGGTGCCCTGGt        >  2:150551/1‑149 (MQ=255)                                                                                                                                             ggAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGa     >  2:1106766/1‑149 (MQ=255)                                                                                                                                              gAGCGTTGGGTGACCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAg    >  2:1260995/1‑149 (MQ=255)                                                                                                                                               aGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgc   <  2:645334/149‑1 (MQ=255)                                                                                                                                                gCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgcg  >  1:1408485/1‑149 (MQ=255)                                                                                                                                                gCTTTGGGTGGCCCAGATTTCCATTGGGGCGTGCTCTCCCCAACATGGCATCGGGTACATATGTTCA‑AATTCGCAGTGCTCCACGCTTTTAGTAATGTGATCCTTGTCCGTTCCAGCCGCGCGGGAGCTGGGGGTGCCCTGGTCGAgcg  >  2:214004/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCA‑AATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCG  >  NC_000913/4116905‑4117196

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