Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	493,778	C→T	D87D (GAC→GAT)	unknown →	LPS‑core synthesis glycosyltransferase PM0512

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	493,778	0	C	T	100.0%	125.1 / NA	47	D87D (GAC→GAT)	unknown	LPS‑core synthesis glycosyltransferase PM0512
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/29);  total (18/29)

TCCTTATATATTAAGGAAAATCCACAATATTTAAGAGAATGTTTGGAAAGTTTGGTGGCACAAACCTGCCCTGCCGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGACTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACTTCTGCGCAAGACATCGTTACATTTACCCAAAAACGCTGCCCGTTTAATCACATGACGGTGGCGTAT  >  1/493545‑494006                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                     tCCTTATATATTAAGGAAAATCCACAATATTTAAGAGAATGTTTGGAAAGTTTGGTGGCACAAACCTGCCCTGCCGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGa                                                                                                                                                                                                                     <  1:348390/251‑1 (MQ=255) tCCTTATATATTAAGGAAAATCCACAATATTTAAGAGAATGTTTGGAAAGTTTGGTGGCACAAACCTGCCCTGCCGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGa                                                                                                                                                                                                                     >  1:261595/1‑251 (MQ=255) tCCTTATATATTAAGGAAAATCCACAATATTTAAGAGAATGTTTGGAAAGTTTGGTGGCACAAACCTGCCCTGCCGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGa                                                                                                                                                                                                                     >  1:76499/1‑251 (MQ=255)  ccTTATATATTAAGGAAAATCCACAATATTTAAGAGAATGTTTGGAAAGTTTGGTGGCACAAACCTGCCCTGCCGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAGTTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATgg                                                                                                                                                                                                                      >  1:452830/1‑249 (MQ=255)    ttATATATTAAGGAAAATCCACAATATTTAAGAGAATGTTTGGAAAGTTTGGTGGCACAAACCTGCCCTGCCGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATAc                                                                                                                                                                                                                  <  1:35720/251‑1 (MQ=255)            tAAGGAAAATCCACAATATTTAAGAGAATGTTTGGAAAGTTTGGTGGCACAAACCTGCCCTGCCGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGata                                                                                                                                                                                                          <  1:354172/251‑1 (MQ=255)                           atatTTAAGAGAATGTTTGGAAAGTTTGGTGGCACAAACCTGCCCTGCCGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGaa                                                                                                                                                                                            >  1:8081/1‑250 (MQ=255)                           atatTTAAGAGAATGTTTGGAAAGTTTGGTGGCACAAACCTGCCCTGCCGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGaa                                                                                                                                                                                            >  1:348874/1‑250 (MQ=255)                           atatTTAAGAGAATGTTTGGAAAGTTTGGTGGCACAAACCTGCCCTGCCGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGaa                                                                                                                                                                                            >  1:148716/1‑250 (MQ=255)                                                                     ccctgccGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGAt                                                                                                                                                  <  1:874509/250‑1 (MQ=255)                                                                      cctggcGATGAAGTCGTTCTGGTGGTTGATGGGGCGGTGGCGGCTTAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATAc                                                                                                                                                <  1:882133/251‑1 (MQ=255)                                                                          ccGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACtat                                                                                                                                             <  1:76633/250‑1 (MQ=255)                                                                           cGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAg                                                                                                                                           <  1:142454/251‑1 (MQ=255)                                                                           cGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAg                                                                                                                                           <  1:560364/251‑1 (MQ=255)                                                                           cGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAg                                                                                                                                           <  1:41617/251‑1 (MQ=255)                                                                                aaGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTc                                                                                                                                      <  1:458373/251‑1 (MQ=255)                                                                                  gTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACAt                                                                                                                                                        >  1:513455/1‑231 (MQ=255)                                                                                  gTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCgg                                                                                                                                    >  1:92803/1‑251 (MQ=255)                                                                                                                        ggAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGaac                                                                                                                                                             <  1:697547/188‑1 (MQ=255)                                                                                                                          aaGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACAt                                                                                             <  1:48796/250‑1 (MQ=255)                                                                                                                          aaGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACAt                                                                                             <  1:567241/250‑1 (MQ=255)                                                                                                                               gtggtTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTgg                                                                                        <  1:59684/250‑1 (MQ=255)                                                                                                                                tggtTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCa                                                                                      <  1:863917/251‑1 (MQ=255)                                                                                                                                           tttGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAAtg                                                                            <  1:899357/250‑1 (MQ=255)                                                                                                                                                    aaGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACtt                                                                   >  1:811813/1‑250 (MQ=255)                                                                                                                                                           cGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGAt                                                                                                                                                  <  1:147477/164‑1 (MQ=255)                                                                                                                                                            gTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAg                                                                                                                                           >  1:218565/1‑170 (MQ=255)                                                                                                                                                            gTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACTTCTgcgc                                                             <  1:313521/248‑1 (MQ=255)                                                                                                                                                                   ttGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACTTCTGCGCAAGACATCGt                                                   <  1:709777/251‑1 (MQ=255)                                                                                                                                                                      gTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAAtg                                                                            >  1:725779/1‑223 (MQ=255)                                                                                                                                                                         aaaTTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACTTCTGCGCAAGACATCGTTACAt                                              <  1:491075/250‑1 (MQ=255)                                                                                                                                                                            ttACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACTTCTGCGCAAGACATCGTTACATTTa                                           >  1:38878/1‑250 (MQ=255)                                                                                                                                                                                  aaCAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAA‑‑ATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGAGATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACTTCTGCGCAAGACATCGTTACATTTACCCAAAAAc                                  <  1:442503/251‑1 (MQ=255)                                                                                                                                                                                     aaTTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACTTCTGCGCAAGACATCGTTa                                                 <  1:633007/235‑1 (MQ=255)                                                                                                                                                                                      aTTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAg                                                                                                                                                                             <  1:784738/110‑1 (MQ=255)                                                                                                                                                                                      aTTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAg                                                                                                                                           <  1:236638/144‑1 (MQ=255)                                                                                                                                                                                      aTTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCat                                                                                     <  1:154316/198‑1 (MQ=255)                                                                                                                                                                                                     aaGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACTTCTGCGCAAGACATCGTTACATTTACCCAAAAACGCTGCCCGTTTAATcac                 <  1:109978/251‑1 (MQ=255)                                                                                                                                                                                                            cAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACTTCTGCGCAAGACATCGTTAc                                                <  1:681186/213‑1 (MQ=255)                                                                                                                                                                                                            cAACGAAGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACTTCTGCGCAAGACATCGTTACATTTACCCAAAAACGCTGCCCGTTTAATCACATGAc            >  1:137573/1‑249 (MQ=255)                                                                                                                                                                                                                 aaGGTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACTTCTGCGCAAGACATCGTTACATTTACCCAAAAACGCTGCCCGTTTAATCACATGACGGTGGCg     <  1:642873/251‑1 (MQ=255)                                                                                                                                                                                                                   ggTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATcc                                                                                                                                                      >  1:812822/1‑104 (MQ=255)                                                                                                                                                                                                                   ggTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATAc                                                                                                                                                >  1:137095/1‑110 (MQ=255)                                                                                                                                                                                                                   ggTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTgg                                                                                        <  1:465417/166‑1 (MQ=255)                                                                                                                                                                                                                    gTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGCAATGTGCCGACTTCTGCGCAAAACATCGTTACATTTACCCAAAAACGCTGCCCGTTTAATCACATGACGGTGGCGTAt  >  1:851780/1‑251 (MQ=255)                                                                                                                                                                                                                    gTTTAAAATATTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACTTCTGCGCAAGACATCGTTACATTTACCCAAAAACGCTGCCCGTTTAATCACATGACGGt         >  1:430550/1‑244 (MQ=255)                                                                                                                                                                                                                           atatTGCTCGCGCGATTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGt                                                                                          >  1:770393/1‑156 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                     TCCTTATATATTAAGGAAAATCCACAATATTTAAGAGAATGTTTGGAAAGTTTGGTGGCACAAACCTGCCCTGCCGATGAAGTCGTACTGGTGTTTGATGGCGCGGTGACGCCTGAATTGGAAGCGGTGGTTGCCGAGTTTGAATCGAAGCTGCCGTTTAATTTGGTGAAATTACCGAACAATTTAGGTTTAGGCAAAGCTCTCAACGAAGGTTTAAAATATTGCTCGCGCGACTGGGTTCTGCGCATGGATACCGATGATATTTGCGTGCCGGAACGTTTTGCCAAACAGGTGGCGTTTATTGAACAACATCCTGATACTATAGTGTTCGGTGGGCAAATTGCCGAATTCGGTGAAAATATCCATGACATAGTGGCATATCGAAATGTGCCGACTTCTGCGCAAGACATCGTTACATTTACCCAAAAACGCTGCCCGTTTAATCACATGACGGTGGCGTAT  >  1/493545‑494006

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