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%	105.8 / NA	39	T38T (ACT→ACC)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/8);  total (31/8)

TGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAAAACTCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTT  >  NC_000913/4617210‑4617673                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                tGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CGGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACAGCCGCTATc                                                                                                                                                                                                                        >  2:43988/1‑251 (MQ=255)                  aaTTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCtt                                                                                                                                                                                                       >  2:215748/1‑251 (MQ=255)                      gtgATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGCCGACACCGACGAGAAAGAAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTAACTGTATCTATCCTCGCTTTATc                                                                                                                                                                                                   >  2:204722/1‑251 (MQ=255)                      gtgATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATc                                                                                                                                                                                                   >  2:243380/1‑251 (MQ=255)                          tgtatATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCg                                                                                                                                                                                                >  2:481760/5‑250 (MQ=255)                            tatATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGAtt                                                                                                                                                                                             >  1:171144/3‑251 (MQ=255)                                   acgtgagtTGCAGATTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATAGGCGACAAGCCAGGAGCATGAAATGACTGATCTGAAAGCAAGAAGCCTGCGAGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCATGCCAAAACCCCGGTC‑GGCAATACCGCCTCTATCTGTCTCTATCCTCGCTTTATCCCTATTGCTCGc                                                                                                                                                                                       >  2:349449/7‑250 (MQ=255)                                                aGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGaa                                                                                                                                                                           >  1:432392/1‑248 (MQ=255)                                                 gTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACAACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAag                                                                                                                                                                         >  2:259405/1‑249 (MQ=255)                                                        ttAGAATACTAACAAACGCGCAAGGTGAATTTTCTTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGACTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAACACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑ACCACTGCAAGAGCAggac                                                                                                                                                                 >  2:403661/1‑248 (MQ=255)                                                           aaaTACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGATAATGAAATGACTGCTCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATAACGACAACACCGACGAAAAAGAAATTGCT‑CTGGGTCATCAGGCCAAACCCCCGGTC‑GGAAAAACCGCCGCTCTCTGTATCTCTCCTCGCATTATCCCGCTTGCTCGCAA‑AAcaatcaaaaaccaggccaccc                                                                                                                                                             >  2:267087/2‑231 (MQ=255)                                                                      gccTCGCAAGGGGTATTTTATTGGCGACAAACCAGGAGAATGAAAAGACTGATCTGAAAGCAAGAAGCCTGCGAGCACTGAAATTGATGGACC‑TGACCACCCTGAATGGCGACGACACCGACGCAAAAGTAATTGCT‑ATGATTCAACAGGCCAAAACCCCGGCC‑GGCAATACCGAACCTATCTGTATCTATCCTCTCTTAATCCCGATTGCTAGCAACGAC‑CTGAAAGAGCAGAGCCCCCCGGAAATcagt                                                                                                                                                  >  2:418260/3‑248 (MQ=255)                                                                             aaGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCt                                                                                                                                            >  1:113940/1‑250 (MQ=255)                                                                              aGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGCCTGATCTGAAAGCACGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTCTCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTAc                                                                                                                                          >  2:398917/1‑251 (MQ=255)                                                                                    aTTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAAc                                                                                                                                    >  1:466497/1‑251 (MQ=255)                                                                                      tttATTGGCGACAAGCAAGGAGAATGAAATGACTGATCTGCAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAAAGAAGGGCACCTCGGAAATCAGTATTGCTACGGTaccc                                                                                                                                   >  2:322958/1‑247 (MQ=255)                                                                                         aggggcgacaggcagggAGAATGAAATGATGGATCGAAAAGCAAGCAGCCTGTGTGCACTAAAATTGTTGGACC‑TGACCACCCTGAATGACGACGACACCGACGAAAAAGTAATTGCT‑CGTTTTCATCAGGCAAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAAc                                                                                                                                <  1:308952/235‑1 (MQ=255)                                                                                                        ggAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACCCGGTAacg                                                                                                                >  2:178662/1‑251 (MQ=255)                                                                                                         gagaAAGAAATGACTGATCGAAAAGCAAGCAGCCTGCGTGCACTGAAATGGATGGACC‑TGACCCCCCTGAATGACGACGCCACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCAAAACCCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTTCGGTAACCAACTTCCCACACGGTAacga                                                                                                               <  1:459037/251‑1 (MQ=255)                                                                                                             agaaaatgatggATCAGAAAGCAGGCAGCATGCGTACACAAAAATTAATGGCCC‑GGCCCCCCCGGAAGGACGACCACCCCGACGAAAAATTATTTGCT‑CTGTGTCAACAGGCAAAACCCCCGGCC‑GTAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACgatgac                                                                                                           <  2:466497/240‑1 (MQ=255)                                                                                                                 aatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGCCGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGAAACCAACTTCCCACACGGTAACGAtgacatcg                                                                                                       >  1:472046/1‑251 (MQ=255)                                                                                                                      atGATCTAAAACCAAGCAGCCTGTGTGCACTGAAATGGA‑GGACCCTCACCCCCCTGAATGACGACGACCCCGCCGAAAAAGTAATTGCT‑CTGTTTCAGCAGGCCAAACCCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATgacatcgaaatc                                                                                                  <  1:285173/250‑1 (MQ=255)                                                                                                                             gAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGCCC‑TGCCCCCCCGAAATGACGCCAACCCCGACGAAAAAGTAATTGCT‑CTGTTTCACCAGACAAAACCCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATcgc                                                                                                <  1:259405/246‑1 (MQ=255)                                                                                                                                          cTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGc                                                                              >  1:139142/1‑251 (MQ=255)                                                                                                                                            gCGTGCACTAAAATTTATGGACC‑TGCCCCCCCTGAATGCCGACGACCCCGACGATAAATTAATTTCT‑CTGTGTCTTTGGGCCAAACCCCCGGTC‑GCCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACTATGACATCGAAATCGCGCTGGCAGAAACCCGTGCgg                                                                            <  1:396346/251‑1 (MQ=255)                                                                                                                                              gTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CCGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGc                                                                           >  2:303330/1‑250 (MQ=255)                                                                                                                                               tGCACTGAGATAGCTGGACC‑TGACCACCCTGAATTCCGACGACACCGACGAGAACGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGCTTGCTCGCAA‑AACCCTGAAAGCGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGCCAGAAACCCGTGCGGcc                                                                          >  2:267937/1‑249 (MQ=255)                                                                                                                                                          tGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACggg                                                              >  1:252660/1‑250 (MQ=255)                                                                                                                                                                          cccTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTaa                                                                                                                                     >  1:346372/1‑165 (MQ=255)                                                                                                                                                                          cccTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTACGGCAATCGCCTACGGTGCCGATGAAGTTGAc                                               >  1:341694/1‑251 (MQ=255)                                                                                                                                                                            cTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAC‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCACTCGCCTCCGGTGCCGATGAAGTTGACGt                                             >  1:301789/1‑251 (MQ=255)                                                                                                                                                                                aggacgacgacACCGACGAAAAAGTAAT‑GCCTCTGTGTCATC‑GGGCAAAACCCCCGTCCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACgtggtg                                         <  1:463292/249‑1 (MQ=255)                                                                                                                                                                                aTGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAAccc                                                                                  >  2:6805/1‑210 (MQ=255)                                                                                                                                                                                    cgacgacACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTcc                                     >  1:121559/1‑251 (MQ=255)                                                                                                                                                                                    cgacgacACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTcc                                     >  1:236254/1‑251 (MQ=255)                                                                                                                                                                                                  gagaaagtatgtgctttttgtCATCACGCAAAACCCCCGGCC‑GACAACACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCt                        <  2:160787/233‑1 (MQ=255)                                                                                                                                                                                                          aattgctcTGTGTCATCAGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCgg                >  1:83253/8‑251 (MQ=255)                                                                                                                                                                                                                      tcatcaGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGACGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAAagagacgg    >  2:254690/1‑247 (MQ=255)                                                                                                                                                                                                                        atcaGGCCAAAACCCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAggtg  >  2:25387/1‑250 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                TGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACC‑TGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAAAACTCCGGTC‑GGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAA‑AACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTT  >  NC_000913/4617210‑4617673

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