Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,593	A→G	A419A (GCA→GCG)	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	1,593	0	A	G	100.0%	91.7 / NA	39	A419A (GCA→GCG)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/26);  total (13/26)

TCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAAC  >  NC_000913/1354‑1819                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                    tCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCCGGCCATTATCTCGGGGGTAGGTGATGGTATGCGCACCTTGCGTGGGGTCTCGGCGGAATTCTTTGCCGCGGTGGCCCccgc                                                                                                                                                                                                                          >  2:216942/1‑251 (MQ=255)   aCGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGGGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGcc                                                                                                                                                                                                                         >  1:183362/1‑250 (MQ=255)     gcgcCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGGATCAGTTTCTGCGTTCCGCAAAGCGACGGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGTCTTTCTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTGTGCCGCGCTGGCCCGCGCCAat                                                                                                                                                                                                                      <  2:143534/251‑1 (MQ=255)           gTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGGATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAAc                                                                                                                                                                                                                <  2:65751/251‑1 (MQ=255)             aTTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑t                                                                                                                                                                                                             >  1:99441/1‑251 (MQ=255)               ttCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑tt                                                                                                                                                                                                            <  1:216942/250‑1 (MQ=255)                tCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTCCTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGGGGTAGGCGTTGGTATGCGCACCTTGCGTGGGTTCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGt                                                                                                                                                                                                          >  2:18652/1‑251 (MQ=255)                        gCTGATTACGCAATCATCTTGCGCATACAGGATCAGTTTCTGCGTTCCGCAAAGCGACTGGGTGCGAGCTGAACGGGCAATGCAGGAATAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCGTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCAGGCCCGCGCCAATATCCACA‑TTGTCGCCATTg                                                                                                                                                                                                  <  2:173914/251‑1 (MQ=255)                               aCGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGACGCTGGCGGTGAAGGAACGGCCGGCCATTATCTCGGTGGTAGGAGATGGTATTAGCACCTTGCGGGGGTGCTCGGAGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCCACA‑TTATCTCCATTGCTCAggg                                                                                                                                                                                           >  1:110444/1‑251 (MQ=255)                                      catcTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGCAACGGCTGGCCCTTATCCCGGTGGTAGGTGCTGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCACTATCAACA‑TTTTCGCCCTTTCTCAGGGAtcttct                                                                                                                                                                                    >  1:139935/1‑251 (MQ=255)                                           tcgcaatacggtaTCAGCTTCTGCGTTCCGCAAGGCGACTGGGTGCGAGCTGAACGGGCAATTCAGGAATAGTTCTCCCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACCGAACGTCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCACATATCAACA‑TTGTCGCCATTGCTCAGTGATCTTCTg                                                                                                                                                                                   <  2:226251/235‑1 (MQ=255)                                                     gtatccgTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGGCTGGCCATTATATAGGTGGAAGGGGATGGTATGCGAACCCTGCGTGGGGTTTAGGCGAAATTTTTTGGCGCGCTGGGACGCGCCAATATCA‑CACTTGTCGCCAATGCTCGGGGCTCTTCTGAACCCTCaaattat                                                                                                                                                                     >  2:264113/7‑246 (MQ=255)                                                         cAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGCAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATATCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCg                                                                                                                                                                 >  1:109980/1‑251 (MQ=255)                                                             ttCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTgg                                                                                                                                                              <  1:68385/250‑1 (MQ=255)                                                              tCTGCGTTCCGCAAAGCGACGGTGTGCGAGCTGAGCGGGTAATGCAGGAAGAGTGCTACCTGGAACTGAAAGAAGGCTTACTGGGGCCGCTGGCGTTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGAGCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCATTAGCAGCA‑TTGTCGCGATTGCGCAGGGATCTTCTGAACGCTCAATCTTTGTCGTGGt                                                                                                                                                             <  2:181184/250‑1 (MQ=255)                                                                               gactttggtcgtgCTGAACGGGCATTGCAGGAAGAGTGATTGCTGGAAGTGAACGAAGGCTTTCTGGAGCCGATGGAGGTGACGGACCGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCTCCTTGCTTGGGAGCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATACCGATGATGCGAc                                                                                                                                           <  2:9183/239‑1 (MQ=255)                                                                                  tgggggCGAGCTGAACCGGCAATGCGGTATGAGGTCTACCTGGCACTGAGAGAAGGCTTACGGGGGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGGCGGCATTGCTCAGGGATCTTCTGAACTCTCAATCCCTGTCGTGGTAAATAACGATGATGCGACCa                                                                                                                                         <  2:128877/245‑1 (MQ=255)                                                                                  gtggtgcgagctaatctggtactGCAGGAGGAGTGCTACCTTGAACTGAAAGAAGGCTTACTTGGTCCGCTGGCGGTGACGGAACGGCTCGCCATTATCTCGGTGGTAGGTGATGGTATGCGCCCCTTGCGTGGCAGCTCGGCGAAAGTCTTTGCCGCGCTGGCCCTCGCCAATATCAACA‑TTGTCGCCATCGCTCAGGGATCTGCTGAACGCTCAATCTCTGTCGCGGTAAATAACGATGATGCGACCa                                                                                                                                         <  2:107179/228‑1 (MQ=255)                                                                                       gCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCGGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCAGTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCg                                                                                                                                   <  1:194185/251‑1 (MQ=255)                                                                                        cGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGATCTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGt                                                                                                                                  <  2:89263/251‑1 (MQ=255)                                                                                                                    tCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCaa                                                                                                       <  1:203279/250‑1 (MQ=255)                                                                                                                       aCCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAAt                                                                                                      <  1:187916/248‑1 (MQ=255)                                                                                                                            gAACTGAAAGAAGGCGTACTGGAGCCGCTGGCGGTGACGTAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGGTACTCAGCAGATGCTGTTCAATACCGATCa                                                                                              <  2:27863/251‑1 (MQ=255)                                                                                                                             aaCTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCGTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGAGGCTGTGCAATTCCGATCAg                                                                                             <  2:52144/251‑1 (MQ=255)                                                                                                                                     gaagGCTTACTGGAGCCGCTGGCGGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGCATGCGCCCCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATTTCAACA‑TTGTCGCCATTGCTCAGGGCTCTTCTGAACGCTCAATCTCTGTCGTGGTCAATAACGGTGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCgc                                                                                     >  1:133814/1‑250 (MQ=255)                                                                                                                                                 gagtcgcttgcggTGAAGGAACGGCTGGCCATTATCTCGGTGGTACGTGATGGTATGCGCCCCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGTCCCTCGACAATATCAACA‑TCGTCGCCATTGCTCTGGGATCTTCTGAACTCTCAAGCTCTGTCGTGGTCAATCACGATGAAGCGCCCACTGGCGTGCGCGTTTCTCATCATATGCTGTTCAATACCGATCAGGGTATCGAAgtgtgtgagat                                                                         <  2:175530/239‑8 (MQ=255)                                                                                                                                                     cgcgggcggTGACGGAACGGCTGGCCATTATCTCGGTTGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTTTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGc                                                                     <  2:183362/243‑1 (MQ=255)                                                                                                                                                         ggcggTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGt                                                                   <  2:99441/245‑1 (MQ=255)                                                                                                                                                                    aaCGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGc                                                      <  1:35778/251‑1 (MQ=255)                                                                                                                                                                       ggCTGGCCATTATCACGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGAGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGACATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATCCCGCTCAGGTGATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCggg                                                   >  2:8281/1‑250 (MQ=255)                                                                                                                                                                                  atcgcGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGa                                        <  1:237292/247‑1 (MQ=255)                                                                                                                                                                                            tAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCCACTGa                                >  2:260880/1‑249 (MQ=255)                                                                                                                                                                                              ggTGATGATATGCGCACCTTGCGAGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGGTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGAGCTGCTAGAGCAACTGAcgcg                            >  2:142238/1‑251 (MQ=255)                                                                                                                                                                                                    ggTATGCGCCCCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCGGCTGGAGCAACTGAAGCGTCaac                       <  1:72732/250‑3 (MQ=255)                                                                                                                                                                                                             aCCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGAGCGCGTTACTCATCAGATGCTGATCAATACCGATCAGGGTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGGTGCTGGAGCAACTGAAGCGTCCACAAAgatggcg             >  2:1273/1‑250 (MQ=255)                                                                                                                                                                                                                 tGCGTGGGATCTCGGCGAAATGCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCTTCAACAAAGCTGGCTgaa          <  1:103020/250‑1 (MQ=255)                                                                                                                                                                                                                 tGCGTGGGATCTCGGCGAAATGCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTgaa          <  2:135329/250‑1 (MQ=255)                                                                                                                                                                                                                     ggggATCTCGGCGAAATTCTTTGCCGCGCTGGCACGCGCCAATATCAACA‑TTGGCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGGCCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTTGGTGGCGGTGGCGGAGCGCTGCTGGAGCAACTGAAGCCTCAACAAATCGGGCTGAAGAATa     <  2:175886/250‑1 (MQ=255)                                                                                                                                                                                                                             cGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCGGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGGGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAAc  <  2:133814/246‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                    TCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAAC  >  NC_000913/1354‑1819

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: ‑