Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,588,108	A→G	Y580H (TAC→CAC)	unknown ←	Lon protease homolog YcbZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,588,108	0	A	G	100.0%	183.6 / NA	52	Y580H (TAC→CAC)	unknown	Lon protease homolog YcbZ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/23);  total (29/23)

CAATAAATCTTCGTAAGTGTAACTGTCTTTGATATTTGGAGTACAAGTGTTCATTATTAATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTAAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGATAAATAAAAAATTCCTGATTTTTAACCGCACTTACCACCGCTTCGGACAAGCTTAGCTGATTAAGCACCGCACTCGGAATCATCACGCCTTGTTTGCCC  >  1/1587865‑1588352                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     cAATAAATCTTCGTAAGTGTAACTGTCTTTGATATTTGGAGTACAAGTGTTCATTATTAATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGAttgt                                                                                                                                                                                                                                               >  1:59735/1‑251 (MQ=255)               aaGTGTAACTGTCTTTGATATTTGGAGTACAAGTGTTCATTATTAATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGAttgt                                                                                                                                                                                                                                               >  1:708879/1‑237 (MQ=255)                aGTGTAACTGTCTTTGATATTTGGAGTACAAGTGTTCATTATTAATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTg                                                                                                                                                                                                                                 >  1:818537/1‑250 (MQ=255)                         gTCTTTGATATTTGGAGTACAAGTGTTCATTATTAATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTGACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACg                                                                                                                                                                                                                        >  1:816146/1‑250 (MQ=255)                         gTCTTTGATATTTGGAGTACAAGTGTTCATTATTAATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACg                                                                                                                                                                                                                        >  1:327803/1‑250 (MQ=255)                                  atTTGGAGTACAAGTGTTCATTATTAATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGCTTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCtgattgat                                                                                                                                                                                                               >  1:434961/1‑250 (MQ=255)                                         gTACAAGTGTTCATTATTAATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGc                                                                                                                                                                                                                                <  1:623413/226‑1 (MQ=255)                                                  ttCATTATTAATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGc                                                                                                                                                                                               <  1:811449/250‑1 (MQ=255)                                                   tCATTATTAATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAg                                                                                                                                                                                             <  1:380666/251‑1 (MQ=255)                                                   tCATTATTAATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAg                                                                                                                                                                                             <  1:851953/251‑1 (MQ=255)                                                      ttatttatCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAgtgt                                                                                                                                                                                          >  1:205846/1‑251 (MQ=255)                                                      ttattaATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTgaga                                                                                                                                                                                                                                                   <  1:205506/194‑1 (MQ=255)                                                      ttattaATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACa                                                                                                                                                                                                    >  1:586876/1‑241 (MQ=255)                                                      ttattaATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAgtgt                                                                                                                                                                                          <  1:366456/251‑1 (MQ=255)                                                      ttattaATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAgtgt                                                                                                                                                                                          >  1:394882/1‑251 (MQ=255)                                                         ttaATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGtgtg                                                                                                                                                                                         >  1:12988/1‑249 (MQ=255)                                                             tcatcCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGaaa                                                                                                                                                                                   <  1:710018/251‑1 (MQ=255)                                                                  cTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAg                                                                                                                                                                                             <  1:836564/236‑1 (MQ=255)                                                                            tGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCCAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATtc                                                                                                                                                                    <  1:430028/251‑1 (MQ=255)                                                                                    cGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATTCGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGc                                                                                                                                                              <  1:256482/249‑1 (MQ=255)                                                                                                        gTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGtaat                                                                                                                                         <  1:881504/250‑1 (MQ=255)                                                                                                           ggATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGAc                                                                                                                                                                                                     <  1:49098/187‑1 (MQ=255)                                                                                                           ggATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATaaaa                                                                                                                                     >  1:607768/1‑251 (MQ=255)                                                                                                            gATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAAt                                                                                                                                    >  1:302365/1‑251 (MQ=255)                                                                                                                 tAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGc                                                                                                                                                                                               >  1:495864/1‑187 (MQ=255)                                                                                                                  aaCAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTcaca                                                                                                                              >  1:396169/1‑251 (MQ=255)                                                                                                                           ttCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGCCAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGAt                                                                                                                      >  1:437725/1‑250 (MQ=255)                                                                                                                           ttCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATc                                                                                                                     <  1:526748/251‑1 (MQ=255)                                                                                                                           ttCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATc                                                                                                                     <  1:618095/251‑1 (MQ=255)                                                                                                                                  aTTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATgg                                                                                                                                                                                                >  1:841164/1‑169 (MQ=255)                                                                                                                                  aTTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGtt                                                                                                              >  1:867221/1‑251 (MQ=255)                                                                                                                                  aTTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGtt                                                                                                              >  1:849769/1‑251 (MQ=255)                                                                                                                                          tAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGc                                                                                                                                                                                               <  1:611222/162‑1 (MQ=255)                                                                                                                                           aGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATgg                                                                                                                                                                                                <  1:489417/160‑1 (MQ=255)                                                                                                                                           aGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTc                                                                                                                                                            >  1:877993/1‑196 (MQ=255)                                                                                                                                                 cacTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGgtaaata                                                                                               >  1:828637/1‑251 (MQ=255)                                                                                                                                                 cacTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGataaata                                                                                               >  1:64167/1‑251 (MQ=255)                                                                                                                                                 cacTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGataaata                                                                                               >  1:373020/1‑251 (MQ=255)                                                                                                                                                         cGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACg                                                                                                                                                                                                         >  1:845759/1‑137 (MQ=255)                                                                                                                                                         cGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGATAAATAAAAAAtt                                                                                        >  1:588025/1‑250 (MQ=255)                                                                                                                                                              aaCTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGATAAATAAAAAATTCCTGAt                                                                                  >  1:694281/1‑251 (MQ=255)                                                                                                                                                              aaCTTGTCCGATCAGTTATTATTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGATAAATAAAAAATTCCTGAt                                                                                  >  1:143537/1‑251 (MQ=255)                                                                                                                                                                      cGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGATAAATAAAAAATTCCTGATTTTTAAcc                                                                          <  1:23856/251‑1 (MQ=255)                                                                                                                                                                       gATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCCCCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGATAAATAAAAAATTCCTGATTTTTAACCg                                                                         >  1:423923/1‑251 (MQ=255)                                                                                                                                                                               aTTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGATAAATAAAAAATTCCTGATTTTTAACCGCACTTacc                                                                 <  1:698172/251‑1 (MQ=255)                                                                                                                                                                                          gAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATgg                                                                                                                                                                                                >  1:590728/1‑113 (MQ=255)                                                                                                                                                                                          gAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGaaa                                                                                                                                                                                   >  1:537012/1‑126 (MQ=255)                                                                                                                                                                                                   aTTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAAtc                                                                                                                                   <  1:734164/165‑1 (MQ=255)                                                                                                                                                                                                            ttttttCTGCCGAAGAAT‑AATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGATAAATAAAAAATTCCTGATTTTTAACCGCACTTACCACCGCTTCGGACAAGCTTAGCTGATTAAGc                                   <  1:186633/251‑1 (MQ=255)                                                                                                                                                                                                               cttCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGATAAATAAAAAATTCCTGATTTTTAACCGCACTTACCACCGCTTCGGACAAGCTTAGCTGATTAAGCAc                                 <  1:92144/250‑1 (MQ=255)                                                                                                                                                                                                                            aTAAATCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGATAAATAAAAAATTCCTGATTTTTAACCGCAc                                                                      <  1:357737/201‑1 (MQ=255)                                                                                                                                                                                                                                 tCCCAAAGACCGTGCTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCg                                                                                                                                                                                <  1:662445/90‑1 (MQ=255)                                                                                                                                                                                                                                               cTTGTGAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGATAAATAAAAAATTCCTGATTTTTAACCGCACTTACCACCGCTTCGGACAAGCTTAGCTGATTAAGCACCGCACTCGGAATCATCACGCCTTGTTTGccc  <  1:675643/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     CAATAAATCTTCGTAAGTGTAACTGTCTTTGATATTTGGAGTACAAGTGTTCATTATTAATCATCCTATAAAATTTGTTTAAGCGATTCAAAATTGAATGCGGGTGGGATTCTAACAGAGCATTCCGGCATTTAGAATAGCTAACACTTGTACGCCTAACTTGTCCGATCAGTTATTTTTTCATGGAAATAGCGATTTAACGATTTTTTCTGCCGAAGAATAAATCCCAAAGACCGTGCTTGTAAGATTGTTGCTCGGAATATTGCTCAATACGCTGATTGATCAAACGGGACAATGGCAGTGTGCTGAAATCGTACTTTTTATTCTCCTGCTCCACCAAATCCCGTTGTAATAAAATCTCACAAGCTTGATCCACGGTTTCCACCGGATAAATAAAAAATTCCTGATTTTTAACCGCACTTACCACCGCTTCGGACAAGCTTAGCTGATTAAGCACCGCACTCGGAATCATCACGCCTTGTTTGCCC  >  1/1587865‑1588352

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: ‑