Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,117,030	G→T	L65M (CTG→ATG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,030	0	G	T	96.2%	158.0 / NA	54	L65M (CTG→ATG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (24/27);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAG‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCA  >  NC_000913/4116737‑4117310                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                          cAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTACCGGTTTCTTTTTTCCCAGACAATCTTGGTTTCACCCTGGTTCGTATTACCGATAGCTGCAATTTGATCGGACCTGTTATCGGCTTTCGCCAT‑CACttt                                                                                                                                                                                                                                                                                    >  1:419430/1‑300 (MQ=255)          accacGACGTGAACGCGCGCGAGAGCATTCCAAATGGGCAGGAATCCACTTAACTTTGGTGCCAGAAAAGTACGGGGCAATCACCGGACCGGTATTGCTGCGGTTATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTg                                                                                                                                                                                                                                                                          <  2:46069/301‑1 (MQ=255)            cacGACGTGCACGCCCGCGAGAGCCTTCCACATGGTCGATGATCTACTCCACTTTGGTGCCAGAAAAGTCCGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCCGATTTCTGCGGTACGACGGCACTTCCAGACAAGTGCGTTATAGATAGGCTTGCCCGTCTC‑TTTTTCCCAGTCCATAGTTGTTTCACGCTGGTTTGTAATTACTATAGCTGACCTTTGAATCGAATTGATATCGGCTTTCGCCCGCCACTTCTCCCacctggt                                                                                                                                                                                                                                                                         >  2:251002/1‑295 (MQ=255)                   ggcccgcgcgAGAGCGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTg                                                                                                                                                                                                                                                                  <  2:409248/296‑1 (MQ=255)                     cccgttCGAGACAGCCTTCTACAAAGTAGAGGAACCACTTCACTTTGGTGCCAGAAAAGTGCGTGTCAATCACCAGACCGGTAATGCTGCGGATATAATCTTCTAAACCGGCACGTTTTAAATGCTCGCAGATTTATGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTATCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGt                                                                                                                                                                                                                                                               <  2:195038/296‑1 (MQ=255)                     cACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATATTATAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACCATGGCGTTCTAGATAGGCTTGCCGGGTTC‑TTTTTCCCTGACCATAGTGGTTTCACGCTGGTTCGTAATACTGATAGCTGCAATTTGTTCGGAACTGATACAGACTTTCGCCCA‑CACTTATACCAGCGTGGACCGTTGGGt                                                                                                                                                                                                                                                               >  2:301565/1‑301 (MQ=255)                             cgAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGACAGACAATGGCGTTATAGAGAGGCTTGCCGGTTTC‑TTTTTCCCAGACAACAGTAGTTTACAGCTGGTTCGTATTACCGATAGCTGCACTTCGTTCGGACATGATATCGGCTGTCGCAAT‑CTCTTCTCCACGAGCGGCGATTGGGTTGGCCCAg                                                                                                                                                                                                                                                        >  2:477564/1‑300 (MQ=255)                                 agCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGAtt                                                                                                                                                                                                                                                     <  1:426577/299‑1 (MQ=255)                                        cccaaggtctgggATTCACTTTCCCTTTGTGCCAGGAAAGTACTTGTCCAAGACCCGACCGGTAATGCTGCCGATATAATTTTCTAAACCGTCACGTTTTTAAAGGTCGCAAATTTTTGCGGTACGACCGCACTGCCAGACAATGGCGTTATAGATAGGCTTACCCGTTTCTTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTgg                                                                                                                                                                                                                                             <  2:419420/290‑1 (MQ=255)                                                    gATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATAGAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTa                                                                                                                                                                                                                                  <  2:206361/299‑1 (MQ=255)                                                         acttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAc                                                                                                                                                                                                                            >  2:312861/1‑300 (MQ=255)                                                             cacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGCGCTTTGGTTTGCCTAGATTTCCATTGGGTCGTTTTCTACCCAACCTgtt                                                                                                                                                                                                                       >  1:1977/1‑299 (MQ=255)                                                               cttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTTGGTGGCCCAGATTTTCATTTGGTCGTGTTCTACCaaaaacacggg                                                                                                                                                                                                                     >  2:281480/1‑293 (MQ=255)                                                                 ttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCTGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCACAGATTTCTGCGGTACGACGTCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGATTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTg                                                                                                                                                                                                                   <  1:403540/301‑1 (MQ=255)                                                                   ggTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGTCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCACACAATAGTGGTTTCACGCTCGTTCGTTATACCGTTAGATTCACTTTGATCGGAACTGTTTTCGGCTTTCGCCAT‑CCCTTCTCCCTTCGTCGACCTTTGGCTTGCCCAGTTTTCCATTGGGTCTTGTTCTCCCCAACCTCGTTTTagg                                                                                                                                                                                                                 >  2:311528/1‑301 (MQ=255)                                                                    gTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGTTTCGTGTTCTACCCAACCTGGTTTTGGGt                                                                                                                                                                                                                >  1:298010/1‑301 (MQ=255)                                                                      gCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTa                                                                                                                                                                                                               <  1:312866/300‑1 (MQ=255)                                                                                tACGGTTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCaaaaccaggttgggtagaacacgac                                                                                                                                                                                                      <  1:281484/299‑24 (MQ=255)                                                                                 aCGGGTTAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaaa                                                                                                                                                                                                   <  1:452206/301‑1 (MQ=255)                                                                                           cACCAGACCGGTATTGCTGCTGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGc                                                                                                                                                                                         <  1:282901/301‑1 (MQ=255)                                                                                              cAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGa                                                                                                                                                                                       >  1:345488/1‑300 (MQ=255)                                                                                                  ccGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGCCAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGCCAATAGTGGTTTCACGCCGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CCCTTCTACTAGCGTGGAGCTTTGGTTGGCCCTGATTTCCATTGGGCCGTGTTCTACCCAAACTGGTTTGAGGTAGTTTTCCTCAAAttgccccgtcaaccac                                                                                                                                                                                   >  2:244635/1‑286 (MQ=255)                                                                                                           ctgcGGATATAATCTTCTAAACCGTCACGTTTTAAATGCCCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGTACTGATATCAGCTTTCGCCAT‑TACTTCTACCAGCGTGGAGCTTTGGCTGGCCGTGATTTTCGTTGGGTCGTTGTTTTCCCATCTTGGTTGTGGGTAGTGTTGCTTTAATTTAGGCTTGGACCCGCTAATGatta                                                                                                                                                                         >  2:198118/1‑299 (MQ=255)                                                                                                                    ttccgttatAAGCCGTAACTGTTTAAATTCTCGCAGATTATTGCGGTACGACGGCACTGCCAGAAAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGt                                                                                                                                                                <  2:154653/293‑1 (MQ=255)                                                                                                                      aTCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCgggg                                                                                                                                                              >  1:499512/1‑298 (MQ=255)                                                                                                                         ttctaatccgtacagagttACATGCTCGCAGATTTCGGCGGTACGAAGGCACGGCCAGACAATGTCGTTATAGATAGGCTTGCGGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACACTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATcc                                                                                                                                                           <  2:461180/284‑1 (MQ=255)                                                                                                                             aaaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGATAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGTCCAGATTTTCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCCACCCGATAATGATATTGGCATCGTGATCCAttt                                                                                                                                                       >  1:224348/1‑300 (MQ=255)                                                                                                                             aaaCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGTTGGCCCAGATTTCCATTGGGTCGTGTTTTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATTTTGGCATCGtagatcggaa                                                                                                                                                       >  1:154655/1‑292 (MQ=255)                                                                                                                                ccGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCACATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGa                                                                                                                                                    >  1:99904/1‑301 (MQ=255)                                                                                                                                           tAAATGCTCGCAGATTTCAGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                                                                                         <  1:301570/301‑1 (MQ=255)                                                                                                                                             aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGCCACGCTATTGATATTGGCATCGTGATCCCTTACGACCGCGCgggcg                                                                                                                                         >  1:213389/1‑297 (MQ=255)                                                                                                                                               tGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGATTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATGGTGATCCATTATGACCGCGCGGGAGCtgg                                                                                                                                     >  2:255912/1‑301 (MQ=255)                                                                                                                                                        tttttTGCGGTACGACGGCACTTCCCGACAATGGCGTTATAGATAGGCTTTCCGGTTTC‑TTTTTCCCAGACAATAGTTGTTTCACGCTGGTCCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCt                                                                                                                            <  2:334405/300‑1 (MQ=255)                                                                                                                                                                         actacaccAAACAATTGCTTTATAAATAGAATTGCCGGTTTC‑TTTTTCCTGAACGTTATGTGTTTCACGCTTGTTCGTAAAACCGATAGCTGCAATTTGATCGTAACTGATATCTGCTTTCGCAAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGata                                                                                                           <  2:53079/295‑1 (MQ=255)                                                                                                                                                                          cACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGatat                                                                                                          <  1:34948/301‑1 (MQ=255)                                                                                                                                                                              gCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTTGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGGTTTGGGTAGATTTGCTCAAATTCGCGGTGCGACAAGCTAATGATATTGGCATCGTGATCAATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGataagata                                                                                                      >  1:53080/1‑296 (MQ=255)                                                                                                                                                                                      aTGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAg                                                                                              <  1:141289/301‑1 (MQ=255)                                                                                                                                                                                                agaaatgattGCCGGTTTC‑TTTTTCCCAGTCACTAATGGTTTCTCGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTAATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTTGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTGTTGGGTAGATTTGCTCAAATTCGCTCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTTCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGttt                                                                                      <  1:477574/291‑1 (MQ=255)                                                                                                                                                                                                 gatagGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTTACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTATTCAACCCGGTTTTGGGTAGATTTGCTCAAATTCGCCCTGCGACACGCTATTGATATTGGCATCGTGATCCATTACGATCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAAt                                                                                   >  1:44430/1‑301 (MQ=255)                                                                                                                                                                                                       cTTGTCGGTTTC‑TTTTTCCCTGACAATAGTTGTGTCACCGTGGTTCGTAATACCGGTAGGTGCAATTTGGTCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTcc                                                                              <  2:345481/300‑1 (MQ=255)                                                                                                                                                                                                         tGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTCATACCGATAGCTGCAATTTGATCGGAACTGATATTGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCACTGGGTCGTGTTCTACACAACCTGGTTTTGGGTAGATTTGCTCACACTCGCGCTGCGCCACGCTTATGATATTGGCATCGTGATCCATCATGTCCGCGCGGGATCTGGTGGTCACCTGGTCGAGGGCAAGGATATTGTTTTTTCTAGTCTGTTTTATTTTTCCCgg                                                                           >  2:302341/1‑300 (MQ=255)                                                                                                                                                                                                                         cccAGACAATAGT‑GTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAAACTGGTTTTGGGTAGATTTGCCCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTTCCACCCCGCGGGAGCTGGTGGTGCCCTGGTCGATCCCAACGATACTTTTTTTTTCAGGCCTGTTTCATTTTCTCGTTGTCATATTaaacag                                                            >  1:177596/1‑295 (MQ=255)                                                                                                                                                                                                                               cAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGCGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGTGCTGCGCCCCGCCAATGATATTGGCATCGTGAACCATTACGACAGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCTCAACGATATATTTTTTTTCAGTCATTTTTTATTTTCCCGTTGTCCTATTCCAGctacgttt                                                      >  2:281271/1‑294 (MQ=255)                                                                                                                                                                                                                                  tAGTGGTTTTGCGCGAGGTCGAAATACATAGAGCTGCCAGTTGATGGGAACGGAGATCGGCGTTCGCCAG‑CTCTTCGACCAGCGTGGGGCATTGGGTGGCCCATATTTCCATTGGGTCGTGTTCTACCCATCCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACAAGCTAAGGATATTGGCATCGTGATCCATGACGACCGCGCGGGGGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCtttt                                                    <  2:1977/300‑1 (MQ=255)                                                                                                                                                                                                                                          tCACGCTGGTTCGTAATACCGAGAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGaa                                           <  1:255915/301‑1 (MQ=255)                                                                                                                                                                                                                                           cACGCGGGTTCGTAATACCGATAGCTGCAATTTGAGCGGAACTGATATCGGCTTGCGCCAT‑CACTTCTACCGTCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAg                                          <  1:88629/301‑1 (MQ=255)                                                                                                                                                                                                                                               cTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGt                                      <  1:394462/301‑1 (MQ=255)                                                                                                                                                                                                                                                 ggTTCGTAATACCGATAGCTGCAATTTGATCGCAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGTGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTTGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCTCTGGTCGAGCGCAACGATATATTTTTTTTCATTCATGTTTAATCGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTtg                                    >  1:457924/1‑301 (MQ=255)                                                                                                                                                                                                                                                  gTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTtgt                                   >  1:130669/1‑301 (MQ=255)                                                                                                                                                                                                                                                          ttacGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGATTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCt                           <  1:311533/298‑1 (MQ=255)                                                                                                                                                                                                                                                               aTAGCTGCAATTTGATTGGAACAGATATCGGCTGTGGCCAT‑CACTTCTACCAGCGTGGGGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGGTATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTttct                      <  1:251005/301‑1 (MQ=255)                                                                                                                                                                                                                                                               aTAGCTGCAATGTGAGCGGAATGGATGTCGGCTTTCGCCAT‑GACTGCGACCAGCTTGGCGTTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCACACTGGGTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTttct                      <  2:298005/301‑1 (MQ=255)                                                                                                                                                                                                                                                                            gATCGGAACTGATATCGGCTTTCGCCAT‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGa         >  1:237422/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                   aCTGATATCGGCGTTCGCCAT‑CAATTATAACAGCGTTGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATAGTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCa  <  1:198120/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                          CAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTC‑TTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAG‑CACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCA  >  NC_000913/4116737‑4117310

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑