Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,211	A→T	A625A (GCA→GCT)	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	2,211	0	A	T	100.0%	183.3 / NA	59	A625A (GCA→GCT)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (26/33);  total (26/33)

CATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGGCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGC  >  NC_000913/1969‑2456                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      catcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGtgatg                                                                                                                                                                                                                                                <  1:136969/250‑1 (MQ=255)            cccGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGt                                                                                                                                                                                                                                    <  1:61987/251‑1 (MQ=255)                   aTTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCgg                                                                                                                                                                                                                             <  2:90775/251‑1 (MQ=255)                         gACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttct                                                                                                                                                                                                                       <  1:377939/251‑1 (MQ=255)                             gCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGAGTTCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctttca                                                                                                                                                                                                                   <  2:278264/251‑2 (MQ=255)                              cACTTCCAGCCAGGCAGTGGGGGATCAATATGCCGACTTCTTG‑CGCGAAGGGTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctttcag                                                                                                                                                                                                                  <  2:73613/251‑3 (MQ=255)                               aCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGCAAAATCGCGGCGTAAATTCCTCTATGACCACCA‑CGTTGGGGCTGGATTACCGGATATTGTGAACCTGCAAAATCTGCTAAATGCTGGTGATGAATTGATTTAGTTCTCCGCCAttctttcagg                                                                                                                                                                                                                 >  2:394727/1‑248 (MQ=255)                                  tCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGt                                                                                                                                                                                                                <  1:197270/249‑1 (MQ=255)                                     agccagGCAGTGGCGGATCAATATGCCGACTTCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATtgat                                                                                                                                                                                                                                         <  2:176907/221‑1 (MQ=255)                                     agccagGCAGTGGCGGATCAATATGCCGACTTCTTG‑CGCGAAGGTTTCCACGGTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATtgat                                                                                                                                                                                                                                         >  1:176907/1‑221 (MQ=255)                                       ccaagCAGTGGCGGATCAATATGCCGAATTCGTG‑CGCGAAGGTGTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCGCGTCGATGGATTACTGCCATCTGTTGCGGCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAGGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCCCAATGCTTGGGATGAATTGATGAAGTTCTCCGTCATTCTTTCTGGCTCGCttt                                                                                                                                                                                                         <  2:329950/247‑1 (MQ=255)                                              gTGGCGGTGCAATATGCCGTCTTCTTGACG‑GAAGGTTTCCCAGTTGTCACGCCTAACAAAAAGGCCAGCACCTCGGCGATGGAGTACTACCATCTTTTGCTTCATGCTGCGGAAAAGTCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGTGGCTGTATTACCGGTTATGTAGATCCTGCAAATTCTGCTCTATGCTTGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGCCTCGCTTTCTTatat                                                                                                                                                                                                  <  2:444020/251‑1 (MQ=255)                                                ggCGGATCAATAGTCCGACTTCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATGCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATc                                                                                                                                                                                                 <  2:422611/250‑1 (MQ=255)                                                 gCGGAACAATATGCCGACTTCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCtt                                                                                                                                                                                               <  2:326550/251‑1 (MQ=255)                                                         atatGCCGACTTCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGCATTGCTGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGt                                                                                                                                                                                       >  2:182059/1‑251 (MQ=255)                                                             gCCGACTTCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAAGCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGa                                                                                                                                                                                   <  1:300740/251‑1 (MQ=255)                                                             gCCGACTGCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGa                                                                                                                                                                                   <  1:248510/251‑1 (MQ=255)                                                                gACTTCTTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACg                                                                                                                                                                                 >  1:404375/1‑250 (MQ=255)                                                                       tG‑CGCGAAGGTTTCAACGTTGTCACGCCGAACAAAAAGGACAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATACCTCTATGAC‑ACCAACGTTGGGACTAGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCAt                                                                                                                                                                          <  2:473497/250‑1 (MQ=255)                                                                             gAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGTTGAAGTTCTCCGGCTTTCTTTCAGGTTCGCTTTCTTATATCTTCGGCATGTTAGACGAAGGCTTGAGttt                                                                                                                                                                    >  2:423138/1‑251 (MQ=255)                                                                             gAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGttt                                                                                                                                                                    >  2:223725/1‑251 (MQ=255)                                                                              aaGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTACGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTtc                                                                                                                                                                   >  1:328600/1‑251 (MQ=255)                                                                                    tCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTTTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCg                                                                                                                                                               >  2:238844/1‑249 (MQ=255)                                                                                            gTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTTGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGGTAGACGAAGGCATGAGTTTCTCCGAGGCGActa                                                                                                                                                      <  1:438963/250‑3 (MQ=255)                                                                                               aCGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATCTCTTCGGCAAGTTAGACGAAGGCATGAGTTTATCCGAGGCGActactct                                                                                                                                                  >  2:246558/1‑245 (MQ=255)                                                                                                cGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGActactctg                                                                                                                                                 <  2:404375/251‑8 (MQ=255)                                                                                                 cccGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCAGGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGAGGAAGTTCTTCGGCATTCTTTCAGGTTCGCTTTCGTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTGCTCCGAGGCGACTACTCTgg                                                                                                                                                <  2:499677/250‑1 (MQ=255)                                                                                                     aaCAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGcgcg                                                                                                                                            >  2:291079/1‑251 (MQ=255)                                                                                                       cAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACAGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGGAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCggg                                                                                                                                          >  2:59206/1‑251 (MQ=255)                                                                                                         aaaaGTCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCATGTTCGCTTTCTTATATCTTCGGCa                                                                                                                                                                                          <  1:279767/201‑1 (MQ=255)                                                                                                         aaaaGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCa                                                                                                                                                                                          >  2:279767/1‑201 (MQ=255)                                                                                                             ggccacagccTCGTCGATGGAGTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATTAC‑ACCATCGTGGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATTTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCGTTTAGGTTCGCTTTCTTATATCTTCGGTAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATgg                                                                                                                                    <  1:300402/243‑1 (MQ=255)                                                                                                                 aaCACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAACCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTta                                                                                                                                <  1:441224/251‑1 (MQ=255)                                                                                                                    acCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGCTGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATAc                                                                                                                             >  1:70531/1‑251 (MQ=255)                                                                                                                             aTGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCgg                                                                                                                     >  1:102928/1‑250 (MQ=255)                                                                                                                                        ccATCTGTTGCGTCAGGCGGCGGAAAAATCGCGGCGTAAATTCCGCTATGAC‑TCCAACGTGGGTGCTGGATTACCGGTTATTGAGATCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGGTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAgatg                                                                                                         <  2:73504/251‑1 (MQ=255)                                                                                                                                             tgTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCtt                                                                                                    >  1:168246/2‑251 (MQ=255)                                                                                                                                              gTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAAGCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCtt                                                                                                    <  2:70853/250‑1 (MQ=255)                                                                                                                                              gTTGCGTCATGCGGCGGAAAAATAGCGGCGTAACTTCCTCTATGAC‑ACCACCGTTGGGGCTGGATTACCGGATATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGCCGCAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATCCCGAACCGGATCCGCGAGATGATCttt                                                                                                   >  2:176078/1‑251 (MQ=255)                                                                                                                                                          ggcggAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGCTTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGCACCGGATCCGCGAGATGATCTTTCTGGTATggctg                                                                                       >  1:39345/1‑248 (MQ=255)                                                                                                                                                             ggAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCggctactctg                                                                                                                                                 >  2:62215/1‑181 (MQ=255)                                                                                                                                                             ggAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCggctactctg                                                                                                                                                 <  1:62215/190‑10 (MQ=255)                                                                                                                                                                        ggcgTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGCAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACAGAACCGGATCCGCGAGATGATCTTTCTGGTATGGCTGTAGCGCGTaagct                                                                          >  1:287142/1‑247 (MQ=255)                                                                                                                                                                           gTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTg                                                                      >  2:100529/1‑251 (MQ=255)                                                                                                                                                                                                 aCGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGa                                                 <  1:204196/251‑1 (MQ=255)                                                                                                                                                                                                         gCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACt                                         <  1:246558/251‑1 (MQ=255)                                                                                                                                                                                                         gCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAGTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACt                                         <  2:102928/251‑1 (MQ=255)                                                                                                                                                                                                            tgATTACCGGGTAGTGATAACCGTCAAAATCTGCTCAATTCTGGGTATGAATTGATGAAGTTCTCCGGTATTCTTTCATGTTCGCTTTCTTATATCTTCGGTAATTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGAGTTATTCCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAATGGGACGTCAACTg                                        <  2:266613/248‑1 (MQ=255)                                                                                                                                                                                                              aTTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGc                                    >  1:413575/1‑251 (MQ=255)                                                                                                                                                                                                                        aTTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGata                           <  1:291079/250‑1 (MQ=255)                                                                                                                                                                                                                           gagaACCTGCAAAAACTGCTCAATGCTTGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCTCTTTCTTATATCTTCGGAAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATattg                        <  1:431825/250‑1 (MQ=255)                                                                                                                                                                                                                              aaCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATattgaagt                    >  2:173796/1‑249 (MQ=255)                                                                                                                                                                                                                              aaCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAAAGGGACGTGAACTGGAGCTGGCGGATattgaaat                    >  1:416639/1‑251 (MQ=255)                                                                                                                                                                                                                                ccTGCAAATGCTGCTCAATGCTGGTGATTAATTGATGTAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTTTTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATattgaaattg                  <  1:483905/251‑1 (MQ=255)                                                                                                                                                                                                                                   gCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTACGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAAc               >  1:369761/1‑251 (MQ=255)                                                                                                                                                                                                                                    cAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGTTTCGCTTTCTGATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAAcc              <  1:305594/251‑1 (MQ=255)                                                                                                                                                                                                                                        aTCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCtgtg          >  1:310961/1‑251 (MQ=255)                                                                                                                                                                                                                                            gCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGtgctgc      <  2:339846/251‑1 (MQ=255)                                                                                                                                                                                                                                                aaTGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGc  >  1:141597/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      CATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTG‑CGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGAC‑ACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGGCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGC  >  NC_000913/1969‑2456

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑