Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,463,952	G→A	S253F (TCC→TTC)	treC ←	trehalose‑6‑phosphate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,463,952	0	G	A	100.0%	131.5 / NA	43	S253F (TCC→TTC)	treC	trehalose‑6‑phosphate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/20);  total (23/20)

CAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGGACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTGCGCGTACCGCTGGATTCTCCC  >  NC_000913/4463710‑4464192                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 cAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGGGATGAAAATTAAAGGTCATCGCCAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCCATGCTCAAGGCTGGTGGAGAACATTTcc                                                                                                                                                                                                                                          >  2:213255/1‑250 (MQ=255)         gaAACCATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTCCGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCCCAAAGCCAGGTTTTGCCCGCGTCCATTTTTCACCCCCGGGATAATCGACCCTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCACTGCTCAAGGCTGGTGGAGAACATTTCACCTACggg                                                                                                                                                                                                                                  >  2:109476/1‑250 (MQ=255)                    cgGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTCGCCCGCGTCCCTTTTTCACCACCGGGATAATCGACCTTCCGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGGCATTAACCcgc                                                                                                                                                                                                                       >  1:186350/1‑251 (MQ=255)                         gATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTGAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAATGTCATCGACAATTCACTGCCTTTCAGAGCTGCGTATCGCTTGCAGTGCTCAATGCGGTTGGAGAACATGTCACCTACGGTCATGTTCCCGCGTgg                                                                                                                                                                                                                   <  2:240558/250‑1 (MQ=255)                            cGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGTAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCTCCACCGGTATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGAAAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTATCCCGCGTGGAGTa                                                                                                                                                                                                               <  1:55869/251‑1 (MQ=255)                            cGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCTGCGTGGCGTa                                                                                                                                                                                                               <  1:11684/251‑1 (MQ=255)                            cGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGGCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGGGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGGCAAAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGt                                                                                                                                                                                                                >  1:201125/1‑250 (MQ=255)                             gTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAAAGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTTCTCGACAATTCACTGCCTGTCAGAGCCGCGTATAGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTa                                                                                                                                                                                                               >  1:64598/1‑250 (MQ=255)                               ggTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAACGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGGCATCGACAATTAACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAac                                                                                                                                                                                                            >  1:51481/1‑251 (MQ=255)                               ggTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAac                                                                                                                                                                                                            <  1:248823/251‑1 (MQ=255)                                    cacCAGAACAAGGCATTCCATGCTACGTTGTGCTGTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGGGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCATAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACAGTCATTAACCCGCGTGGTGTaacctcatc                                                                                                                                                                                                       <  2:198357/251‑8 (MQ=255)                                                   ttCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCAGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTg                                                                                                                                                                                        >  1:95200/1‑251 (MQ=255)                                                            aCGTTGTGCAGTCCTTGTTGCCAGTGGCGGAACAATGTTTGCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTTATGAAAATTAAAGGTCATCGACAATTCACGGCCGGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTTGGGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGGAAACACATCGCGGTTCATCTCGTGCAAAAACt                                                                                                                                                                                <  2:133115/250‑1 (MQ=255)                                                             cGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCg                                                                                                                                                                              >  1:209950/1‑251 (MQ=255)                                                                 gtgCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGtgtg                                                                                                                                                                          >  1:147939/1‑251 (MQ=255)                                                                                     ggCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCATCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGTCCCGTCTGTGTa                                                                                                                                                       <  2:24735/250‑1 (MQ=255)                                                                                                gTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGTGATAATCGGCCTTCAGGTGATGAAAATTAAAGGTCATCTACAAGTCACTGCCTGTCGGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAGCATTTCACCTACGTTCATTAACCCGCGTGGCGTAAGCACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGcc                                                                                                                                            <  1:37259/250‑1 (MQ=255)                                                                                                gTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGcc                                                                                                                                            <  2:123697/250‑1 (MQ=255)                                                                                                           cacaAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCACGTCGCCGTCCa                                                                                                                                <  2:147939/251‑1 (MQ=255)                                                                                                            acaAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTTCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAACGACGCCCGTCGCCGTCCAg                                                                                                                               <  2:241896/251‑1 (MQ=255)                                                                                                             ccaagacAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGCAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTCTCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTCAACACATCGCGGTTCATCTCGAGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGCCGCCCGTCGCCGTCCAgg                                                                                                                              >  2:23716/7‑251 (MQ=255)                                                                                                                         tAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGaa                                                                                                                   >  2:189259/1‑250 (MQ=255)                                                                                                                             cAGCGTCCATTTTTCACCACCGGGATAATCGACCTTAAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCCATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGGGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAAcgcg                                                                                                              >  2:54882/1‑251 (MQ=255)                                                                                                                              aGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCgg                                                                                                             >  2:28584/1‑251 (MQ=255)                                                                                                                              aGCGTCCATTTTTCAACACAGGGATAATCGACCTTCAGGTGATGAAAATTAACGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCgg                                                                                                             >  2:31314/1‑251 (MQ=255)                                                                                                                                      tttttCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCTGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTAACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGAGGTTCATATCGAGCAAAAACTAGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTAGCCGTCCAGGTCTTCAGGGGAACGCGGGTCTTTgg                                                                                                     >  1:114807/1‑251 (MQ=255)                                                                                                                                           caccacCTGGATAATCGACCGTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATc                                                                                                <  2:30253/251‑1 (MQ=255)                                                                                                                                             ccaccGGGATAATCGACCGTCAGGTGATGAAAATTAAAGTTCATCGTCAATTCGCTGCCTGTCAAAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAGCTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCATGTCTGCAGGGAAACGCGGGTCTTTGGAGATCAg                                                                                              <  2:192333/251‑1 (MQ=255)                                                                                                                                                ccGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACGCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGAtt                                                                                           <  1:54882/251‑1 (MQ=255)                                                                                                                                                                  aGGTGATGAAACTTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCa                                                                          >  1:209738/1‑250 (MQ=255)                                                                                                                                                                      gatgaAAATGAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTTGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACAGCGCGGTTCAGCTCGTGCAAAAACTCGTGTGCTCGTGGCCCTTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCGTCAGGGAAACGCGGGTGTTTGGAGATCAGATTCACCACATCCAGGCTCAACCCg                                                                     <  2:179481/251‑1 (MQ=255)                                                                                                                                                                                   ggTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGAt                                                        >  2:134398/1‑251 (MQ=255)                                                                                                                                                                                    gTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATc                                                       <  1:213255/251‑1 (MQ=255)                                                                                                                                                                                                tCACTGCCTGACAGAGCCGCGTATCGCTGGCACTGCTCAAGGCTGGCGGAGAACATTTCACCTACGGTCATTAACCAGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCACCCCGTCGACCCCACGATCGGCCCAGAACTc                                           >  2:120837/1‑251 (MQ=255)                                                                                                                                                                                                 cACTGCCTGTCAGAGCCGCGCATCGCTGGCAATGCGCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACCCATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTGGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGGCCCCACGATCGTCCCAGAACTca                                          <  2:43496/251‑1 (MQ=255)                                                                                                                                                                                                   cTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTcac                                                                                        >  1:74007/1‑203 (MQ=255)                                                                                                                                                                                                   cTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTcac                                                                                        <  2:74007/203‑1 (MQ=255)                                                                                                                                                                                                                    gTATCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTg                       >  2:127758/1‑251 (MQ=255)                                                                                                                                                                                                                      aTCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTgcg                     >  2:126967/1‑251 (MQ=255)                                                                                                                                                                                                                      aTCGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTgc                      <  1:235727/250‑1 (MQ=255)                                                                                                                                                                                                                        cGCTGGCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCCCGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTgctcg                   >  2:215466/1‑248 (MQ=255)                                                                                                                                                                                                                            ggCAATGCTCAAGGCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACACCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGAGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGACCCAGAACTCACAGACTTTTTTCAGCTCTGCGCGTAcc               >  1:214205/1‑251 (MQ=255)                                                                                                                                                                                                                                         gCTGGTGGAGAACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTGCGCGTACCGCTGGATTCTccc  <  2:9841/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 CAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGCGGAACAATGTTTTCAACGCCACAAAGTCAGGTTTAGCCAGCGTCCATTTTTCACCACCGGGATAATCGACCTTCAGGTGATGAAAATTAAAGGTCATCGACAATTCACTGCCTGTCAGAGCCGCGTATCGCTGGCAATGCTCAAGGCTGGTGGAGGACATTTCACCTACGGTCATTAACCCGCGTGGCGTAAACACATCGCGGTTCATCTCGTGCAAAAACTCGTGTGCTCGTGGCCCGTCGGTGTAGAAGCGACGCCCGTCGCCGTCCAGGTCTTCAGGGAAACGCGGGTCTTTGGAGATCAGATTCACCACATCCAGGCGCAACCCGTCGACCCCACGATCGGCCCAGAACTCACAGACTTTTTTCAGCTCTGCGCGTACCGCTGGATTCTCCC  >  NC_000913/4463710‑4464192

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