Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,141,445	T→C	K65K (AAA→AAG)	unknown ←	Minor pilin of type IV secretion complex, VirB5

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,141,445	0	T	C	100.0%	176.2 / NA	63	K65K (AAA→AAG)	unknown	Minor pilin of type IV secretion complex, VirB5
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/34);  total (29/34)

GTAATTTCACCTTGACGTAAATTCTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTAGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAGCCATTTATTTCACCATACAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCTGTCTGTGAGGCTAATTGGTCGGCATACGCTTTCAGCTCATTTTTATAATGCGAAACTGTATCCGTCCACTGTTTAGCCTCTTTTGCCCATTGTTGCACCGTTTGAAGCATGCTCATTTGACGTTGACCGGCAGCTGCGGCATCAACAACCGGAATACCGCTGGCTTGTACCAATACCGGCATAGATAAAACAGCAATTAAAGTGAGTTTTTTAAACATTTTCATCACGTTCTCCTTATACACAACGTTGTTTATAAAT  >  1/1141218‑1141673                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                     gtgaTTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATgaga                                                                                                                                                                                                               >  1:142724/4‑248 (MQ=255)   gaTTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATgagata                                                                                                                                                                                                             >  1:350373/2‑246 (MQ=255)   gaTTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATgagata                                                                                                                                                                                                             >  1:239945/2‑246 (MQ=255)                                 gTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTg                                                                                                                                                                               >  1:897788/1‑251 (MQ=255)                                          gTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCATTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgt                                                                                                                                                                       <  1:492531/250‑1 (MQ=255)                                              tATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaac                                                                                                                                                                  <  1:808919/251‑4 (MQ=255)                                              tATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaac                                                                                                                                                                  <  1:422696/251‑4 (MQ=255)                                              tATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaac                                                                                                                                                                  <  1:253734/251‑4 (MQ=255)                                                tttttCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTGCATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaacag                                                                                                                                                                <  1:889162/251‑6 (MQ=255)                                                 ttttCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaacagt                                                                                                                                                               >  1:740101/1‑245 (MQ=255)                                                           ggCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGttttgctt                                                                                                                                                                                     >  1:827483/1‑214 (MQ=255)                                                            gCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATg                                                                                                                                                                                                                  <  1:393345/189‑1 (MQ=255)                                                            gCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                                        <  1:85713/247‑5 (MQ=255)                                                            gCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                                        <  1:707896/247‑5 (MQ=255)                                                             cTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacata                                                                                                                                                     >  1:737084/1‑242 (MQ=255)                                                                 tggCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatactgatt                                                                                                                                               <  1:889212/249‑14 (MQ=255)                                                                    cATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGttt                                                                                                                                                             >  1:341961/1‑234 (MQ=255)                                                                          cgcTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                                        <  1:848010/233‑5 (MQ=255)                                                                          cgcTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGa                                                                                                                                      >  1:703525/1‑251 (MQ=255)                                                                                                    aCTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTCGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCgg                                                                                                            <  1:447852/251‑1 (MQ=255)                                                                                                      tAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGaa                                                                                                          <  1:803527/251‑1 (MQ=255)                                                                                                         aTCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATAc                                                                                                       <  1:110855/251‑1 (MQ=255)                                                                                                                    gggTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTga                                                                                            >  1:245680/1‑250 (MQ=255)                                                                                                                                aTTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATcaactac                                                                                                               >  1:620605/1‑217 (MQ=255)                                                                                                                                aTTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCggt                                                                                >  1:687568/1‑250 (MQ=255)                                                                                                                                aTTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCggt                                                                                >  1:244693/1‑250 (MQ=255)                                                                                                                                       accATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGaca                                                                         <  1:25655/250‑3 (MQ=255)                                                                                                                                        ccATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                                        <  1:93089/171‑5 (MQ=255)                                                                                                                                            ttatttCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatactga                                                                                                                                                 <  1:341936/174‑12 (MQ=255)                                                                                                                                              atttCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATcaact                                                                                                                 >  1:695142/1‑203 (MQ=255)                                                                                                                                                    ccATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatac                                                                                                                                                    >  1:707682/1‑155 (MQ=255)                                                                                                                                                     cATATAAACTATGCGCTCCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGcagc                                                                                                                           <  1:435093/187‑4 (MQ=255)                                                                                                                                                      atataaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAg                                                          <  1:490461/247‑1 (MQ=255)                                                                                                                                                      atataaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAg                                                          <  1:260477/247‑1 (MQ=255)                                                                                                                                                       tataaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTgaactaatgc                                                                                    <  1:809715/221‑10 (MQ=255)                                                                                                                                                       tataaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTgaac                                                                                          <  1:809530/215‑4 (MQ=255)                                                                                                                                                       tataaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGt                                                         >  1:735545/4‑251 (MQ=255)                                                                                                                                                        ataaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                                        >  1:783842/3‑151 (MQ=255)                                                                                                                                                          aaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatac                                                                                                                                                    >  1:463491/1‑149 (MQ=255)                                                                                                                                                            aCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatac                                                                                                                                                    <  1:305237/155‑9 (MQ=255)                                                                                                                                                                 gcgcTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTa                                               >  1:655244/1‑251 (MQ=255)                                                                                                                                                                    cTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatactgatt                                                                                                                                               <  1:114637/152‑14 (MQ=255)                                                                                                                                                                         ccAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGAtttt                                       >  1:828561/1‑251 (MQ=255)                                                                                                                                                                           aaTAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaac                                                                                                                                                                  >  1:113527/1‑123 (MQ=255)                                                                                                                                                                           aaTAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                                        >  1:296677/1‑132 (MQ=255)                                                                                                                                                                                 ggCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatac                                                                                                                                                    >  1:66591/1‑126 (MQ=255)                                                                                                                                                                                  gCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTGTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGt                               <  1:592353/250‑1 (MQ=255)                                                                                                                                                                                            cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCtt                        <  1:142079/247‑1 (MQ=255)                                                                                                                                                                                            cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATac                    <  1:691574/251‑1 (MQ=255)                                                                                                                                                                                            cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATac                    <  1:294045/251‑1 (MQ=255)                                                                                                                                                                                            cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATac                    <  1:725014/251‑1 (MQ=255)                                                                                                                                                                                             aCGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATaca                   <  1:337627/251‑1 (MQ=255)                                                                                                                                                                                              cGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATacac                  <  1:58906/251‑1 (MQ=255)                                                                                                                                                                                                aaTTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGacaagac                                                                     >  1:177883/1‑192 (MQ=255)                                                                                                                                                                                                 aTTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACaa                >  1:342920/1‑250 (MQ=255)                                                                                                                                                                                                 aTTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAAc               <  1:799367/251‑1 (MQ=255)                                                                                                                                                                                                   tccggtctgtGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACgt             <  1:364165/247‑1 (MQ=255)                                                                                                                                                                                                      ggtctgtGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCt                                                                                               >  1:463962/2‑166 (MQ=255)                                                                                                                                                                                                      ggtctgtGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGacaagac                                                                     >  1:872881/2‑186 (MQ=255)                                                                                                                                                                                                         ctgtGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATacac                  >  1:336927/1‑240 (MQ=255)                                                                                                                                                                                                              aGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATaaa   <  1:567908/250‑1 (MQ=255)                                                                                                                                                                                                               ggCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTGGCTTGAACTAATGCCGGTATGGACAAGACAGCAATTAAAGTGAGTTTTTTAAACATTTTCATCACGTTCTCCTTATACACAACgt             <  1:679189/239‑1 (MQ=255)                                                                                                                                                                                                                            gCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAAt  >  1:237594/1‑237 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                     GTAATTTCACCTTGACGTAAATTCTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTAGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAGCCATTTATTTCACCATACAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCTGTCTGTGAGGCTAATTGGTCGGCATACGCTTTCAGCTCATTTTTATAATGCGAAACTGTATCCGTCCACTGTTTAGCCTCTTTTGCCCATTGTTGCACCGTTTGAAGCATGCTCATTTGACGTTGACCGGCAGCTGCGGCATCAACAACCGGAATACCGCTGGCTTGTACCAATACCGGCATAGATAAAACAGCAATTAAAGTGAGTTTTTTAAACATTTTCATCACGTTCTCCTTATACACAACGTTGTTTATAAAT  >  1/1141218‑1141673

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