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%	147.9 / NA	45	C12C (TGC→TGT)	aphA	acid phosphatase/phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/19);  total (26/19)

GTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGT‑ACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGC  >  NC_000913/4269210‑4269689                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  gTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGtgcg                                                                                                                                                                                                                                         <  2:382639/250‑4 (MQ=255)   aTCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGctc                                                                                                                                                                                                                                      >  2:178337/1‑251 (MQ=255)     ccTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAACAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGctct                                                                                                                                                                                                                                     >  2:144967/1‑250 (MQ=255)          cATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAAc                                                                                                                                                                                                                                 >  1:241968/1‑249 (MQ=255)          cATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAg                                                                                                                                                                                                                               <  1:466115/251‑1 (MQ=255)          cATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAg                                                                                                                                                                                                                               <  1:144967/251‑1 (MQ=255)                          ttttttATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCATTACTTTCTCCGTCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAATGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGcattg                                                                                                                                                                                                               >  2:13383/1‑247 (MQ=255)                                     aaaaaaGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTcc                                                                                                                                                                                                    >  1:451454/1‑251 (MQ=255)                                       aaaaGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCtt                                                                                                                                                                                                  >  1:575703/1‑251 (MQ=255)                                         aaGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTc                                                                                                                                                                                                 >  2:277928/1‑250 (MQ=255)                                                   cGTATTTACTGGGGTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGATTTTGCAACAAATCTCACAATAAATAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTaccc                                                                                                                                                                                      >  2:21537/1‑248 (MQ=255)                                                    gTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAAccc                                                                                                                                                                                     >  1:531382/1‑251 (MQ=255)                                                                  aaaTATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCTGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTg                                                                                                                                                                       >  2:520101/1‑251 (MQ=255)                                                                  aaaTATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTg                                                                                                                                                                       <  1:395145/251‑1 (MQ=255)                                                                        aaTCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCCCAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGc                                                                                                                                                                 >  2:434608/1‑251 (MQ=255)                                                                         atcatcCTGCTTTTCATCACAAAAACCGCCGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCt                                                                                                                                                                >  1:496772/1‑251 (MQ=255)                                                                         atcatcCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCt                                                                                                                                                                >  1:613/1‑251 (MQ=255)                                                                          atatcCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGGTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAAGCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTACCGTTGCCAGGCtg                                                                                                                                                               <  2:88354/249‑2 (MQ=255)                                                                              ccTGCTTTTCATCACAAAAACCGCAGATAATCCGTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCt                                                                                                                                                            <  1:106848/250‑1 (MQ=255)                                                                                 gCTTTTCATAACAAAAACCGCAGATAATCATTAATTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACACATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCAAACCGGCAATCAGTGACGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGACCTGGCCTAATATCCTTAACCGCTTAACCCTGGGACTAACGTTGCCAGGATTGATGaa                                                                                                                                                         >  2:46339/1‑250 (MQ=255)                                                                                        atcaCAAAACCCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACCACTCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCAcc                                                                                                                                                 >  2:48395/1‑251 (MQ=255)                                                                                                  ccGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGTAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTgg                                                                                                                                       <  1:25065/251‑1 (MQ=255)                                                                                                       gATAATCCTTCCTTTCCCCGGCAGCTGGCGGTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATGGGGTTTc                                                                                                                                  <  1:488997/251‑1 (MQ=255)                                                                                                                tccttTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGcac                                                                                                                          >  1:8917/1‑250 (MQ=255)                                                                                                                                tGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTc                                                                                                         >  1:72947/1‑251 (MQ=255)                                                                                                                                 ggCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTc                                                                                                         <  1:277928/250‑1 (MQ=255)                                                                                                                                      tATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAggg                                                                                                   >  1:179074/1‑251 (MQ=255)                                                                                                                                      tATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAgg                                                                                                    <  2:143336/250‑1 (MQ=255)                                                                                                                                         ggTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCg                                                                                                 <  2:451454/250‑1 (MQ=255)                                                                                                                                             aGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCAAACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAACCAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGACAATAGCCTCGCAGGGCGTCCGc                                                                                            >  1:585546/1‑251 (MQ=255)                                                                                                                                                         caacaaATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTgg                                                                                <  1:10723/251‑1 (MQ=255)                                                                                                                                                                   cacaAGAAAAATTTACAACATTCTGACTAGTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAGCAGTTCCGCTGTTGCCCTGTCCTCATCGCCCTCACCGCTTACCCCTGGGACTAACGTAGCCAGGCTTGCTGAACAGGCACCCATACAGTGGGTTTCGGTCGCACAAATTTAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGat                                                                        <  1:558258/249‑1 (MQ=255)                                                                                                                                                                             aGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCGCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACg                                                            <  1:367951/251‑1 (MQ=255)                                                                                                                                                                                           gACTCTTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGACCAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGTAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGCGGGGGTTGATATCGCTGACACGGTGATTTTTTTcc                                               >  2:77124/1‑250 (MQ=255)                                                                                                                                                                                               aTTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGT‑ACTTTTTTCCAGTCCg                                         >  1:323510/1‑251 (MQ=255)                                                                                                                                                                                                     ggAAAAATATGATCAAGATCACTCTGGCAATCAGTTCCGTTTGTTTATTGTTCGCTCTATACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTATCCCTGTGACTCTCGTTGCCAGGCTTTCTGACCATGCCCCCATTCATTGGGTTTCGGGCGCACAAATTGAAAATAGCCTCGCAGTGCGTTCGCCAATGGCGTTGGGGTTTGATATCGATGCCACGGT‑ACTCTTTCCCAGTCCCGGCTTc                                   <  1:72200/251‑1 (MQ=255)                                                                                                                                                                                                        aaaaTATGCGCACGATCACACAGGCAATCAGAGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGGTTGCACGGT‑ACTTTTTTCCCGTCCGGGCTTCTgg                                >  1:197555/1‑251 (MQ=255)                                                                                                                                                                                                            tatGCGCAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACACATTGAAAATAGCCTCGCAGGGCGTCCGCCACTGGCGGTGGGGGTTGATATCGATGACACGGA‑ACTTTTTTCCCGTAGGGGCTTCTGgcggg                            >  1:54808/1‑249 (MQ=255)                                                                                                                                                                                                            taggcgcAAGATCTCAACGGCATTCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCGCATCGCCCTCCCCGCTTAACCCTGGGACTAACGGTGCCAGGCTTGCTGAACAGCCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGGAGGGCGTCCGCCAATGGCGGTGGGGGTTGATAGCGAAGACACGGT‑ACGTTTTTCCAGGCCGGGCTTCTGgcg                              <  2:54808/246‑1 (MQ=255)                                                                                                                                                                                                                cgcAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGGTGACACGGA‑ACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCa                         >  1:147000/1‑250 (MQ=255)                                                                                                                                                                                                                cgcAAGATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGT‑ACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCa                         >  1:125244/1‑250 (MQ=255)                                                                                                                                                                                                                     gaacacacaTGCAATCAGTGCCGGTTGTTTATTGTTCGCTCTAACCAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCGGGGACTAACGTTGCCAGGCTTGCGGAACAGGCCCCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAAGGGCGGTGGGGTTTGATATCGATGAAACGGT‑ACTTTTTTCCAGGCCGGGCTTCTGGCGCGGCaaaaaaa                   <  2:448977/248‑1 (MQ=255)                                                                                                                                                                                                                     gATCACACAGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATTACACGGG‑ACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCaaaaaa                    >  1:32956/1‑250 (MQ=255)                                                                                                                                                                                                                             aGGCAATCAGTGCCGTTTGTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGT‑ACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCg           <  2:179074/251‑1 (MQ=255)                                                                                                                                                                                                                                      gTGCCGTTTTTTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTACCGTTGCCCGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGT‑ACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGc  <  2:320016/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  GTATCCTTACATCGAGTAATAAACATTTTTTATACAAAAAAAGACAGGAACGTATTTACTGGGTTAAATATAATCATCCTGCTTTTCATCACAAAAACCGCAGATAATCCTTCCTTTCCCCGGCAGCTGGCGTTATGGTCAGATGGTTTTTGCAACAAATCTCACAATAAAAAGTTTCAACATACTGACTATTTAGGGAAAAATATGCGCAAGATCACACAGGCAATCAGTGCCGTTTGCTTATTGTTCGCTCTAAACAGTTCCGCTGTTGCCCTGGCCTCATCTCCTTCACCGCTTAACCCTGGGACTAACGTTGCCAGGCTTGCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGT‑ACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGC  >  NC_000913/4269210‑4269689

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