Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,269,449	C→T	C12C (TGC→TGT)	aphA →	acid phosphatase/phosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,269,449	0	C	T	100.0%	134.3 / NA	40	C12C (TGC→TGT)	aphA	acid phosphatase/phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/15);  total (25/15)

ATGTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCGAAGA  >  NC_000913/4269208‑4269694                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       aTGTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTc                                                                                                                                                                                                                                               >  2:137053/1‑251 (MQ=255)             aTCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGt                                                                                                                                                                                                                                   >  2:130466/1‑251 (MQ=255)                gAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTcc                                                                                                                                                                                                                                <  1:512148/251‑1 (MQ=255)                 aGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCg                                                                                                                                                                                                                               >  1:243351/1‑251 (MQ=255)                      taaACATTTTTTATACAAAAAAAGACAGGAACGTATTTCCTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAACGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGtt                                                                                                                                                                                                                          >  2:101779/1‑251 (MQ=255)                       aaaCATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAACAGTTTCAACATACTGACTATTTAGGGAAAAATCTGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTg                                                                                                                                                                                                                         >  2:435958/1‑251 (MQ=255)                                 tataCAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTc                                                                                                                                                                                                               >  1:226898/1‑251 (MQ=255)                                      aaaaaaaGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCAtctc                                                                                                                                                                                                          >  1:581035/1‑251 (MQ=255)                                                                   tAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGtt                                                                                                                                                                             >  1:387118/1‑251 (MQ=255)                                                                            tcatcCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAACTCTCACAATCAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTACCGTTGCCAGGCtt                                                                                                                                                                    >  2:578881/1‑251 (MQ=255)                                                                              atcCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGc                                                                                                                                                                  <  1:126131/251‑1 (MQ=255)                                                                                               aaaatCCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTGGCCCTGGCCTCATCTCCGTCACCGCTTAACCCTGGGACTAACGTGGCCAGGCTTGCTGA‑ACAGGCGCCcattc                                                                                                                                                <  2:225777/251‑1 (MQ=255)                                                                                                         gATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGTAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTc                                                                                                                                      <  1:6061/251‑1 (MQ=255)                                                                                                            aaTCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGt                                                                                                                                   >  2:587027/1‑251 (MQ=255)                                                                                                            aaTCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGt                                                                                                                                   >  1:473807/1‑251 (MQ=255)                                                                                                                 ttccttTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACCGGCACCCATTCATTGGGTTTCGGTCGcac                                                                                                                              >  2:181992/1‑251 (MQ=255)                                                                                                                 ttccttTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGcac                                                                                                                              >  1:106611/1‑251 (MQ=255)                                                                                                                        ccccGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTg                                                                                                                        >  1:452066/1‑250 (MQ=255)                                                                                                                               aGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCTCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑CCAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGc                                                                                                                >  1:131192/1‑251 (MQ=255)                                                                                                                                           ggTCAGATGGGTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGt                                                                                                    >  1:313160/1‑251 (MQ=255)                                                                                                                                              cAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCg                                                                                                 <  2:132876/251‑1 (MQ=255)                                                                                                                                                             acaaATATCACAATAAACAATTTATACATACTGACTATTTAGGGATAAATATGCGCAAGATCACACAGGCAATCAGTGCGGTTTGTTTATTGTCCGCTCTAAACAGATCCGCTGTTGACCTGGCCTCATCTCCTTCACCGCTTAAACCTGGGACTAACGTTGCACGGCTTGCTGATAC‑GGCACCAATTCATTGGGGTTAGGTCGCACAAATTGAAAATAGACTCGCTGGGCGTCAGCCAATGGAGGCgggg                                                                                  >  2:525598/1‑251 (MQ=255)                                                                                                                                                                      acaATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATc                                                                         <  1:428206/251‑1 (MQ=255)                                                                                                                                                                                gTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACgg                                                               <  2:166477/251‑1 (MQ=255)                                                                                                                                                                                  ttCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGt                                                              >  1:10722/1‑250 (MQ=255)                                                                                                                                                                                        aTACTGACTATTTAGGTAAAAATATGCGCAAGACCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTGCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACttttt                                                       <  2:10512/251‑1 (MQ=255)                                                                                                                                                                                               cTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTTTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCCGGGCGTCCGCCAATGGCGGTGGGGGTTGATTTCCATTACACCGTACTTTTTTCCAGTc                                                >  2:114483/1‑251 (MQ=255)                                                                                                                                                                                                 aTTTAGGGAAAAATATGCGCAAGATCACAAAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGGTAACGATGACACGGGACTTTTTTCCAGTccc                                              >  2:341945/1‑250 (MQ=255)                                                                                                                                                                                                   ttAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTTCAGTCCggg                                            >  1:151302/1‑251 (MQ=255)                                                                                                                                                                                                           aaaTATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACACATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGGAACGGTACTTTTTTCCAGTCCGGGCTTCTGgc                                    >  1:429898/1‑251 (MQ=255)                                                                                                                                                                                                           aaaTATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGTGCTTCTGgc                                    <  1:130466/251‑1 (MQ=255)                                                                                                                                                                                                                      aGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCaaaaaa                         <  1:76198/251‑1 (MQ=255)                                                                                                                                                                                                                            cacaGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTtc                   <  1:543839/251‑1 (MQ=255)                                                                                                                                                                                                                              caGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCGGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATTACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTtctc                 <  1:341945/251‑1 (MQ=255)                                                                                                                                                                                                                               aGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCg                <  2:452066/251‑1 (MQ=255)                                                                                                                                                                                                                                      caggtccGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTACCCCGGGGACTAACGTTGCCAGGCGTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGGCGCAAAAAGTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTGTGAGATCGATGACACGGTACTTTTGTCCAGTCCGGGCTTCGGGCGCGGCAAAAAAGCCTTCTCTCCAGaa          <  2:547261/245‑1 (MQ=255)                                                                                                                                                                                                                                      cAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCCGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGaa          >  1:483215/1‑250 (MQ=255)                                                                                                                                                                                                                                         ggCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCCCCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCGCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCg      <  2:226898/250‑1 (MQ=255)                                                                                                                                                                                                                                          gCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTTGCCAGAAAGcc      >  1:390505/1‑249 (MQ=255)                                                                                                                                                                                                                                          gCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATATCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCCTTCATTGGGGTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGAGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGGTTCTGGCGACGCAAAAAACCCCTCTCGTCAGGAAGCg      >  1:515817/1‑250 (MQ=255)                                                                                                                                                                                                                                             gtttgtttATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTTATATCGATGACACGGTACTTTTTTCCAGTCCGGGCGTCTGGCGCGGCAAAAAAACCTTTTCGCCAGAAAACGa     >  1:445814/1‑248 (MQ=255)                                                                                                                                                                                                                                             gtttgtttATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCGAAGa  <  1:557970/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       ATGTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGA‑ACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCGAAGA  >  NC_000913/4269208‑4269694

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