Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,625,533	C→T	Y213Y (TAC→TAT)	serB →	phosphoserine phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,625,533	0	C	T	100.0%	117.8 / NA	38	Y213Y (TAC→TAT)	serB	phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/22);  total (16/22)

ACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAA  >  NC_000913/4625290‑4625774                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  aCTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGGGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGcgcg                                                                                                                                                                                                                                            <  2:3875/251‑1 (MQ=255)           gctggcgtggtgttggCGGATGTAATCGAGCGGGCGGTGCGGTGTGAACCCGATTTTCCCGGCAGCCTGCGCAGCCGGGGGGCGACGCTGAAAGTCGCTGCAGCCAATATTATGCAACAGGTGCGTGAAAATCTGCCGCTGATGACAGGCTTAACGCCACTGGTGCTCAAGCCGAAAACGCTGGGCTGGAAAGTTGCGATTGCCGCTGGCGGCTTCACTTTCTTTGCTGAATATCTGCCCGACAAGTTgcg                                                                                                                                                                                                                                  <  2:218505/238‑1 (MQ=255)                 gagaTGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCTACAAGTTGCGCCTGAc                                                                                                                                                                                                                            <  2:193435/251‑1 (MQ=255)                 gagaTGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGAc                                                                                                                                                                                                                            >  1:60692/1‑251 (MQ=255)                             gaatgaaCCGAACGGGCGGTGCGCGGCCAACCCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGGGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAg                                                                                                                                                                                                                 <  2:216138/245‑1 (MQ=255)                             gAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACCAGTTGCGCCTGACCGCCGTGGt                                                                                                                                                                                                                   >  1:216138/1‑248 (MQ=255)                             gAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGc                                                                                                                                                                                                                <  1:144730/251‑1 (MQ=255)                                   aCCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATg                                                                                                                                                                                                           >  1:212627/1‑250 (MQ=255)                                   aCCGAACGGGCGATGCGAGGCGAACTCGATTTTACCGCCATCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCCGACGCCAATATTCTGCAAAAGGTGCGTGAAAATCTGCCGCTGATGCCAAGCTTAAAGCACCTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGACGTGGTAGCCAATg                                                                                                                                                                                                           >  1:108327/1‑250 (MQ=255)                                         cGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGga                                                                                                                                                                                                    <  1:4583/251‑1 (MQ=255)                                               aTGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAAGCTGCCGCTGAGGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCAt                                                                                                                                                                                              <  1:13375/251‑1 (MQ=255)                                                 gcgcGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATgg                                                                                                                                                                                            <  1:264641/251‑1 (MQ=255)                                                   ggggggAACTCGATTTTACCGCCAGCCTGCGCAGCGGGTTGGCGACGCTGAAATGCGCTGTCGCCAATATTCTGCAACAGGTGCGTTAAAAGCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGTCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTATGCTGAAGATCGTCGCGAAAATTTTCGCCTTACCGCGGTGGTAGCCAATGAACTTGAGATCATGGAc                                                                                                                                                                                          <  2:170967/246‑1 (MQ=255)                                                             cGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGCCGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCTCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAAttt                                                                                                                                                                                 >  2:79639/1‑250 (MQ=255)                                                                            ccTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCg                                                                                                                                                                 <  2:49195/251‑1 (MQ=255)                                                                            ccTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCg                                                                                                                                                                 >  1:211966/1‑251 (MQ=255)                                                                                   agccgcgaggCGACGATGAACGGCGCTGACGCCAATATTCTGCAAAAGGTGCGTAAAAATCTGACGCTGATGCCAGACTTAACGCAACTGGTGCTCACGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCACCGTGGTAGCCAATGAACTGGAGATAATGGACGGTCAATTTACCGGCAATGTGATCGGCGACAt                                                                                                                                                          >  2:256155/9‑251 (MQ=255)                                                                                   aGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACAt                                                                                                                                                          <  1:79639/251‑1 (MQ=255)                                                                                      cgcgtgGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGCTCGGCGACATCGt                                                                                                                                                       >  1:188225/4‑251 (MQ=255)                                                                                      cgcgtgGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGt                                                                                                                                                       >  2:5706/4‑251 (MQ=255)                                                                                                             tgccGCCAATATTCTGCAACAGGAGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGCTTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTCGCCAATGTACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAACGCGaaa                                                                                                                                 >  2:88437/4‑250 (MQ=255)                                                                                                              gACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGTCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGaaa                                                                                                                                 <  1:141747/249‑1 (MQ=255)                                                                                                                                          agaaacTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCGTAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCGTGACCGCCGTGGTAGCCAATGAACTGGAGATCAGGGACGGTAAATTTACCGGCAATGTGATCGGCGACAGCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCTCGCAGGGGTATGaa                                                                                                   <  2:39308/246‑1 (MQ=255)                                                                                                                                               cTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATccc                                                                                              >  1:22810/1‑251 (MQ=255)                                                                                                                                                 gCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATccc                                                                                              <  1:19748/249‑1 (MQ=255)                                                                                                                                                       gATGCCAGCCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCGGGAAAGTGGCGATTGCCTCCGGCTGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTGCCGGCAATGTGATCGGCGACATCGTAGACTCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGgcg                                                                                       <  1:140344/250‑1 (MQ=255)                                                                                                                                                         tGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGgcgc                                                                                      >  1:73190/1‑249 (MQ=255)                                                                                                                                                                    aCGCAACTGGTGCTCAAGCTGGAAACGCTGGGCGGGAAAGTGGCGAGTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgat                                                                         <  2:73190/251‑1 (MQ=255)                                                                                                                                                                        aggtggggCTCAAGCTGGAAACGCGGGGGTGGAAAGTGGCGATTGCCCCCGGGGGCTTTACTTTCTTGGCTGAATATCTGCGCGACAAGGTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgtttggc                                                                     <  2:225600/244‑1 (MQ=255)                                                                                                                                                                          cTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTACGGGGGCTTTACTTTCTTTGCTATATATATGCGCGACAAGTTGCGCCTGACCGCCGAGGTAGCCACTGATCTGGAGATACCGGACGGGACATTTACCGGCAATCTGATCGGGGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTTCGAAATCCCGCTGgcgacgcccggggggcttggcgg                                                                   >  1:34049/1‑231 (MQ=255)                                                                                                                                                                                       gggaagcccgtggcggtaaaGTGGCGATGCCCTCTCGCGGCTTTACTTTCTTGGGTGAATTTCGTCGCGACAAGTTGCGGCTGACCGCCGTGGTAGCCACTGAACGGGAGATCATGGACGGGAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAATCTCTGACGCGCCTCTCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGACAATGAc                                                      <  2:89431/234‑1 (MQ=255)                                                                                                                                                                                        ggAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGTCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAAGCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGAcc                                                     <  1:253651/251‑1 (MQ=255)                                                                                                                                                                                                       gAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAg                                      <  1:251909/251‑1 (MQ=255)                                                                                                                                                                                                         aaGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGCTTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCg                                    >  1:228051/1‑251 (MQ=255)                                                                                                                                                                                                                        ccGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGAt                      <  1:5706/250‑1 (MQ=255)                                                                                                                                                                                                                                      ttttttGCTGAAGATCGGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgcc         <  2:211966/246‑1 (MQ=255)                                                                                                                                                                                                                                         tttGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgccaagcc    >  2:62668/1‑251 (MQ=255)                                                                                                                                                                                                                                            gCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgccaagccaa  >  2:100135/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  ACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAA  >  NC_000913/4625290‑4625774

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