Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	393	T→G	R19R (CGT→CGG)	thrA →	fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	393	0	T	G	100.0%	116.0 / NA	38	R19R (CGT→CGG)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (21/17);  total (21/17)

GGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGGTTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTG  >  NC_000913/145‑631                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               ggCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAAGTTTTCTGCGGGtt                                                                                                                                                                                                                                             >  1:296474/1‑251 (MQ=255)      aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGGTTTCTGCGGGTTGCCg                                                                                                                                                                                                                                         >  1:176071/1‑250 (MQ=255)           gcagacagaTAACAATTACAGAGTACACAACATCCATGAAACACATTAGCACCAACATTACCACCACAATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACCCAGAAAAAAGCCAGCACCTGACGGTGCGGGCTTTTTTT‑‑CGACCAAAGGTAACAAGGAACCAAGCATGCGTGAGTTGAAGTTCGGCGGTACATCAGTGGCAAATGAAGAACGTTTTCTGCGGGTTggctctattct                                                                                                                                                                                                                                  >  2:39779/2‑241 (MQ=255)             agaaaaaTAAAAATTACAGAACACACAACTACCAGGAAACACATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCAGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGACCGTTTTCTGCGGGTTGCCGATATTCTg                                                                                                                                                                                                                                 <  2:244405/245‑1 (MQ=255)                                  tACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg                                                                                                                                                                                                            >  2:217253/1‑251 (MQ=255)                                                                                tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAaccacc                                                                                                                                                              >  1:105873/1‑251 (MQ=255)                                                                                  accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCtg                                                                                                                                                            >  2:83674/1‑251 (MQ=255)                                                                                        aCCACAGGTACGGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCGGCCAGTGCGGGCTTTTTTTT‑CGACCAAAGTTAACGAGGTAACAACCATGCGAGTGTGGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCg                                                                                                                                                      <  2:446003/251‑1 (MQ=255)                                                                                        aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCACCCACCtggtgg                                                                                                                                                        >  1:125604/1‑249 (MQ=255)                                                                                              gatAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACAATGCGAGTGTTGAAGTTCGGCGGTTCATCAGTGGCAAATGCAGAACTTTTTCTGCGGGTTGCCGATATTCTGTAAAGCAATTCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTTCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGAtt                                                                                                                                                <  1:21230/249‑1 (MQ=255)                                                                                                  aCGGTGCGGGGTGACGCGCAAGGGAAACACAAAAAAAAGCCCGCACCTGACAGTGCGGTCTTTTTTTT‑CGACCAAAGTTAACGAGGTAACAACCATGCGAGTGTTGAAGTTGGGCGGTACATCAGTGGCAAATGCAGACCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCTTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGaaa                                                                                                                                            <  2:492670/251‑1 (MQ=255)                                                                                                              gaggcgTACAGGAAACACAGAAAAAAGCCCGTACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTATCGAGTTAACAACCTTGCGAGTGTTGAAGTTCGGCGGTACATCTTTGTCAAATGCAGAACTTTTTCTTCGGGTTTCCGATATTCTGTAAAGCAATGCCAGGAAGGGGCAGGTGGCCACCTTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGc                                                                                                                                 <  1:468293/247‑1 (MQ=255)                                                                                                                           aaaacaaaaaaaaGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACAAAAGGTAACGAGGTAACAACCATGCGAGTGTGGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCGGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGGGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCttt                                                                                                                   <  2:305598/244‑1 (MQ=255)                                                                                                                            accaaGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATCCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTa                                                                                                                  <  1:42340/247‑1 (MQ=255)                                                                                                                                          gcGCACCTGAAAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAAAAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGAAATGCCAGGCAGGGGCAGGGGGGCAGCGTAATCTCTGCACCCGCCAAAATCAACAAACACATGGTGGCGATGGTTGAAAAAAACATTAGCGGCAAGGGTGCTTTAAACAatatggat                                                                                                      >  2:221102/2‑245 (MQ=255)                                                                                                                                          ccGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAatat                                                                                                          <  1:221102/245‑1 (MQ=255)                                                                                                                                             ccccTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATGAGCGTCCAGGATGCTTTACCCAATAGCAGCGATGc                                                                                                 <  1:391447/249‑1 (MQ=255)                                                                                                                                             cACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCCATAACAGCGGTGc                                                                                                 >  2:215490/1‑251 (MQ=255)                                                                                                                                               ccTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCg                                                                                               <  1:83674/251‑1 (MQ=255)                                                                                                                                                                 ttttttCGACCAAAGGAAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGCTATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAACTCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAAGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGa                                                                               >  1:449023/1‑250 (MQ=255)                                                                                                                                                                  tttttCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGCTGCTTGAAAAAACCATTAGCGGCCCGGAATCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGACGAAc                                                                             >  2:68595/1‑251 (MQ=255)                                                                                                                                                                   ttttCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAAc                                                                             <  1:217253/250‑1 (MQ=255)                                                                                                                                                                         aCCAAAGGTAACGAGGTAACAAACATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAAACACCAACCACCTGGAGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCACTAACAGCGATGGCGAACGGATTTTTCCCGAACtttttg                                                                       >  2:337001/1‑248 (MQ=255)                                                                                                                                                                           cAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACAATTACGGGCCAGGAAGCTTTACCCAATATCAGCGATGCCAAACGGATTTTTGCCTAACtttgaacgg                                                                    >  1:272452/1‑245 (MQ=255)                                                                                                                                                                            aaaGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACggg                                                                   >  1:217260/1‑251 (MQ=255)                                                                                                                                                                                 tAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTcg                                                              >  1:268263/1‑251 (MQ=255)                                                                                                                                                                                          aacaacCATGCGAGATTTGAATTTCGGCGTTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGTACTTTTGACGGGACTCgccgccgcc                                                      <  1:270624/250‑1 (MQ=255)                                                                                                                                                                                           acaacCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGccc                                                     >  1:206802/1‑250 (MQ=255)                                                                                                                                                                                                 aTGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCggg                                              >  1:103520/1‑251 (MQ=255)                                                                                                                                                                                                     gatggtTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTc                                          <  1:367865/247‑1 (MQ=255)                                                                                                                                                                                                                   ggcggTACATCAGTGGCAAATGCAGACCTTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAAtt                            <  1:411495/251‑1 (MQ=255)                                                                                                                                                                                                                       gTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGaa                         >  1:464860/1‑250 (MQ=255)                                                                                                                                                                                                                           aTCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCCCCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATAAGCGATGCCGAACGTATTTTTGCCGAACTTTTGCCGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGaaaaatt                    >  2:428184/1‑248 (MQ=255)                                                                                                                                                                                                                             aaGTGGCAAATGCAGAACGTTTTCTGGGGGTTGCCGATTTTCTGGAAAGCAATGCAAGGCAGGGGCATGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGTCGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTtc                  <  1:293461/250‑1 (MQ=255)                                                                                                                                                                                                                                     aaTGCAGAACGTTTTCTGAGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGACCGTATTTTTGCCGACCTTTTGACGGGACGCGACGCCGACAAGCAGGGATTCCAGCTGGCGCAATTGAACACTTTCGTCGATc           >  1:222856/1‑250 (MQ=255)                                                                                                                                                                                                                                          aGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCAAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAACCCATTACCGTCCCGGATGCTTTACCCAAAATCAGCGATGCCGAACGTATTTTTGCCTCACTTTTGACGGGACACGCCCCCGCCCCGCCGGGCTTCCATCTGACGCACTTGAAAACTTTCGTCGATCAGGAAt     >  1:216798/1‑251 (MQ=255)                                                                                                                                                                                                                                             ccGTTTTCTGCGGGTTGCCGTTATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTg  <  1:337001/250‑1 (MQ=255)                                                                                                                                                                                                                                             aCGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGATTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTg  <  1:462414/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               GGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGAACTTTTGACGGGACTCGCCGCCGCCCAGCCGGGGTTCCCGCTGGCGCAATTGAAAACTTTCGTCGATCAGGAATTTG  >  NC_000913/145‑631

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