Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	285,274	G→A	T3T (ACC→ACT)	unknown ←	UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase (EC 6.3.2.‑)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	285,274	0	G	A	100.0%	144.7 / NA	54	T3T (ACC→ACT)	unknown	UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase (EC 6.3.2.‑)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/30);  total (24/30)

TAACATATACTCGATGCAGGGATTGCCGCGTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTGGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGTAAATTCATTAATGAAGTTAAGGAATAAAAGATGAAAACATTAGGTCAATTTATCGTAGAAAAACAAGCGGAATACCCACATGCTAAAGGAGAATTGAGCGG  >  1/285028‑285510                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             tAACATATACTCGATGCAGGGATTGCCGCGTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTg                                                                                                                                                                                                                                          <  1:159421/251‑1 (MQ=255) tAACATATACTCGATGCAGGGATTGCCGCGTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTg                                                                                                                                                                                                                                          <  1:510524/251‑1 (MQ=255)    cATATACTCGATGCAGGGATTGCCGCGTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCa                                                                                                                                                                                                                                       >  1:630032/1‑251 (MQ=255)    cATATACTCGATGCAGGGATTGCCGCGTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCa                                                                                                                                                                                                                                       >  1:224294/1‑251 (MQ=255)    cATATACTCGATGCAGGGATTGCCGCGTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCa                                                                                                                                                                                                                                       >  1:283359/1‑251 (MQ=255)              aTGCAGGGATTGCCGCGTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTc                                                                                                                                                                                                                              >  1:355016/1‑250 (MQ=255)                      aTTGCCGCGTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTtatat                                                                                                                                                                                                                     <  1:805079/251‑1 (MQ=255)                      aTTGCCGCGTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTtatat                                                                                                                                                                                                                     <  1:679226/251‑1 (MQ=255)                      aTTGCCGCGTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTtatat                                                                                                                                                                                                                     <  1:406643/251‑1 (MQ=255)                      aTTGCCGCGTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTCGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTtatat                                                                                                                                                                                                                     <  1:880064/251‑1 (MQ=255)                             cgTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGtt                                                                                                                                                                                                              >  1:9075/1‑251 (MQ=255)                             cgTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGtt                                                                                                                                                                                                              <  1:273958/251‑1 (MQ=255)                                   aTGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTc                                                                                                                                                                                                        <  1:314392/251‑1 (MQ=255)                                      gCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTCAttt                                                                                                                                                                                                         >  1:638445/1‑247 (MQ=255)                                                   tGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAAt                                                                                                                                                                                        >  1:179671/1‑251 (MQ=255)                                                    gACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAAtg                                                                                                                                                                                       >  1:24066/1‑251 (MQ=255)                                                     aCCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAAtg                                                                                                                                                                                       <  1:407377/250‑1 (MQ=255)                                                               gcgcTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAAt                                                                                                                                                                             <  1:496154/250‑1 (MQ=255)                                                                      ttGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGaa                                                                                                                                                                      <  1:477421/250‑1 (MQ=255)                                                                        gTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATcac                                                                                                                                                                          <  1:394987/244‑1 (MQ=255)                                                                          aaCTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAgg                                                                                                                                                                 <  1:398658/251‑1 (MQ=255)                                                                           aCTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGa                                                                                                                                                                <  1:721545/251‑1 (MQ=255)                                                                                   ggctCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTtat                                                                                                                                                                                                                       <  1:365712/188‑1 (MQ=255)                                                                                   ggctCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAg                                                                                                                                                                  <  1:7966/241‑1 (MQ=255)                                                                                      tCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAAt                                                                                                                                                     <  1:672821/251‑1 (MQ=255)                                                                                      tCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAAt                                                                                                                                                     <  1:571383/251‑1 (MQ=255)                                                                                      tCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAAt                                                                                                                                                     <  1:800859/251‑1 (MQ=255)                                                                                                gTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTtatat                                                                                                                                                                                                                     >  1:597866/1‑177 (MQ=255)                                                                                                gTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAAtc                                                                                                                                                             >  1:550988/1‑233 (MQ=255)                                                                                                      tttCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGc                                                                                                                                     <  1:42093/251‑1 (MQ=255)                                                                                                                aTTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTataat                                                                                                                            <  1:633075/250‑1 (MQ=255)                                                                                                                           gAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGa                                                                                                                                                                >  1:656954/1‑203 (MQ=255)                                                                                                                                    tCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGt                                                                                                       <  1:590520/251‑1 (MQ=255)                                                                                                                                          gtggGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCgg                                                                                                                                                                                                                >  1:242071/1‑140 (MQ=255)                                                                                                                                          gtggGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGTAAAtt                                                                                                  <  1:503004/250‑1 (MQ=255)                                                                                                                                            gggTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTcc                                                                                                                                                                                                                             >  1:855238/1‑125 (MQ=255)                                                                                                                                                  aCGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTc                                                                                                                 >  1:668981/1‑227 (MQ=255)                                                                                                                                                        gtgtCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTtatat                                                                                                                                                                                                                     <  1:823452/121‑1 (MQ=255)                                                                                                                                                                  cTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGTAAATTCATTAATGAAGTTAAGGAATAAAAg                                                                         >  1:260524/1‑251 (MQ=255)                                                                                                                                                                    gTTACATGATAACCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGTAAATTCATTAATGAAGTTAAGGAATAAAAGAt                                                                       >  1:454720/1‑251 (MQ=255)                                                                                                                                                                     ttACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAAGTGAGTAAATTCATTAATGAAGTTAAGGAATAAAAGATg                                                                      <  1:351059/251‑1 (MQ=255)                                                                                                                                                                          tGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGTAAATTCATTAATGAAGTTAAGGAATAAAAGATGaaaa                                                                  >  1:81608/1‑250 (MQ=255)                                                                                                                                                                                cccATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAAtg                                                                                                                                                                                       >  1:194971/1‑127 (MQ=255)                                                                                                                                                                                             cGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGTAAATTCATTAATGAAGTTAAGGAATAAAAGATGAAAACATTAGGTCAATTTATCGt                                               >  1:364458/1‑250 (MQ=255)                                                                                                                                                                                                   tcGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTg                                                                                                                                                                                                                                  >  1:483887/1‑65 (MQ=255)                                                                                                                                                                                                      cAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAg                                                                                                                    >  1:39187/1‑172 (MQ=255)                                                                                                                                                                                                             gccCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTtatat                                                                                                                                                                                                                     <  1:712409/68‑1 (MQ=255)                                                                                                                                                                                                               ccATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGTAAATTCATTAATGAAGTTAAGGAATAAAAGATGAAAACATTAGGTCAATTTATCGTAGAAAAACAAGCGGAATAc                            <  1:177360/251‑1 (MQ=255)                                                                                                                                                                                                                  tAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGTAAATTCATTAATGAAGTTAAGGAATAAAAGATGAAAACATTAGGTCAATTTATCGTAGAAAAACAAGCGGAATACCca                         <  1:400089/251‑1 (MQ=255)                                                                                                                                                                                                                             cAAATGCCGAGAATATGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGTAAATTCATTAATGAAGTTAAGGAATAAAAGATGAAAACATTAGGTCAATTTATCGTAGAAAAACAAGCGGAATACCCACATGCTAAAgg              <  1:173832/251‑1 (MQ=255)                                                                                                                                                                                                                                         atatGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAAt                                                                                                                                                     >  1:739590/1‑104 (MQ=255)                                                                                                                                                                                                                                         atatGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGTAAATTCATTAATGAAGTTAAGGAATAAAAGATGAAAACATTAGGTCAATTTATCGTAGAAAAACAAGCGGAATACCCACATGCTAAAGGAGAATTGAGCgg  >  1:145558/1‑251 (MQ=255)                                                                                                                                                                                                                                          tatGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGTAAAtt                                                                                                  <  1:773579/154‑1 (MQ=255)                                                                                                                                                                                                                                           atGAAGATGTTTAGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGTAAATTCATTAATGAAGTTAAGGAATAAAAGATGAAAACATTAGGTCAATTTATCGTAGAAAAACAAGCGGAATACCcac                        >  1:290383/1‑227 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             TAACATATACTCGATGCAGGGATTGCCGCGTTTCATGGCATTTCCCACGATGACCATATCCGGCGCTGGTTGTAACTGGCTTGGCTCAAAGTGCGGTATGATTTCAATGTGATTTTCGGCAAGAAAAGTGCTCATTGGTGGGTAAACGTTGGTGTCGGAGCCTGTTACATGATAGCCCATTTGCTTGGCGATCATCGCAATGCCGCCCATAAATGTGCCGCAAATGCCGAGAATATGAAGATGTTTGGTTGTCATGTTGAGTCCTTTTATATACCGGTTTATTTCTGTTTCCAATAAAAATGTGTTGTGAATCACAGAAATAGGAATCTCTTTAATCATTCAAAATTGAAGCTGTTATAATAAAATAAGTTCAGAATTGAGTAAATTCATTAATGAAGTTAAGGAATAAAAGATGAAAACATTAGGTCAATTTATCGTAGAAAAACAAGCGGAATACCCACATGCTAAAGGAGAATTGAGCGG  >  1/285028‑285510

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