Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,735,191	T→C	Y253H (TAC→CAC)	mepH →	murein DD‑endopeptidase, space‑maker hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,735,191	0	T	C	100.0%	252.3 / NA	63	Y253H (TAC→CAC)	mepH	murein DD‑endopeptidase, space‑maker hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/34);  total (29/34)

GAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGC  >  NC_000913/1735047‑1735338                                                                                                                                                 |                                                                                                                                                    gAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCGCTTCTCTCAGTGAAGACCACTGGc                                                                                                                                               <  2:235299/151‑1 (MQ=255)   aCTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCa                                                                                                                                              <  1:238910/150‑1 (MQ=255)      gAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAgcgc                                                                                                                                          >  1:326325/1‑151 (MQ=255)           aCGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACt                                                                                                                                      >  2:1164103/1‑150 (MQ=255)                     ggTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGgcgc                                                                                                                            <  2:478155/150‑1 (MQ=255)                      gTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCtc                                                                                                                          <  1:925584/151‑1 (MQ=255)                         tttttCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCtcgtc                                                                                                                       >  1:1112180/1‑151 (MQ=255)                          ttttCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCtcgtcg                                                                                                                      >  1:100351/1‑151 (MQ=255)                          ttttCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCtcgtcg                                                                                                                      >  2:770122/1‑151 (MQ=255)                          ttttCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCtcgt                                                                                                                        <  1:44090/149‑1 (MQ=255)                              ccGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGt                                                                                                                   <  2:393644/150‑1 (MQ=255)                              ccGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTa                                                                                                                  <  1:203256/151‑1 (MQ=255)                                gTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAAt                                                                                                                <  2:913886/151‑1 (MQ=255)                                  aCTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATg                                                                                                               >  2:199221/1‑150 (MQ=255)                                  aCTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGa                                                                                                              <  1:163781/151‑1 (MQ=255)                                   cTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGa                                                                                                              <  1:84486/150‑1 (MQ=255)                                    tCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCggg                                                                                                                    <  2:881717/143‑1 (MQ=255)                                    tCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCggg                                                                                                                    >  1:881703/1‑143 (MQ=255)                                      aGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTaa                                                                                                                 >  1:378202/1‑144 (MQ=255)                                      aGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTaa                                                                                                                 <  2:378205/144‑1 (MQ=255)                                      aGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGAcccc                                                                                                          <  2:310297/151‑1 (MQ=255)                                         gACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCaaa                                                                                                       <  2:674663/151‑1 (MQ=255)                                          aCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCaaa                                                                                                       >  1:348833/1‑150 (MQ=255)                                           cgcgGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCaaaaa                                                                                                     <  1:794392/151‑1 (MQ=255)                                           cgcgGCACAGCCGATCATGACGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTAAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCaa                                                                                                        >  1:892436/1‑148 (MQ=255)                                             cgGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAac                                                                                                    >  2:834498/1‑150 (MQ=255)                                                 acaGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTc                                                                                               >  2:844837/1‑151 (MQ=255)                                                  caGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCg                                                                                              >  1:27256/1‑151 (MQ=255)                                                       gATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTaaa                                                                                         <  1:492582/151‑1 (MQ=255)                                                         tCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTaaaa                                                                                        <  1:590159/150‑1 (MQ=255)                                                             gTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACttt                                                                                    <  2:1185322/150‑1 (MQ=255)                                                                  cGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCtg                                                                              >  1:418198/1‑151 (MQ=255)                                                                  cGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCtg                                                                              >  2:1193929/1‑151 (MQ=255)                                                                  cGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCtg                                                                              <  1:1045009/151‑1 (MQ=255)                                                                    tgtATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCtgtt                                                                            >  2:387385/1‑151 (MQ=255)                                                                      tATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCtgttgt                                                                          >  1:391519/1‑151 (MQ=255)                                                                       aTGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCtgttgt                                                                          <  1:1082652/150‑1 (MQ=255)                                                                       aTGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTgttgtt                                                                         <  1:767714/151‑1 (MQ=255)                                                                           cggcaacggcaaATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACgg                                                                     <  1:98337/151‑1 (MQ=255)                                                                             gcaacggcaaATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCa                                                                   <  1:60564/151‑1 (MQ=255)                                                                                acggcaaATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAAc                                                                 >  2:987612/1‑150 (MQ=255)                                                                                   gcaaATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAggg                                                             <  2:468948/151‑1 (MQ=255)                                                                                         ttATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGtt                                                       >  1:317911/1‑151 (MQ=255)                                                                                              cAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCt                                                  <  1:221221/151‑1 (MQ=255)                                                                                              cAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCt                                                  <  2:1052510/151‑1 (MQ=255)                                                                                                gTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCttt                                                <  1:379530/151‑1 (MQ=255)                                                                                                      gcgcACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGtctc                                          >  2:1016574/1‑151 (MQ=255)                                                                                                            aGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCttt                                    >  2:510335/1‑151 (MQ=255)                                                                                                             ggTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCtttt                                   >  1:572570/1‑151 (MQ=255)                                                                                                                         cAAATCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTAcc                        <  1:56171/150‑1 (MQ=255)                                                                                                                          aaaTCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCt                      <  1:67416/151‑1 (MQ=255)                                                                                                                          aaaTCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCt                      >  1:435070/1‑151 (MQ=255)                                                                                                                          aaaTCACTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCt                      <  2:669130/151‑1 (MQ=255)                                                                                                                               aCTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATcc                  >  2:922932/1‑150 (MQ=255)                                                                                                                                cTTCTCTCAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCAt                >  2:504488/1‑151 (MQ=255)                                                                                                                                  tctctcAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCAt                <  1:611033/149‑1 (MQ=255)                                                                                                                                  tctctcAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATa               <  1:1067244/150‑1 (MQ=255)                                                                                                                                       cAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGcaca         <  1:849560/151‑1 (MQ=255)                                                                                                                                       cAGTGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGcac          >  2:388350/1‑150 (MQ=255)                                                                                                                                          tGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAAt       >  1:478779/1‑150 (MQ=255)                                                                                                                                          tGAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAAt       >  2:250136/1‑150 (MQ=255)                                                                                                                                           gAAGACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATaa     <  1:66277/151‑1 (MQ=255)                                                                                                                                              gACCACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGc  >  2:1001333/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GAACTGAAAAACGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGC  >  NC_000913/1735047‑1735338

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