Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	2,066,318	A→G	L511L (CTA→CTG)	unknown →	Peptidoglycan lipid II flippase MurJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	2,066,318	0	A	G	100.0%	234.0 / NA	78	L511L (CTA→CTG)	unknown	Peptidoglycan lipid II flippase MurJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (40/38);  total (40/38)

TGTCGGCGACACTCAACGCGTTTCTGCTTTATCGCGGCTTGGCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTATTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATACCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATGGCGTGCATCTGGGGCATCAAGCTATTTTGCGTCATTTGCGGGAAAAA  >  1/2066071‑2066563                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      tGTCGGCGACACTCAACGCGTTTCTGCTTTATCGCGGCTTGGCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTg                                                                                                                                                                                                                                                    <  1:661059/251‑1 (MQ=255)      gCGACACTCAACGCGTTTCTGCTTTATCGCGGCTTGGCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGAt                                                                                                                                                                                                                                               >  1:351682/1‑251 (MQ=255)       cGACACTCAACGCGTTTCTGCTTTATCGCGGCTTGGCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGTCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGAt                                                                                                                                                                                                                                               <  1:802917/250‑1 (MQ=255)       cGACACTCAACGCGTTTCTGCTTTATCGCGGCTTGGCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGAt                                                                                                                                                                                                                                               <  1:283079/250‑1 (MQ=255)        gACACTCAACGCGTTTTTGCTTTATCGCGGCTTGGCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGAtt                                                                                                                                                                                                                                              <  1:641653/250‑1 (MQ=255)            cTCAACGCGTTTCTGCTTTATCGCGGCTTGGCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCaa                                                                                                                                                                                                                                         >  1:140699/1‑251 (MQ=255)                aCGCGTTTCTGCTTTATCGCGGCTTGGCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGCGGCATGTTGGTGCTGTTGGGGATTCGCAAACGt                                                                                                                                                                                                                                     >  1:217073/1‑251 (MQ=255)                     tttCTGCTTTATCGCGGCTTGGCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGAtt                                                                                                                                                                                                                                 <  1:787141/250‑1 (MQ=255)                                   cgGCTTGGCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACg                                                                                                                                                                                                                            <  1:410437/241‑1 (MQ=255)                                   cgGCTTGGCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCAc                                                                                                                                                                                                                   <  1:633658/250‑1 (MQ=255)                                     gCTTGGCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTc                                                                                                                                                                                                                 <  1:898344/250‑1 (MQ=255)                                         ggCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGt                                                                                                                                                                                                            >  1:729233/1‑251 (MQ=255)                                                        gtATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGAttt                                                                                                                                                                                                                                >  1:208005/1‑216 (MQ=255)                                                        gtATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACg                                                                                                                                                                                                                            >  1:417275/1‑220 (MQ=255)                                                        gtATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACg                                                                                                                                                                                                                            <  1:762989/220‑1 (MQ=255)                                                        gtATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTtgt                                                                                                                                                                                             >  1:666709/1‑251 (MQ=255)                                                        gtATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTtgt                                                                                                                                                                                             >  1:347687/1‑251 (MQ=255)                                                                            gTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGc                                                                                                                                                                            >  1:491233/1‑248 (MQ=255)                                                                                            gAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTAttt                                                                                                                                                         >  1:674508/1‑251 (MQ=255)                                                                                                 ttttGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTAtttttt                                                                                                                                                      <  1:772378/249‑1 (MQ=255)                                                                                                 ttttGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTc                                                                                                                                                     <  1:716375/250‑1 (MQ=255)                                                                                                     gtttgCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGAtt                                                                                                                                                                                 >  1:151610/1‑218 (MQ=255)                                                                                                        tgCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAAt                                                                                                                                             <  1:141568/251‑1 (MQ=255)                                                                                                             ccgTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTAt                                                                                                                                         <  1:51790/250‑1 (MQ=255)                                                                                                               gTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCAc                                                                                                                                                                                                                   >  1:12537/1‑174 (MQ=255)                                                                                                                 catcatGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAac                                                                                                                                    <  1:177148/251‑1 (MQ=255)                                                                                                                   tcatGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAaca                                                                                                                                   >  1:567521/1‑250 (MQ=255)                                                                                                                   tcatGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGAGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAaca                                                                                                                                   >  1:120178/1‑250 (MQ=255)                                                                                                                            tgttgGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGc                                                                                                                         <  1:1881/251‑1 (MQ=255)                                                                                                                                ggtttggtACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTc                                                                                                                                                                                                                 <  1:846961/159‑1 (MQ=255)                                                                                                                                    tggtACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGc                                                                                                                                                                            <  1:571673/192‑1 (MQ=255)                                                                                                                                    tggtACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAac                                                                                                                                    >  1:144860/1‑232 (MQ=255)                                                                                                                                    tggtACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGAtt                                                                                                                 >  1:290199/1‑251 (MQ=255)                                                                                                                                     ggtACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGAtt                                                                                                                 >  1:57865/1‑250 (MQ=255)                                                                                                                                     ggtACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGAtt                                                                                                                 >  1:410918/1‑250 (MQ=255)                                                                                                                                      gtACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATg                                                                                                                                        <  1:528840/226‑1 (MQ=255)                                                                                                                                       tACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCg                                                                                                               <  1:2855/250‑1 (MQ=255)                                                                                                                                         caacaGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAgg                                                                                                            <  1:416757/251‑1 (MQ=255)                                                                                                                                              gCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCAc                                                                                                                                                                                                                   <  1:563834/143‑1 (MQ=255)                                                                                                                                              gCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGc                                                                                                                         <  1:746613/233‑1 (MQ=255)                                                                                                                                                    gCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAAtt                                                                                                 >  1:474765/1‑251 (MQ=255)                                                                                                                                                                    cGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCa                                                                                 >  1:152167/1‑251 (MQ=255)                                                                                                                                                                     gATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACATATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCa                                                                                 >  1:225351/1‑250 (MQ=255)                                                                                                                                                                      atatGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACg                                                                               >  1:123472/1‑251 (MQ=255)                                                                                                                                                                               ttCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGTTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACTGAAATTTATGGAACAAATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGt                                                                      >  1:31601/1‑251 (MQ=255)                                                                                                                                                                               ttCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGt                                                                      <  1:94178/251‑1 (MQ=255)                                                                                                                                                                                   cACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAAc                                                                  <  1:46243/251‑1 (MQ=255)                                                                                                                                                                                    aCCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAAcc                                                                 <  1:67323/251‑1 (MQ=255)                                                                                                                                                                                     cccGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGAt                                                                                                                                                                                  <  1:337934/137‑1 (MQ=255)                                                                                                                                                                                     cccGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTa                                                                                                                                                            >  1:500803/1‑159 (MQ=255)                                                                                                                                                                                     cccGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCa                                                                >  1:76263/1‑251 (MQ=255)                                                                                                                                                                                     cccGTGTACATTAGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCa                                                                >  1:719827/1‑251 (MQ=255)                                                                                                                                                                                       cGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGc                                                                                                                                                                            <  1:885827/141‑1 (MQ=255)                                                                                                                                                                                       cGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATgg                                                                                                                                       >  1:647289/1‑178 (MQ=255)                                                                                                                                                                                                 ggCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTgggg                                                                                                                                                                                                                                                 <  1:405182/62‑1 (MQ=255)                                                                                                                                                                                                 ggCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTg                                                    >  1:513700/1‑251 (MQ=255)                                                                                                                                                                                                 ggCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTg                                                    >  1:568719/1‑251 (MQ=255)                                                                                                                                                                                                   cTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGAt                                                  <  1:126848/251‑1 (MQ=255)                                                                                                                                                                                                      aCGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATgg                                                >  1:571103/1‑250 (MQ=255)                                                                                                                                                                                                      aCGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATGGc                                               <  1:246896/251‑1 (MQ=255)                                                                                                                                                                                                            cTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATACAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATGGCGTGCAt                                         >  1:429350/1‑251 (MQ=255)                                                                                                                                                                                                               aTCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATGGCGTGCATCTg                                      <  1:106191/251‑1 (MQ=255)                                                                                                                                                                                                                 cGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAg                                                                                                                                                              >  1:238867/1‑129 (MQ=255)                                                                                                                                                                                                                   tgtTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAAttt                                                                                                <  1:735557/189‑1 (MQ=255)                                                                                                                                                                                                                               tGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATGGCGTGCATCTGGGGCATCAAGCTAttt                      <  1:99390/251‑1 (MQ=255)                                                                                                                                                                                                                               tGGTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATGGCGTGCATCTGGGGCATCAAGCTAttt                      <  1:626992/251‑1 (MQ=255)                                                                                                                                                                                                                                ggTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTataatgc                                                                                                                                                                   <  1:188145/109‑3 (MQ=255)                                                                                                                                                                                                                                ggTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGt                                                                      >  1:607844/1‑202 (MQ=255)                                                                                                                                                                                                                                ggTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGt                                                                      <  1:280713/202‑1 (MQ=255)                                                                                                                                                                                                                                 gTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATGGCGTGCGTCTGGGGCATCAAGCTATTTTg                    >  1:597775/1‑251 (MQ=255)                                                                                                                                                                                                                                 gTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATGGCGTGCATCTGGGGCATCAAGCTATTTTg                    >  1:416020/1‑251 (MQ=255)                                                                                                                                                                                                                                 gTTTACGGTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATGGCGTGCATCTGGGGCATCAAGCTATTTTg                    >  1:881925/1‑251 (MQ=255)                                                                                                                                                                                                                                        gTGGCATGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAGGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATGGCGTGCATCTGGGGCATCAAGCTATTTTGCGTCAtt             >  1:818306/1‑251 (MQ=255)                                                                                                                                                                                                                                              tGTTGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTa                                                          <  1:737219/200‑1 (MQ=255)                                                                                                                                                                                                                                                ttGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGc                                                                                                                                                                            >  1:108141/1‑84 (MQ=255)                                                                                                                                                                                                                                                ttGGTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATGGCGTGCATCTGGGGCATCAAGCTATTTTGCGTCATTTGCGGGaa     >  1:636804/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATgg                                                                                                                                       <  1:732576/118‑1 (MQ=255)                                                                                                                                                                                                                                                   gTGCTGTTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATGCCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATGGCGTGCATCTGGGGCATCAAGCTATTTTGCGTCATTTGCGGGaaaaa  >  1:172248/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      TGTCGGCGACACTCAACGCGTTTCTGCTTTATCGCGGCTTGGCGAAAGCCGATGTGTATCATTTCAGCCGACAAAGTGCGGTCTTTTTTCTGAAAGTTTTGTTTGCCGCCGTCATCATGGGATTGTTGGTTTGGTACAACAGCCCAAGCCTGACGGAATGGAACGATATGCGCTTTCTCACCCGTGTACATTGGCTGACGTGGCTCATCGTGTTGGCAGCGATGGTTTACGGTGGCATGTTGGTGCTATTGGGGATTCGCAAACGTGATTTGACGGGCGCGCACTCCGCGTGAAATTAGTTTTTGTGGCGCATGGATTTATGCGTATAATACCGAAGTATTTTTTCACCGAAATTTATGGAACACATTGAATGCAGCTGATTCGAGGCCTGGGAAATTTACCGCAAAATTTGCACGCTTGCGCGTTAACCATCGGTAATTTTGATGGCGTGCATCTGGGGCATCAAGCTATTTTGCGTCATTTGCGGGAAAAA  >  1/2066071‑2066563

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