Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,310,568	G→A	E145K (GAA→AAA)	unknown →	Putative transposase InsK for insertion sequence element IS150

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,310,568	0	G	A	100.0%	182.4 / NA	68	E145K (GAA→AAA)	unknown	Putative transposase InsK for insertion sequence element IS150
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/39);  total (29/39)

GCGAAATTATCGCCTATGACTTGTCGCCAAGTCCGAATTTTGAACAGATAACAAGAATGATGGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCGAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTACGCGCAAACCCAGAAAAATTATTCA  >  1/1310323‑1310703                                                                                                                                                                                                                                                      |                                                                                                                                        gCGAAATTATCGCCTATGACTTGTCGCCAAGTCCGAATTTTGAACAGATAACAAGAATGATGGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATg                                                                                                                                     <  1:334562/250‑1 (MQ=255)      aTTATCGCCTATGACTTGTCGCCAAGTCCGAATTTTGAACAGATAACAAGAATGATGGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGttatta                                                                                                                               <  1:565646/251‑1 (MQ=255)       ttATCGCCTATGACTTGTCGCCAAGTCCGAATTTTGAACAGATAACAAGAATGATGGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTtattat                                                                                                                              <  1:693369/251‑1 (MQ=255)          tCGCCTATGACTTGTCGCCAAGTCCGAATTTTGAACAGATAACAAGAATGATGGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGc                                                                                                                           <  1:652036/251‑1 (MQ=255)            gCCTATGACTTGTCGCCAAGTCCGAATTTTGAACAGATAACAAGAATGATGGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCaa                                                                                                                         >  1:224851/1‑251 (MQ=255)             ccTATGACTTGTCGCCAAGTCCGAATTTTGAACAGATAACAAGAATGATGGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGAATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAg                                                                                                                        >  1:104807/1‑251 (MQ=255)                        tCGCCAAGTCCGAATTTTGAACAGATAACAAGAATGATGGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAATGACCAAATGTTATTATGGCAAGCGGTTTGAAAc                                                                                                             <  1:55793/251‑1 (MQ=255)                         cGCCAAGTCCGAATTTTGAACAGATAACAAGAATGAGAGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAAcc                                                                                                            <  1:793152/251‑1 (MQ=255)                                               gATAACAAGAATGATGGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAAc                                                                                                             >  1:478599/1‑228 (MQ=255)                                               gATAACAAGAATGATGGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACg                                                                                       >  1:195624/1‑250 (MQ=255)                                                aTAACAAGAATGATGGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCt                                                                                     <  1:783680/251‑1 (MQ=255)                                                            atgGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGt                                                                         >  1:8208/1‑251 (MQ=255)                                                                    gACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATta                                                                   >  1:746460/1‑249 (MQ=255)                                                                     aCGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATtac                                                                  <  1:756127/249‑1 (MQ=255)                                                                      cGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATtactac                                                               <  1:764772/251‑1 (MQ=255)                                                                      cGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATtactac                                                               <  1:688469/251‑1 (MQ=255)                                                                      cGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATtactac                                                               <  1:199013/251‑1 (MQ=255)                                                                      cGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATtactac                                                               <  1:386807/251‑1 (MQ=255)                                                                        aTGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGaa                                                                                           <  1:39402/221‑1 (MQ=255)                                                                        aTGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACaa                                                             >  1:295809/1‑251 (MQ=255)                                                                        aTGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACaa                                                             >  1:200577/1‑251 (MQ=255)                                                                        aTGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACaa                                                             >  1:60726/1‑251 (MQ=255)                                                                         tGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAAc                                                            <  1:818011/251‑1 (MQ=255)                                                                          gTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCAtt                                                                    <  1:562293/242‑1 (MQ=255)                                                                              aaGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATtac                                                                  >  1:305943/1‑240 (MQ=255)                                                                                 ttGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATg                                                    <  1:209549/251‑1 (MQ=255)                                                                                      cGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTc                                               <  1:760066/251‑1 (MQ=255)                                                                                        gCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGttatta                                                                                                                               >  1:302747/1‑169 (MQ=255)                                                                                        gCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCa                                                          <  1:237340/238‑1 (MQ=255)                                                                                        gCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTc                                               <  1:691596/249‑1 (MQ=255)                                                                                             aCGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAAcc                                                                                                            <  1:772524/183‑1 (MQ=255)                                                                                             aCGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATg                                                    >  1:249462/1‑239 (MQ=255)                                                                                             aCGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTaaa                                        >  1:70947/1‑251 (MQ=255)                                                                                             aCGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTaaa                                        >  1:345302/1‑251 (MQ=255)                                                                                                    ttcttcATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTaaaa                                 >  1:288610/1‑251 (MQ=255)                                                                                                          attcaGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGt                                                                         <  1:5852/205‑1 (MQ=255)                                                                                                               gACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATAc                                                                           <  1:498557/198‑1 (MQ=255)                                                                                                               gACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATtactac                                                               >  1:386324/1‑210 (MQ=255)                                                                                                               gACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGt                                                   >  1:371922/1‑222 (MQ=255)                                                                                                                aCCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagccc                      <  1:653105/250‑6 (MQ=255)                                                                                                                aCCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagcc                       <  1:437419/249‑5 (MQ=255)                                                                                                                aCCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagc                        <  1:873995/248‑4 (MQ=255)                                                                                                                     ggATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGaa                                                                                           >  1:771704/1‑176 (MQ=255)                                                                                                                       aTGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagccctgtggaat              <  1:839662/251‑14 (MQ=255)                                                                                                                        tGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAGTCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGGGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagcactgtggaata             >  1:437838/1‑237 (MQ=255)                                                                                                                          gCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGaa                                                                                           >  1:482179/1‑171 (MQ=255)                                                                                                                          gCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGt                                                   <  1:664557/211‑1 (MQ=255)                                                                                                                          gCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagccctgtggaatac            <  1:509261/250‑16 (MQ=255)                                                                                                                          gCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagccctgtggaatac            <  1:866914/250‑16 (MQ=255)                                                                                                                          gCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagc                        <  1:415198/238‑4 (MQ=255)                                                                                                                              tatCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagccctgtggaatacagaac       <  1:35943/251‑21 (MQ=255)                                                                                                                                   aaTGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagccctgtggaatacagaactcagt  <  1:268034/251‑26 (MQ=255)                                                                                                                                    aTGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATg                                                    <  1:897250/200‑1 (MQ=255)                                                                                                                                          ggTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAg                                >  1:676020/1‑214 (MQ=255)                                                                                                                                           gTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATg                                                    >  1:435297/1‑193 (MQ=255)                                                                                                                                           gTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATg                                                    <  1:658638/193‑1 (MQ=255)                                                                                                                                            ttATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagccctgtggaatacag          >  1:565221/1‑217 (MQ=255)                                                                                                                                                 aGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATtac                                                                  >  1:517322/1‑173 (MQ=255)                                                                                                                                                   cGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagccctgtggaat              <  1:228708/223‑14 (MQ=255)                                                                                                                                                              ggAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagc                        <  1:421185/202‑4 (MQ=255)                                                                                                                                                              aggaaaGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTc                                                                       >  1:564268/4‑155 (MQ=255)                                                                                                                                                               gAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGa                                                                               >  1:609610/1‑146 (MQ=255)                                                                                                                                                                           aGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagccctgtggaatacagaac       <  1:747505/206‑21 (MQ=255)                                                                                                                                                                            ggCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTc                                               <  1:186824/165‑1 (MQ=255)                                                                                                                                                                                  gTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagccc                      <  1:901921/184‑6 (MQ=255)                                                                                                                                                                                          ggAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATtactac                                                               >  1:638374/1‑135 (MQ=255)                                                                                                                                                                                                       gCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTAtt                                                >  1:102279/1‑137 (MQ=255)                                                                                                                                                                                                           tGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCAAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTaagccctgtggaatac            >  1:663244/1‑154 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                        GCGAAATTATCGCCTATGACTTGTCGCCAAGTCCGAATTTTGAACAGATAACAAGAATGATGGAACAGACGATGTTAAAGTTGGACGGCACTACGCCAATTCTTCATTCAGACCAAGGATGGCAATATCAAATGCAAGGTTATCAGCGTTTATTGAAGGAAAAGGCCATTAGGCAAAGTATGTCGGGAAAAGGTAACTGCCTTGATAATCGTGCAATGGAGCGCTTTTTCGGTCGATTAAAGACCGAATGTTATTATGGCAAGCGGTTTGAAACCTTCGCGCAACTTGAGAAAACGCTCCACGAATACGTTCATTACTACAACCATGAATGTATTCAAGTAAAATTAAAAGGACTACGCGCAAACCCAGAAAAATTATTCA  >  1/1310323‑1310703

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