Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,206,024	G→A	I503I (ATC→ATT)	purH ←	bifunctional AICAR transformylase/IMP cyclohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,206,024	0	G	A	100.0%	148.8 / NA	47	I503I (ATC→ATT)	purH	bifunctional AICAR transformylase/IMP cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/16);  total (31/16)

CAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACGATAGCGTTGGATTTCACAAACTTCGCCACCTTC  >  NC_000913/4205787‑4206270                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        cAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGTCGGAAGTGGCGCATGTCGGTGAATAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGcc                                                                                                                                                                                                                                              >  2:326986/1‑249 (MQ=255)        aaCGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCCTGTCGGTGAAGAGCATCGCAATAACGTGCTCTT‑CCGCGGCGGCAATCACTTCGTCATAACGAATAGAACCGCCAGGCTggct                                                                                                                                                                                                                                     >  2:361311/1‑249 (MQ=255)                      cgGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGAtt                                                                                                                                                                                                                                    >  1:46882/1‑238 (MQ=255)                                                      gCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCgg                                                                                                                                                                                       >  2:86210/1‑251 (MQ=255)                                                        aacaaAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGa                                                                                                                                                                                      >  1:333600/1‑250 (MQ=255)                                                                       aaCCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACAGCCAGGCTGGATTACGCAAGTCACGCCCGCAGCTGCGGCGGCACCAATACCGTCGCGGAACGGGAAAAACGCGTc                                                                                                                                                                      >  2:160842/1‑251 (MQ=255)                                                                        aCCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATTTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTAATTACGAATAGAACCGCCAGGCTGGCTTACGCAGGTCACGCCCCCAGCGGCGGCGGCATCAATAACGTCGCGGGACGGGCAGAAcggggc                                                                                                                                                                      >  2:55702/1‑246 (MQ=255)                                                                              ggcgACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACgg                                                                                                                                                                                  >  1:72400/1‑232 (MQ=255)                                                                                 gACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCCCGGGACGGGAAGAACGCGTCAGAAGCCATc                                                                                                                                                            >  2:349415/1‑251 (MQ=255)                                                                                            tttCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCCGAAGCCATTGAGGAACCttt                                                                                                                                                 >  1:327373/1‑251 (MQ=255)                                                                                                       gcgcgtgccccgccccgccGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCATTCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGc                                                                                                                                      <  1:74412/238‑1 (MQ=255)                                                                                                         gtgtgcccccccccgccGTTTCCAATCACTAATTCTTCCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGAAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCt                                                                                                                                    <  2:288334/240‑1 (MQ=255)                                                                                                             gtccacccccgccTTTACCAACCGCGAATTCTTTCATCTATTGCTCCATTAGTGGCGGAAGTGGCGCATGTCGGTGAAGATAACCGCAATGCCGTGCACGT‑CGGCAGCGGCAATCATTTCGTCATCACGAAAAGAACCGTCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATCCCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCa                                                                                                                                 <  1:188584/243‑1 (MQ=255)                                                                                                              cTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGACGAGCATCGCAATACCGTGCTCGT‑CGGCGGTGGCAACCACTTCGTCATCACGACTAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCAGCAACAATACCGTCGCGGAACGGGGAGAACGCCGCAGAAGCCACTGAGGGACCTTTCACTTCCAGGCCTTCacc                                                                                                                               >  2:177048/1‑249 (MQ=255)                                                                                                                  cgcccgccGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGGAACTTTCCATTCCAGGCCTTCATCCGc‑c                                                                                                                          >  2:113218/1‑251 (MQ=255)                                                                                                                       gccgTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCttt                                                                                                                     >  2:279968/1‑251 (MQ=255)                                                                                                                        ccgTTACCAATCACTAATACTTTCATCTCTTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGACAGGCTGGATTACGCAGGTCACGCACGCAGCAGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGTCCTTCATCGGC‑CGCttt                                                                                                                     >  1:339031/1‑250 (MQ=255)                                                                                                                         cgTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTCTCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGAATACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGCAGAACGCGTCAGAAGCCATCAAGGAACCTTTCACTTCCAGGCCTTCATAGGC‑CGCTTTa                                                                                                                    >  1:46388/1‑250 (MQ=255)                                                                                                                               cAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGCACGGGAAGAACGCGCCAGAAGCCAACGAGGAACCTTTCCCCTCCAGGCCTTCATCGGCACG‑TTTAATACCgg                                                                                                             >  1:308110/1‑251 (MQ=255)                                                                                                                                          aCTTTCATCTATTGCTCCATTAAGGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGc                                                                                                   <  1:175879/250‑1 (MQ=255)                                                                                                                                              tCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCCGGCCTTAATCGGC‑CGCTTTAATACCGGCGAATTTTGCGGAGt                                                                                              >  1:201932/1‑251 (MQ=255)                                                                                                                                                    aTTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACAcg                                                                                         <  1:228544/250‑1 (MQ=255)                                                                                                                                                        cTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCt                                                                                    >  2:297701/1‑251 (MQ=255)                                                                                                                                                          ccATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTcc                                                                                  >  1:7928/1‑250 (MQ=255)                                                                                                                                                           cATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAAAGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCAt                                                                                 >  1:60506/1‑251 (MQ=255)                                                                                                                                                              tAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCt                                                                               <  1:247466/250‑1 (MQ=255)                                                                                                                                                                           ggCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATg                                                                 >  1:252559/1‑251 (MQ=255)                                                                                                                                                                             cgcATGTCGGTGAGGAGCATCGCAATACCGTGCT‑GTGCGGCGGCGGCAATCCCTTCGTCATCACGAATGGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGcc                                                               <  2:156135/251‑1 (MQ=255)                                                                                                                                                                                  gTCGTTGAAGATCATCGCAATCCCGTGCTCGC‑CGGCGGCGCCAACCCCTTCGTCATCCCGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATCCCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATa                                                           <  1:66166/250‑1 (MQ=255)                                                                                                                                                                                     ggTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAACAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTc                                                        >  1:155062/1‑250 (MQ=255)                                                                                                                                                                                      gTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCat                                                      >  1:172837/1‑251 (MQ=255)                                                                                                                                                                                            agCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCAGCGCCAATGCCGATCGTCATATTGtt                                                >  2:298926/1‑251 (MQ=255)                                                                                                                                                                                              cATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGtttt                                              >  1:216272/1‑251 (MQ=255)                                                                                                                                                                                                       aCCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAg                                      <  2:333600/250‑1 (MQ=255)                                                                                                                                                                                                        ccGTGCTCTG‑CGGAGGCGTCAATCATTTCGTCATCACGAATAGACCCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGAc                                    <  1:7985/251‑1 (MQ=255)                                                                                                                                                                                                         cGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACg                                   >  1:333561/1‑251 (MQ=255)                                                                                                                                                                                                               ctt‑cgGCGGCGGCAATCACTTCGTCATCACAAACAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACGATAg                               <  2:155062/247‑1 (MQ=255)                                                                                                                                                                                                                gt‑cgGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTCATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACGATAGCGt                            >  1:248011/1‑251 (MQ=255)                                                                                                                                                                                                                 t‑cgGCGGCGGCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACGATAGCGtt                           >  1:138188/1‑251 (MQ=255)                                                                                                                                                                                                                   cggcggcggCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCTGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACGATAGCGTTg                          <  2:67695/251‑1 (MQ=255)                                                                                                                                                                                                                   cggcggcggCAATCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACGATAGCGTTg                          <  2:60506/251‑1 (MQ=255)                                                                                                                                                                                                                   cggcggcggCAATCACTTCGTAATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACGATAGCGtt                           <  1:68343/250‑1 (MQ=255)                                                                                                                                                                                                                      gggnggcAATCCCTCCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACGATAGCGTTGGAt                       <  2:373499/250‑1 (MQ=255)                                                                                                                                                                                                                               tCACTTCGTCATCACGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACGATAGCGTTGGATTTCACAAAc              >  1:100974/1‑251 (MQ=255)                                                                                                                                                                                                                                cACTTCGTCATCGCGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGAGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACGATAGCGTTGGATTTCACAAACt             >  2:143507/1‑251 (MQ=255)                                                                                                                                                                                                                                      gtaatcaCGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACGATAGCGTTGGATTTCACAAACTTCGCCa       <  2:184695/248‑1 (MQ=255)                                                                                                                                                                                                                                           caCGAATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACGATAGCGTTGGATTTCACAAACTTCGCCACCTTc  >  2:103922/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        CAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGT‑CGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTTTCACTTCCAGGCCTTCATCGGC‑CGCTTTAATACCGGCGATTTTTGCGGAGTACACGCGGCTCATCTGGCCCGCGCCAATGCCGATAGTCATATTGTTTTTGGCATAGACGATAGCGTTGGATTTCACAAACTTCGCCACCTTC  >  NC_000913/4205787‑4206270

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