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	97.4%	119.0 / ‑6.0	39	L65M (CTG→ATG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (20/17);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGA  >  NC_000913/4116747‑4117323                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                      ccacGCCGTGCACGCTCGCGAGAGCCTTCCACATGGTCGCGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACATTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTgg                                                                                                                                                                                                                                                                                      <  2:315599/301‑1 (MQ=255)          gCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATACCGGCTTTCGCCATCACTTCTACCCGCGTGGAGCATTggg                                                                                                                                                                                                                                                                             >  1:90969/1‑301 (MQ=255)                        gCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTc                                                                                                                                                                                                                                                                >  2:186521/1‑300 (MQ=255)                           ttCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGATCTTTGGGTGGCCCAGATTTCCAtt                                                                                                                                                                                                                                                            >  2:97909/1‑301 (MQ=255)                            tCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTg                                                                                                                                                                                                                                                           <  2:90972/301‑1 (MQ=255)                                           caccACTTCTATTTGGGGCCAATAAAGTGCGAGTCAATATCCATACCGATGTTCCTGCGGATAGTATCTTCTACACCGTCACGTTTTAAATGCTCGCAGATGTCTGCGGTACGACGTCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTAccc                                                                                                                                                                                                                                            <  1:302786/299‑1 (MQ=255)                                                     cttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGCTCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGTTCGTGTTCTACCCAACCTGGttt                                                                                                                                                                                                                                  >  2:39446/1‑301 (MQ=255)                                                         ggTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTggg                                                                                                                                                                                                                              >  1:77994/1‑301 (MQ=255)                                                          gTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTggg                                                                                                                                                                                                                              <  2:77996/300‑1 (MQ=255)                                                                 aaaaGGACGGGTCAATCACCAGACCGTGATTGCTGCGGATATAATCTTCTAAACCGTCACGGTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCGACCTGGTTTTGGGTAGAt                                                                                                                                                                                                                         <  2:305542/298‑1 (MQ=255)                                                                  aaaGTACGGGTCAATCCCCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCaaaaccaggttgggtagaacacg                                                                                                                                                                                                                     <  2:32361/301‑22 (MQ=255)                                                                             cAATCGCCAGACCGGTATGTCTGCGGATTTTACCTTCTAAACCGTCACGTTTTAAATGCTCTAAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTcgcg                                                                                                                                                                                                          <  1:253668/301‑1 (MQ=255)                                                                                        cttccgatcttgcGGATATAATCTTCTAATCCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGc                                                                                                                                                                                               <  2:225860/291‑1 (MQ=255)                                                                                                  tgcGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGTTTTGCTCAAATTCGCGCTGCGCCACGcagcacggaaa                                                                                                                                                                                     >  1:225852/1‑291 (MQ=255)                                                                                                                       cGTCACGTTTTAAATGCTCGCAGATTTCTGTGGTTTGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCAGCAATTTGATCGGTA‑CTGATATCAGTTTTCGCCATCACTTATACCAGCGTGTCGCTTTGGCTCGCTCAGATTTCCATTGTGTCGAGTTGCATCCTAATTGGTTTGGGGTAGCTTTGCCCAAATTCGGTTTGCGCCAGTCCAATAGTATTAGCATCGttaccacttccaa                                                                                                                                                                 >  1:299313/1‑288 (MQ=255)                                                                                                                                   aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                                                                                                      >  2:73435/1‑299 (MQ=255)                                                                                                                                       cTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATTGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGAGGCCCAGATTTCCATTGGGTCGTGTCCTACCCAACCTGGTTTTGGGTAGTTTTGCTTACATTCCCGTTTCGCCACGCTAATCAGATTTGCATCGTGATCCATTACGccccccggtgagtgggtg                                                                                                                                                >  1:124909/1‑286 (MQ=255)                                                                                                                                       cTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTTATCCATTACGACCGCGCGGGAGCtggtg                                                                                                                                                >  2:89403/1‑301 (MQ=255)                                                                                                                                        tCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTTCCGGTTTCTTTTTCCCAGACAATAGTGGTTTAACGCTGGTTCGTAATACCGATAGCTGCAATTTGAAAGGAA‑CTGATATCGGTTTTCGCCAACACTTCTACCAGCGTGGAGCTTTGGGTGGCCTAGATTTCCTTGGGCTCGTGTTCTACCTACACTGGTTTTGGGAAGATTTTCTCAAATTCGCGCTGTGACACGCTAATGATATGGTCATTGTGATCAAttcacaccggacaggcggtggtgg                                                                                                                                               >  1:204938/1‑279 (MQ=255)                                                                                                                                                  cTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCCATTTCGCGCTGCGACACGCTAATGATATTGGATTTGTGATCCATTACGACCGCGCGGGAACTGGTGGTGCCCTGGt                                                                                                                                      >  1:277127/1‑300 (MQ=255)                                                                                                                                                            ccgGCACTGCCAGGTAATGGCGTTATAGTTAGGCTTGCCGGTTTCTTTTTCCCAGAAAATAGTGGTTTCCGGCTGGTTCGTAATACCGATGGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACg                                                                                                                           <  1:186515/300‑1 (MQ=255)                                                                                                                                                                      cAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTTTTCTACCCAACCTGGTTTTGGGTAGCTTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATTGTGTTCCATTACGACCGACCGGGAGCTGGTGGTTCCGTGTTCTAGCTCAACGATATAttttt                                                                                                                 >  1:243192/1‑301 (MQ=255)                                                                                                                                                                                                        cTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATAtt                                                                               >  2:56826/1‑301 (MQ=255)                                                                                                                                                                                                           tttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACa                                                                            <  2:27159/301‑1 (MQ=255)                                                                                                                                                                                                            ttCCCAGACAATAGTGTTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAg                                                                           <  2:284009/301‑1 (MQ=255)                                                                                                                                                                                                              cccAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCTATTGGGCCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATCTTTAATTGTACCGTAGTCATATTAAAGc                                                                          >  2:119458/1‑300 (MQ=255)                                                                                                                                                                                                                   aCAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTACATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATTTTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGTTCGAGCGCAACGATATATTTTTTTTCAGTATTGTTTAATTTTTCCGTAGTATTTTTACAGCTAAg                                                                      >  1:205211/1‑299 (MQ=255)                                                                                                                                                                                                                                 aCGCTGGTTCGTAATACCGATAGCTTCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAgg                                                      <  1:73433/301‑1 (MQ=255)                                                                                                                                                                                                                                            tAATACCTATAGCTGCAATTTGTTCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGtt                                            <  1:56824/300‑1 (MQ=255)                                                                                                                                                                                                                                            tAATACCGATAGCTGCAATTTGATCGGAA‑CGGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGtt                                            <  2:124912/300‑1 (MQ=255)                                                                                                                                                                                                                                                    gaagcgtgaaTTTGATCGGAATCTG‑TTGCTGCTTTCGCCATTACGTCTACCAGTGTGGAGCTTTGGGTTGCCCAGATTTCCATTTGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTttc                                    <  1:89400/292‑1 (MQ=255)                                                                                                                                                                                                                                                    aTAGCTGCAATATGATCGGAA‑CTGTTATCGGCTTTTCACATCACTACGTCCAGCTCGGAGCGTTGTGTGGCCCAGATTTCCATTGGGTCGTGTACTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATGCGCGCTGCAACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTAGTCGAGCGCAACGATATATTTTTTTTCAGGCATGTTTAATTTTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTttct                                   <  2:204944/301‑1 (MQ=255)                                                                                                                                                                                                                                                        cTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCTTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTTGTTTCCTTttcttc                                 >  1:271080/1‑299 (MQ=255)                                                                                                                                                                                                                                                        cTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCATTTTCTTCca                               >  2:231648/1‑301 (MQ=255)                                                                                                                                                                                                                                                         ttcAATTTTCTCGAAA‑CTTATATCGGCTTTCGCCATCACTTCTTCCAGCGTATATCTTTGATTGGCCCAGATTTCCATCGGGTCGTTTTCTACCCACCCTGGTTTTAGGTATATTTGCTCAAATTCGCACTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGATCGAGCACAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCcac                              <  2:204485/299‑1 (MQ=255)                                                                                                                                                                                                                                                                          cTGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCaa              >  2:223026/1‑301 (MQ=255)                                                                                                                                                                                                                                                                           tGATATCGGCTTTCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTTTTCTGAAGGAGTTGTGGTTTCCTTTTTTTCCACAACACAGATATCGCaat             >  1:314411/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                    gttgcgacAGCACTTCTAGAAGCGTGGAGCTTTTGGTGGCCCAGTGTTCCTTTGAGGCGTGATCTACCCCAGCTGGTTTTGGGTAGATTTACTCAAATTCGCGCTGCGACCCGCTAATGATATTGGCATCCTGTTCCATTACGACCGCGCGGGGGATGGTGATGCCGTGGTCGAGCGCTACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAACGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTGTTCTTCCACAACACAGATATCGCAAGGCAAATgg     <  2:186188/293‑1 (MQ=255)                                                                                                                                                                                                                                                                                      ttCGCCATCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTTAAGGAGTTTTGGTTTCCTTTTCTTCCACAACACCGATATCGAAAGGCAAATGGCGa  >  2:308637/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                      CCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAA‑CTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGA  >  NC_000913/4116747‑4117323

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