Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,436	T→C	T38T (ACT→ACC)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,436	0	T	C	100.0%	126.8 / NA	44	T38T (ACT→ACC)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/13);  total (31/13)

TACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTG‑TCATCAGGCCAAAACTCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGA  >  NC_000913/4617194‑4617681                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       tACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGCCGAGAAAGTAATTGCTCTGTG‑TCCTCAGGCCAAACCCCCGG‑TCgg                                                                                                                                                                                                                                               >  2:152945/1‑251 (MQ=255)     gTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCCCCCTGAATGACGACGACACCGACGATAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCaag                                                                                                                                                                                                                                           <  1:174703/251‑2 (MQ=255)        aTGCAAACTTGTAAGTAGATTTCCTTAATTTTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATAcc                                                                                                                                                                                                                                        <  2:385390/251‑1 (MQ=255)           cAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGTAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACAACCCTGATTGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TAATAAGGCCAAAACCCCGG‑TCGGCAATAccgcc                                                                                                                                                                                                                                     >  2:281002/1‑251 (MQ=255)                            ttCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATATATcc                                                                                                                                                                                                                    >  2:455910/1‑251 (MQ=255)                             tCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGAACTGACCACCCTGAATGACTACGACACCGACGCGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTTTATCTATCCt                                                                                                                                                                                                                   >  1:435031/1‑251 (MQ=255)                              cctaaattgggatttatATCGAGGTGTGTTGCGGAGTAGATTTTAGAATACTAAAAAACTCGAAAGGTAAATTTTATGGGCGACAAGCCAGGAGAATAAAATGTCTAATCTGAAAGCAAGCACCCGGCGTGCACTGAAATTGAGGGCCCTACCCCCCCTGAATGACGACGCCCCCGACGAAAAAGAATTTGCTCTGTT‑TCTTCAGCCCAAACCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCt                                                                                                                                                                                                                   <  2:205971/234‑1 (MQ=255)                              ccTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAAAATACTAACAAACTCGAAGGGTGATTTTTATTGGCGACAAGCCAGGAGAATGAAATGATGGATCGGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGCCCGCACCACCCTGAATGACGACGACCCCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCAAAACCCCCGG‑TCGTCAATACCGCCGCTATCTGTATCTATCCt                                                                                                                                                                                                                   <  1:239559/250‑1 (MQ=255)                                    ttGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCttt                                                                                                                                                                                                             >  1:70633/1‑250 (MQ=255)                                         aTGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCg                                                                                                                                                                                                       >  1:191879/1‑251 (MQ=255)                                                                   aGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAACCACTGAAAGAGc                                                                                                                                                                              >  2:504424/1‑250 (MQ=255)                                                                                  aaCAAACTCGCAAGGTGAATTTTATTGGCGCCAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAAt                                                                                                                                                              >  2:287274/1‑251 (MQ=255)                                                                                             aaGGTGGATTTTATTGGCGACAAGCCAGGAGAATGCCATGACTGATCTGAAAGAAAGCAGCCTGCGTGCACTGAAATTGATGGACCAGACCACCCTGAATGACGACGACCCCGCAGAGAAAGTAATTGATCTGTG‑TCATCAGGCCAAAACCCCGG‑ACGGCACTACCGCCGCTATCTGTCTCTCACCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACACCGGAAATCCGTATTGCt                                                                                                                                                    >  2:196459/1‑250 (MQ=255)                                                                                               ggTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAAc                                                                                                                                                                                         >  2:135888/1‑211 (MQ=255)                                                                                               ggTGAATTTTATGGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGCCCCCCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTT‑TCATCAGGCAAAACCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAAc                                                                                                                                                                                         <  1:135888/211‑1 (MQ=255)                                                                                                           tGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTc                                                                                                                                     >  1:106792/1‑251 (MQ=255)                                                                                                               gACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTaa                                                                                                                                             >  1:236908/1‑239 (MQ=255)                                                                                                                    gCCAGGCGACTGCAATGACTGCTCTGAAAGCCAGCAGCCTGCGCGCACTGAAATTGCTGGACCTGACCACCCTGAATGACGACGACACCGACGAGAACGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCACTACCGCCGCTATCTGTATCTATCATCGCTTTATCCCGATTGCTCGCAAACAACTGAAAGAGCAGGGCGCCCCGGAAATCCGTATTGCTACGGTACCCAACTTCCCGCACGGt                                                                                                                            >  2:233540/1‑251 (MQ=255)                                                                                                                                aaatgaCTGATCTGAAAGCAAGCAGCCGGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCCATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCCAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCACCTTCCCACCCGGTAACGCtgacatc                                                                                                                >  2:325566/1‑251 (MQ=255)                                                                                                                                     attGATCTGAAAGCAAGCAGCCTGCGTGCACTGAATTTTATGGCCCTGACCCCCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCCGGCCAAACCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGTTAACGAtgacatcgaaat                                                                                                           <  1:152945/249‑4 (MQ=255)                                                                                                                                       ggATCTAAAAGCAAGCAGCCTGCGTGCACTAAAATTGATGGACCTGCCCCCCCTGAATGACGACGACACCGAGGAAAAAGTAATTGCTCTGTG‑TCATCAGGCAAAACCCCCGG‑TCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGAATGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGacatcgaaatcg                                                                                                         <  1:9279/250‑1 (MQ=255)                                                                                                                                               aagcaagcAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGc                                                                                                  >  2:388693/1‑250 (MQ=255)                                                                                                                                                   aggcacccggCGTGCACTAAAATTGATGGACCGGACCACCCTAAATGACGACGACCCCGACGAAAAAGTAATTGCTCGGTT‑TCATCAGGCCAAACCCCCGT‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                                                                                             <  1:300934/242‑1 (MQ=255)                                                                                                                                                      cagcCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAg                                                                                                                                                                            >  2:207519/1‑169 (MQ=255)                                                                                                                                                      cagcCTGCGTGCACTAAAATTGATGGCCCTGACCCCCCTGAATGACGACGACCCCGACGAGAAAGTAATTGCTCGTTG‑TCATCAGGCCAAACCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAg                                                                                                                                                                            <  1:207519/169‑1 (MQ=255)                                                                                                                                                        cccggCGTGCACTGAAAGTGATGGACCGGCCCCCCCTAAAGGACGACGACCCCGAGGAGAATGTAATTGTT‑TGTGGTCATCGGGCCAAACCCCCGG‑TCGGCAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGt                                                                                        <  2:106792/247‑1 (MQ=255)                                                                                                                                                             cGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGCCATCGAAATCGCGCTGGCAGAAACCCGTGCGGc                                                                                   >  2:235767/1‑251 (MQ=255)                                                                                                                                                             cGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGc                                                                                   >  2:86634/1‑251 (MQ=255)                                                                                                                                                                gCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAAt                                                                                >  1:52627/1‑251 (MQ=255)                                                                                                                                                                           gATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTg                                                                     >  1:195841/1‑251 (MQ=255)                                                                                                                                                                           gATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTg                                                                     >  1:190287/1‑251 (MQ=255)                                                                                                                                                                                    gACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCCTCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGCTTGCTCGCAAAACACTGACAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTACCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCcgatgc                                                              >  2:397292/1‑248 (MQ=255)                                                                                                                                                                                      ccaccCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGtt                                                          >  2:29222/1‑251 (MQ=255)                                                                                                                                                                                          ccTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACg                                                      >  1:33925/1‑251 (MQ=255)                                                                                                                                                                                                       gacACCGACGAGAAAGTAATTGCTCTGTG‑TCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGt                                          >  2:271430/1‑250 (MQ=255)                                                                                                                                                                                                                         attgctcTGTG‑TCATCAGGCCAAACCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCggg                       <  1:378491/245‑1 (MQ=255)                                                                                                                                                                                                                            gctcTGTG‑TCATCAGGCCAAA‑CCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGAGGTGTTCCCGTAACGCGCGCTGATGGCGGGGAAc                   >  2:453024/4‑251 (MQ=255)                                                                                                                                                                                                                               cctggatcagcacgCCAAACCCCCGG‑CCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACg                  <  2:195841/238‑1 (MQ=255)                                                                                                                                                                                                                                   gtgacaaaGGCCAAAACCCCGGGTCGGAAAAACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAg                <  1:271430/243‑1 (MQ=255)                                                                                                                                                                                                                                      catcaGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAgagggg           >  2:171057/1‑246 (MQ=255)                                                                                                                                                                                                                                             ccAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGtttt    >  1:169116/1‑251 (MQ=255)                                                                                                                                                                                                                                               aaaaCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGCCATCGAAATCGCGCTGGCAGAAACCCGTGCGGCCATCGCCTCCGGTGCCGATGAAGTTGACGTGGTGTTCCCGTCCCGCGCGCTGATGGCGGGTAACGAGCCGGGTGGTTTTgt  >  2:222495/1‑250 (MQ=255)                                                                                                                                                                                                                                               aaaaCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGACCGGGCAGGGTGGGTTTgt  >  2:140474/1‑250 (MQ=255)                                                                                                                                                                                                                                               aaaaCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGAACGAGCGGGGTGGTTTTg   >  1:290030/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       TACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTG‑TCATCAGGCCAAAACTCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGA  >  NC_000913/4617194‑4617681

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