Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,141,493	C→T	E49E (GAG→GAA)	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,493	0	C	T	100.0%	166.3 / NA	62	E49E (GAG→GAA)	unknown	Minor pilin of type IV secretion complex, VirB5
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/35);  total (27/35)

GTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTAGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAGCCATTTATTTCACCATACAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCTGTCTGTGAGGCTAATTGGTCGGCATACGCTTTCAGCTCATTTTTATAATGCGAAACTGTATCCGTCCACTGTTTAGCCTCTTTTGCCCATTGTTGCACCGTTTGAAGCATGCTCATTTGACGTTGACCGGCAGCTGCGGCATCAACAACCGGAATACCGCTGGCTTGTACCAATACCGGCATAGATAAAACAGCAATTAAAGTGAGTTTTTTAAACATTTTCATCACGTTCTCCTTATACACAACGTTGTTTATAAATCGGAATCCAGTCATTAGCGGCTTCGCCGACTTCTTCGATAATTTGATGTAGGA  >  1/1141250‑1141726                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          gTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTg                                                                                                                                                                                                                                    >  1:897788/1‑251 (MQ=255)          gTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCATTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgt                                                                                                                                                                                                                            <  1:492531/250‑1 (MQ=255)              tATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaac                                                                                                                                                                                                                       <  1:808919/251‑4 (MQ=255)              tATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaac                                                                                                                                                                                                                       <  1:253734/251‑4 (MQ=255)              tATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaac                                                                                                                                                                                                                       <  1:422696/251‑4 (MQ=255)                tttttCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTGCATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaacag                                                                                                                                                                                                                     <  1:889162/251‑6 (MQ=255)                 ttttCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaacagt                                                                                                                                                                                                                    >  1:740101/1‑245 (MQ=255)                            gCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                                                                                             <  1:707896/247‑5 (MQ=255)                            gCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                                                                                             <  1:85713/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:260477/247‑1 (MQ=255)                                                                                                                      atataaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAg                                                                                                               <  1:490461/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:725014/251‑1 (MQ=255)                                                                                                                                                            cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATac                                                                         <  1:294045/251‑1 (MQ=255)                                                                                                                                                            cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATac                                                                         <  1:691574/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)                                                                                                                                                                                                                        gagatacggtttctgTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCaataat             >  1:176362/15‑251 (MQ=255)                                                                                                                                                                                                                             acggtttctgTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGAt        <  1:483430/242‑1 (MQ=255)                                                                                                                                                                                                                                   tctgTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCAATCATTAGCGGCTTCGCCAACTTCTTCAATAATTTGATGTaaga  <  1:51409/248‑4 (MQ=255)                                                                                                                                                                                                                                        ccattGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAACGTTGTTTATAAATTGGAATCCaat                                            >  1:329920/5‑202 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          GTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTAGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAGCCATTTATTTCACCATACAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCTGTCTGTGAGGCTAATTGGTCGGCATACGCTTTCAGCTCATTTTTATAATGCGAAACTGTATCCGTCCACTGTTTAGCCTCTTTTGCCCATTGTTGCACCGTTTGAAGCATGCTCATTTGACGTTGACCGGCAGCTGCGGCATCAACAACCGGAATACCGCTGGCTTGTACCAATACCGGCATAGATAAAACAGCAATTAAAGTGAGTTTTTTAAACATTTTCATCACGTTCTCCTTATACACAACGTTGTTTATAAATCGGAATCCAGTCATTAGCGGCTTCGCCGACTTCTTCGATAATTTGATGTAGGA  >  1/1141250‑1141726

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