Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	149,939	A→C	I41L (ATT→CTT)	unknown →	Phage protein (ACLAME 1194)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	149,939	0	A	C	100.0%	126.0 / NA	43	I41L (ATT→CTT)	unknown	Phage protein (ACLAME 1194)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (20/23);  total (20/23)

GGACTACCTGCTACTCCCGCCGATGTATCGTATCTAGACAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGACAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTATTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGAACAGTTTGATTTAGATCGGGCGCTGGGAGGCGAAGCTTTTGTGCTCGAAGGCGGAAACA  >  1/149693‑150083                                                                                                                                                                                                                                                       |                                                                                                                                                 ggACTACCTGCTACTCCCGCCGATGTATCGTATCTAGACAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTGCGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTtcttat                                                                                                                                              >  1:445063/1‑251 (MQ=255)      aCCTGCTACTCCCGCCGATGTATCGTATCTAGACAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCa                                                                                                                                         >  1:669634/1‑251 (MQ=255)      aCCTGCTACTCCCGCCGATGTATCGTATCTAGACAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCa                                                                                                                                         >  1:801944/1‑251 (MQ=255)       ccTGCTCCTCCCGCCGATGTATCGTATCTAGACAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCag                                                                                                                                        <  1:98062/251‑1 (MQ=255)        cTGCTACTCCCGCCGATGTATCGTATCTAGACAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCaga                                                                                                                                       >  1:726462/1‑251 (MQ=255)             aCTCCCGCCGATGTATCGTATCTAGACAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGa                                                                                                                                  >  1:682918/1‑251 (MQ=255)                 ccgccgATGTATCGTATCTAGACAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGAt                                                                                                                              <  1:400355/251‑1 (MQ=255)                    ccgATGTATCGTATCTAGACAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCag                                                                                                                                        >  1:508912/1‑238 (MQ=255)                           atcgtatcTAGACAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGc                                                                                                                    <  1:117586/251‑1 (MQ=255)                                  cTAGACAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGa                                                                                                              >  1:736135/1‑250 (MQ=255)                                  cTAGACAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGAc                                                                                                             <  1:530863/251‑1 (MQ=255)                                       cAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGAc                                                                                                         >  1:98372/1‑250 (MQ=255)                                          aTACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAgg                                                                                                     >  1:784652/1‑251 (MQ=255)                                            aCAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCGCCAGAGGTGGCGATTGGTCGTTGCAGTGGACTGACCAGGAc                                                                                                   <  1:107583/251‑1 (MQ=255)                                                      gAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGAc                                                                                                             >  1:473504/1‑231 (MQ=255)                                                      gAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCAtgtg                                                                                         >  1:832155/1‑251 (MQ=255)                                                       aaaTAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATgtgt                                                                                        <  1:403018/251‑1 (MQ=255)                                                          tAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGtgttgt                                                                                     <  1:313394/251‑1 (MQ=255)                                                                acCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTc                                                                                                                         <  1:333373/209‑1 (MQ=255)                                                                acCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGAcc                                                                                                        <  1:708760/226‑1 (MQ=255)                                                                acCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCg                                                                                 >  1:833726/1‑249 (MQ=255)                                                                     aaccaaccCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGAc                                                                                                   >  1:222117/1‑226 (MQ=255)                                                                     aaccaaccCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGAc                                                                                                   <  1:93855/226‑1 (MQ=255)                                                                       ccaaccCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCaa                                                                        <  1:869862/251‑1 (MQ=255)                                                                          accCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGa                                                                      >  1:782192/1‑250 (MQ=255)                                                                           ccccAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACa                                                                    >  1:587067/1‑251 (MQ=255)                                                                           ccccAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACa                                                                    >  1:540849/1‑251 (MQ=255)                                                                                  cccTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGa                                                             >  1:602270/1‑251 (MQ=255)                                                                                    cTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGAc                                                                                                   <  1:700122/211‑1 (MQ=255)                                                                                    cTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGaa                                                            >  1:347519/1‑250 (MQ=255)                                                                                    cTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGAAc                                                           <  1:115129/251‑1 (MQ=255)                                                                                    cTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGAAc                                                           <  1:7423/251‑1 (MQ=255)                                                                                                gCCTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTTGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGt                                                                                                                       <  1:587807/179‑1 (MQ=255)                                                                                                  cTTTTTATTTCCGAAAAAAGGAAAATAAAAAGTTGAATGAAAAAGAAAAAAACCCGTTTTTCAATGGCAATTACGGTATGAGGCAATCCGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCCCCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGAACCGTTTGATTTAGAt                                             <  1:291177/251‑1 (MQ=14)                                                                                                  cTTTTTATTACTCAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGAACCGTTTGATTTAGAt                                             <  1:75066/251‑1 (MQ=255)                                                                                                              cAAAAAAGGAAAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGAAc                                                           >  1:773311/1‑225 (MQ=255)                                                                                                                     ggaaaaTAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGAACCGTTTGATTTAGATCGGGCGCTGAGAGGCGAAg                          >  1:155340/5‑251 (MQ=255)                                                                                                                                 gTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCAcc                                                                                                                                          <  1:837761/127‑1 (MQ=255)                                                                                                                                   tgaatgaaCAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGAACCGTTTGATTTAGATCGGGCGCTGAGAGGCGAAGCTTTtgtactggaa            <  1:83093/251‑8 (MQ=255)                                                                                                                                   tgaatgaaCAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGAACCGTTTGATTTAGATCGGGCGCTGAGAGGCGAAGCTTTtgtactggaa            <  1:401627/251‑8 (MQ=255)                                                                                                                                             aGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGAACCGTTTGAt                                                   <  1:19921/202‑1 (MQ=255)                                                                                                                                             aGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGAACCGTTTGATTTAGAt                                             <  1:248483/208‑1 (MQ=255)                                                                                                                                             aGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGAACCGTTTGATTTAGATCGGGCGCTGAGAGGCGAAGCTTTTGTACTGGAAGGCGGAAACa  >  1:138690/1‑251 (MQ=255)                                                                                                                                                                        ttACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTCTTATCACCAGAGGTGACGa                                                                                                                               <  1:462870/99‑1 (MQ=17)                                                                                                                                                                                                                                                       |                                                                                                                                                 GGACTACCTGCTACTCCCGCCGATGTATCGTATCTAGACAAATACAAGCCAAAGAAATAAACAACCCTAACCAACCCCAAGCCCTCCACGTGAGGGCCTTTTTATTACTCAAAAAAGGACAATAAAATGTTGAATGAACAAGAAAAACACGCGTTTTTAAATGGCAATTACGGTATGACGCAATCGGGGGATAAGTGTAAACTGGTGTTAGTGGGCGAAACTCACCGAGATAATCCGTATTTGGTTATTATCACCAGAGGTGACGATTGGTCGTTGCAGTGGACTGACCAGGACCTATCATGTGTTGTCTCGCTTTGGCAAGACAAGCCCGAACAGTTTGATTTAGATCGGGCGCTGGGAGGCGAAGCTTTTGTGCTCGAAGGCGGAAACA  >  1/149693‑150083

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