Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,115,968	G→A	L419L (CTG→TTG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,115,968	0	G	A	97.9%	159.8 / NA	48	L419L (CTG→TTG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (25/21);  minor base C (1/0);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCAC  >  NC_000913/4115723‑4116214                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       gtggtcttcCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCGGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGAGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGGCTGGAACTGCATCAAGAAAt                                                                                                                                                                                                                                                   >  2:79942/5‑251 (MQ=255)    ttcttcCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTg                                                                                                                                                                                                                                                 <  1:10788/250‑1 (MQ=255)    ttcttcCAACGCCATCTCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATCTCGGACTGGAACTGCATCAAGAAATTGt                                                                                                                                                                                                                                                >  2:314972/1‑251 (MQ=255)           cccgccttgccgcgTTTAACCGCTTTTTTCCAGCCTGCGTACCGGTAATTACGCTCAGTGGTTTTTATGCCTGGTCGGAATTCGCGCTTAATCCCCGCTTTCTCTTGCATCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGTGCGCTCAACGCGGGTGCCTAGAATATCGGACTGGAACTGCATCAATAAATTGTTTGCTAc                                                                                                                                                                                                                                         <  1:179930/241‑1 (MQ=255)             cGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTg                                                                                                                                                                                                                                       >  1:111352/1‑251 (MQ=255)              gCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCCAGGTTCTGCCAGAAGCCAACGGCCAGGCCTGCGAGATAGGCCGCACCCAACGCGGTGACTTCACGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCCAGACATTGTTTGCTACTgc                                                                                                                                                                                                                                      >  2:39822/1‑251 (MQ=255)                   cgcgcgTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGAGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGATCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGccac                                                                                                                                                                                                                                 >  1:401536/1‑251 (MQ=255)                                     tttCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCATCATCCACGCGCAGGGCg                                                                                                                                                                                                                <  2:136315/250‑1 (MQ=255)                                           gcctgcgaacgggaaaTTACGCTCAAAGGATTCGATGCACGGACGGCACCAGCGCCCACTCCCCGCTTTCTCTTGTAGCTCGACGAGGTTCTGCCAGAAGCCAAAAGACAGGCCTGCGAGATAGGCCGCACACAATGCGGTGACTCCGCGCACTTCCCGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTACATCCGGAAATTGTTTGCTACCGCGCCACCATCCACGCGACGGGCGAGCAGAc                                                                                                                                                                                                         >  2:53818/15‑251 (MQ=255)                                              tGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACgg                                                                                                                                                                                                       >  2:390207/1‑250 (MQ=255)                                                      gttAATTCCGCTTAGTGGTTTAGGTGCTTCGACGCAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCTTCGAGAATCTACCATGAGCCCACCGCCAGGCCTGCGAGCTAGGTCGCACCCAACGCGGTGCCTTCCCCCACTTTCGCGCGCCTAACGCCCGCCCCGAGAATATCGGACTGGCACTGCGTCAAGAAATTGTTTGCTACTGCGCCACCCTCCACGGGTAGGGCGtttattcggatccccgcg                                                                                                                                                                                              >  2:185125/3‑234 (MQ=255)                                                          attccGCTCAGTGTTTTCGATGCCTGGACGGAACTCGCGCTCAATCCCCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCAAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCgg                                                                                                                                                                                          <  2:396512/247‑1 (MQ=255)                                                          aTTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCgg                                                                                                                                                                                          >  1:372359/1‑251 (MQ=255)                                                           ttACGCTCTGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCgg                                                                                                                                                                                          <  2:81096/250‑1 (MQ=255)                                                                       gTTTCGATGCCTGGAGGGAACTCGCGCTCAACCACCGCTTTCTCTCGCACCTCGTCGAGGTTCTGCCAGAAGCAAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTc                                                                                                                                                                             <  2:161703/251‑1 (MQ=255)                                                                            gATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGc                                                                                                                                                                         >  1:362499/1‑250 (MQ=255)                                                                               gCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTATGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGAGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGt                                                                                                                                                                     >  1:13327/1‑251 (MQ=255)                                                                                     aCGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCAcgcg                                                                                                                                                               >  1:59474/1‑251 (MQ=255)                                                                                                  ccAATCACCGCTTCCTCTTGCGGCGCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAAt                                                                                                                                                  <  2:13327/250‑1 (MQ=255)                                                                                                       cACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACt                                                                                                                                             >  2:248414/1‑251 (MQ=255)                                                                                                        aCCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTc                                                                                                                                            >  1:251955/1‑251 (MQ=255)                                                                                                          cGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCAcg                                                                                                                                                                 <  1:91036/228‑1 (MQ=255)                                                                                                          cGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCAcg                                                                                                                                                                 >  2:91036/1‑228 (MQ=255)                                                                                                              ttCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGc                                                                                                                                        <  2:264066/249‑1 (MQ=255)                                                                                                                 gctTGCACCTCGTGGGGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATGGTTGGCTACTGCGCCACCATCCGCGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCGGCATCGTTGCCAGCTCGTCACGCGTCTGATAAGCAATAGACTCCAGCGTcg                                                                                                                                    <  2:41984/249‑1 (MQ=255)                                                                                                                      cAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGt                                                                                                                               >  1:273566/1‑250 (MQ=255)                                                                                                                            gtcgAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAAtgt                                                                                                                        >  1:65382/1‑251 (MQ=255)                                                                                                                               gaggtttggCCAGAAGCCACCCGCCAGGCCTGCGAGAGAGGCCTCACCCAATGCGGTGACTTCGAGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAAGATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCTCGCGCATGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGTTTCCAGCATGTCACGCGTCTGATAAGCAATAGCCTCCAGAGTCGCGCGTATCAtgtgtt                                                                                                                     <  2:304131/243‑3 (MQ=255)                                                                                                                                            aaGCCACCCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGATGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCGTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCAcccc                                                                                                        >  1:118395/1‑251 (MQ=255)                                                                                                                                                       cAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGAc                                                                                             <  1:378523/251‑1 (MQ=255)                                                                                                                                                                 agaTAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGCGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCCGACCGACAATCGc                                                                                   >  2:337178/1‑251 (MQ=255)                                                                                                                                                                     aGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCcg                                                                               <  2:239492/251‑1 (MQ=255)                                                                                                                                                                       gCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCcgcg                                                                             <  1:126289/251‑1 (MQ=255)                                                                                                                                                                        ccGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCcgcgc                                                                            <  1:102745/251‑1 (MQ=255)                                                                                                                                                                              ccAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCCGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATAc                                                                      >  1:336686/1‑251 (MQ=255)                                                                                                                                                                                             cgcgCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCATGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGgc                                                       >  1:170669/1‑251 (MQ=255)                                                                                                                                                                                                          gcgcTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCgg                                          <  1:40423/251‑1 (MQ=255)                                                                                                                                                                                                          gcgcTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTGTGCTGCTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCgg                                          <  2:169483/251‑1 (MQ=255)                                                                                                                                                                                                                   gcgGGTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCg                                 <  1:248414/251‑1 (MQ=255)                                                                                                                                                                                                                      ggTGCCGAGAATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACAGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACTGAAAATCGCCCCGCGCGCATACGGGCCCCAGAACGGCGCACCCAGCCCGGTCAATGCCggc                               >  1:169963/1‑249 (MQ=255)                                                                                                                                                                                                                            gagaATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACAt                        >  1:289628/1‑251 (MQ=255)                                                                                                                                                                                                                            gagaATATCACACAGGAACTGCATCACGCAATTGTTTGCTACTACGCCAACATCCACGCGCAGGGCGCGCAGACGGATACCAGCGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATACGCCATAGACTCCAGCGTCGCGCGTATCATGTGGTTAGCGTTCACCCCACGAGTCAGACGGAACATCGCCCCGCGAGCATACGGGTCCCAGTACGGCGCACCCAGCCCCGCAAATCCCGGTACCACAt                        >  2:56475/1‑251 (MQ=255)                                                                                                                                                                                                                             agaATATCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCCCCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTTCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCACAGTACGGCGCCCCCAGCCCGGTAAATGCCGGAACCACATa                       <  1:117447/251‑1 (MQ=255)                                                                                                                                                                                                                                atatCGGACTGGAACTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATacac                    <  2:303389/251‑1 (MQ=255)                                                                                                                                                                                                                                            aaCTGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGtt        >  2:139443/1‑251 (MQ=255)                                                                                                                                                                                                                                               tGCATCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTg     >  2:267416/1‑251 (MQ=255)                                                                                                                                                                                                                                                 cATCAAGAAATTGTTTGCGACTGCGACACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATATACTCCAGCGGCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGTCCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACTGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCa   <  1:282473/251‑1 (MQ=255)                                                                                                                                                                                                                                                  aTCAAGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCAc  >  1:342144/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       GTGTTCTTCCCACGCCATCGCGCGTTTAACCGCTTTTTTCCAGCCTGCGTAACGGTAATTACGCTCAGTGGTTTCGATGCCTGGACGGAACTCGCGCTCAATCACCGCTTTCTCTTGCAGCTCGTCGAGGTTCTGCCAGAAGCCAACCGCCAGGCCTGCGAGATAGGCCGCACCCAATGCGGTGACTTCGCGCACTTCCGGGCGCTCAACGCGGGTGCCGAGAATATCGGACTGGAACTGCATCAGGAAATTGTTTGCTACTGCGCCACCATCCACGCGCAGGGCGTGCAGACGGATACCAGAGTCGGCCTGCATCGCTTCCAGCACGTCACGCGTCTGATAAGCAATAGACTCCAGCGTCGCGCGTATAATGTGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCAC  >  NC_000913/4115723‑4116214

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