Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,456,338	A→G	N185S (AAC→AGC)	mpl →	UDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,456,338	0	A	G	100.0%	126.8 / NA	38	N185S (AAC→AGC)	mpl	UDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/17);  total (21/17)

GGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAACAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAACCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCG  >  NC_000913/4456095‑4456580                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   gggTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAACAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGa                                                                                                                                                                                                                                              <  2:294308/250‑1 (MQ=255)        ggCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAACAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGTTCCTCAACAGCCTTGAGTTCGAt                                                                                                                                                                                                                                       <  1:106280/250‑1 (MQ=255)                     acacacGGCAAAACCCCCCCCGCGGGAATGGCGACCTGGATTCTGGAACAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCTATCACGCCGATATCtt                                                                                                                                                                                                                         <  2:121576/251‑1 (MQ=255)                      cactcGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAACAGTGTGGTTACAAACCGTGCTTTGTAAGCGGCGGTGTGCCGGGGAACTTTGAGGTTTCTGCTCATCTGGGCGAAAGCGACTTCTTTGGTATCGAAGCGGATGAGTATGACGGCGCCTTCTTCGACAAACGCTCGAAATGTGGTCATGACTGCCCGCGTACGCTGACCCTCAACAGCCTTGAGTTCCACCACGCCTATATCtga                                                                                                                                                                                                                        <  2:30322/247‑3 (MQ=255)                      cacaaGGCAAAATCACCACCGCGGGAATGGCGACCTGGATGCTGGAACAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTTTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGAACGTTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCttt                                                                                                                                                                                                                        <  1:58045/251‑1 (MQ=255)                        cacGGCAAAACCACCCACGCGGGAATGGCGACCTGGATTCTGGAACAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGCTCTCTTTg                                                                                                                                                                                                                       >  2:327261/1‑250 (MQ=255)                                          gcgGGAATGGCGACCTGGATTCTGGAACAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACCAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATcc                                                                                                                                                                                                     >  1:287064/1‑250 (MQ=255)                                                       ccTGGATTCTGGCACAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTccac                                                                                                                                                                                       >  2:33425/1‑251 (MQ=255)                                                           gATTCTGGAACAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATc                                                                                                                                                                                   >  1:361815/1‑251 (MQ=255)                                                              tCTGGAACAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTg                                                                                                                                                                                 <  1:199013/250‑1 (MQ=255)                                                                       gtgtGGTTACAAACCGGGCTTTGTAAGCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCg                                                                                                                                                                       <  2:106884/251‑1 (MQ=255)                                                                                     cGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGc                                                                                                                                                         >  1:5943/1‑251 (MQ=255)                                                                                      gggCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGcc                                                                                                                                                        <  2:125801/251‑1 (MQ=255)                                                                                             gTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGtattat                                                                                                                                                 >  1:255514/1‑251 (MQ=255)                                                                                                    gcggTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCg                                                                                                                                          <  2:82483/251‑1 (MQ=255)                                                                                                         gtgCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGaaa                                                                                                                                      >  1:59085/1‑250 (MQ=255)                                                                                                                         aGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCt                                                                                                                         >  2:118987/1‑247 (MQ=255)                                                                                                                         aGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCt                                                                                                                         <  1:118987/247‑1 (MQ=255)                                                                                                                            tttCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAAc                                                                                                                    >  2:280745/1‑249 (MQ=255)                                                                                                                               cGGCTCATTTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGACACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGCCCTCCACCTGAAAcggccc                                                                                                               >  2:299990/1‑246 (MQ=255)                                                                                                                                 gCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCAt                                                                                                             >  2:206096/1‑251 (MQ=255)                                                                                                                                              aaaGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGgctgc                                                                                                >  1:85977/1‑251 (MQ=255)                                                                                                                                                 gCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTATGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTgg                                                                                             <  2:255514/251‑1 (MQ=255)                                                                                                                                                   gACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTggcg                                                                                           >  2:321952/1‑249 (MQ=255)                                                                                                                                                         tttGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACACACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCCTGGCGATGGGCTGCTGgagcgagc                                                                                      >  1:148032/1‑250 (MQ=255)                                                                                                                                                             ttATCGAAGCGGATGAGTATGACTGCGACTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAgg                                                                                   >  1:129945/1‑249 (MQ=255)                                                                                                                                                                 cGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAGGAGCtgg                                                                             >  2:68136/1‑251 (MQ=255)                                                                                                                                                                                    tNCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTg                                                          >  1:98747/1‑251 (MQ=255)                                                                                                                                                                                             ttcGACAAACGCTCTAAATTTGTTCATTACTGACCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCAAGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCAGGCAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGa                                                  >  1:224776/1‑250 (MQ=255)                                                                                                                                                                                                                    tCATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCg                          >  1:75168/1‑251 (MQ=255)                                                                                                                                                                                                                     cATTACTGCCCGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTGCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCGGCGGGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGa                         <  2:113097/251‑1 (MQ=255)                                                                                                                                                                                                                            tcccGCGTCCGATGACCATCAACAGCCTTGAGTGCGATCACGCCGATATCTTTGACGACCGGAAAGCGATCCAGAAACAGTTCCGCCAGCTTGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAGGAGCTGGTGGGTGAGCATGGTCACTGGCAGGCGAAAAAGCTGCCCACCGATGCTTCCGATTGGGa                   <  2:48602/249‑1 (MQ=255)                                                                                                                                                                                                                              ccGCGTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGt                >  2:154701/1‑251 (MQ=255)                                                                                                                                                                                                                               cgcgTACGCTGATCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGtt               <  2:70663/251‑1 (MQ=255)                                                                                                                                                                                                                                          aTCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATgg    <  2:129486/251‑1 (MQ=255)                                                                                                                                                                                                                                           tCCTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGc   >  2:207172/1‑251 (MQ=255)                                                                                                                                                                                                                                            ccTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAGGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCg  <  1:88317/251‑1 (MQ=255)                                                                                                                                                                                                                                            ccTCAACAGCCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGc   <  1:242748/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   GGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAACAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGCTCTAAATTTGTTCATTACTGCCCGCGTACGCTGATCCTCAACAACCTTGAGTTCGATCACGCCGATATCTTTGACGACCTGAAAGCGATCCAGAAACAGTTCCACCATCTGGTGCGTATCGTTCCGGGGCAGGGCCGTATTATCTGGCCGGAAAATGACATCAACCTGAAACAGACCATGGCGATGGGCTGCTGGAGCGAGCAGGAGCTGGTGGGTGAGCAGGGTCACTGGCAGGCGAAAAAGCTGACCACCGATGCTTCCGAATGGGAAGTTTTGCTGGATGGCG  >  NC_000913/4456095‑4456580

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