Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,486,551	C→A	A160A (GCC→GCA)	unknown →	Flp pilus assembly protein TadD, contains TPR repeat

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,486,551	0	C	A	100.0%	162.2 / NA	58	A160A (GCC→GCA)	unknown	Flp pilus assembly protein TadD, contains TPR repeat
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/35);  total (23/35)

ATTCATCTCTACTCTATTTACAGCCGTTGCTAACAAACAGTGGGCAGTTAAGTGAAAAGGCCAAAGTTCTACAAGCAAGAAATTTAAATCAACTCAAGCGCTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCCAAAAACGGGGAAGTTTATAATCTTCGTGGTGTGACTTATGCCTTAATGGGAAATCCTAATAAAGCCAACGAAGATATAAATAAAGCCCGCGAGTATTTTTTAAACGATGCTGTAGCGGTAAACAATCTTGCTATGCTAAGTATCGTCAACGGAGATTATCGTAATGCCGTTTCATTGTTATTACCTCAGTATTTAAATGGGGTTAGGGAACAGCGTTTAGTACACAATCTAATCTTTGCTTTAATCAAAAATAATGACATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAATACATCTCCTGATGATTTAATTAA  >  1/1486316‑1486787                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             gttCATCTCTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTa                                                                                                                                                                                                                               >  1:627529/2‑251 (MQ=255)     aTCTCTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTaaatg                                                                                                                                                                                                                           <  1:631311/251‑3 (MQ=255)      tctcTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTaaatgg                                                                                                                                                                                                                          >  1:628643/1‑248 (MQ=255)      tctcTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTaaatga                                                                                                                                                                                                                          >  1:89663/1‑248 (MQ=255)         ctactctactTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTaaatg                                                                                                                                                                                                                           <  1:590627/247‑3 (MQ=255)         ctactctactTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAatgatg                                                                                                                                                                                                                        >  1:671526/1‑250 (MQ=255)           actctaCTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGcag                                                                                                                                                                                                                     <  1:173912/251‑3 (MQ=255)                acttACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGcagtggca                                                                                                                                                                                                                >  1:652437/3‑244 (MQ=255)                 cttACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGcagtggcag                                                                                                                                                                                                               >  1:674347/2‑243 (MQ=255)                        ccGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTaaat                                                                                                                                                                                                          <  1:854547/249‑2 (MQ=255)                        ccGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGCAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTaaat                                                                                                                                                                                                          <  1:784591/249‑2 (MQ=255)                            tGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatg                                                                                                                                                                                                    >  1:23537/1‑248 (MQ=255)                             gCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTaaat                                                                                                                                                                                                          <  1:352841/244‑2 (MQ=255)                              cTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatggc                                                                                                                                                                                                  <  1:620879/251‑6 (MQ=255)                              cTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatggc                                                                                                                                                                                                  <  1:364789/251‑6 (MQ=255)                                aacaaacaGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatggcta                                                                                                                                                                                                <  1:568789/251‑8 (MQ=255)                                aacaaacaGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatggc                                                                                                                                                                                                  >  1:691556/1‑244 (MQ=255)                                         cgggCAGTTAAGTGAAAGGGCCAAAAGTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATggctatgcta                                                                                                                                                                                           <  1:157384/246‑1 (MQ=255)                                           ggCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTAt                                                                                                                                                                                      <  1:826444/250‑1 (MQ=255)                                           ggCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTAt                                                                                                                                                                                      <  1:625293/250‑1 (MQ=255)                                               gTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatggct                                                                                                                                                                                                 >  1:362741/1‑229 (MQ=255)                                               gTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATcatta                                                                                                                                                                                 >  1:535635/1‑247 (MQ=255)                                               gTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATcatta                                                                                                                                                                                 >  1:704594/1‑247 (MQ=255)                                                ttAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTCCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATcatta                                                                                                                                                                                 <  1:18251/250‑5 (MQ=255)                                                           ggCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTa                                                                                                                                                                                       <  1:522092/233‑1 (MQ=255)                                                           ggCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTAt                                                                                                                                                                      <  1:620741/250‑1 (MQ=255)                                                             ccAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGt                                                                                                                                                                   <  1:117556/251‑1 (MQ=255)                                                             ccAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGt                                                                                                                                                                   <  1:151741/251‑1 (MQ=255)                                                                aaattCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAAt                                                                                                                                                                >  1:297629/4‑251 (MQ=255)                                                                  attCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGc                                                                                                                                                              <  1:858571/250‑1 (MQ=255)                                                                  attCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGc                                                                                                                                                              <  1:557515/250‑1 (MQ=255)                                                                                    ttAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCattattattgcc                                                                                                                                            <  1:589520/251‑4 (MQ=255)                                                                                    ttAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCattattattgcc                                                                                                                                            <  1:27851/251‑4 (MQ=255)                                                                                            acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatggcta                                                                                                                                                                                                <  1:270611/187‑8 (MQ=255)                                                                                            acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATggctatgct                                                                                                                                                                                            <  1:585531/191‑1 (MQ=255)                                                                                            acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATtattat                                                                                                                                                <  1:566666/235‑1 (MQ=255)                                                                                            acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTAtt                                                                                                                                    >  1:572086/5‑251 (MQ=255)                                                                                            acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTAtt                                                                                                                                    >  1:266487/5‑251 (MQ=255)                                                                                            acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTAtt                                                                                                                                    >  1:802420/5‑251 (MQ=255)                                                                                             cttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTAttt                                                                                                                                   <  1:829723/248‑1 (MQ=255)                                                                                                aaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTaaa                                                                                                                                <  1:844866/251‑1 (MQ=255)                                                                                                       aTCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGcagtgg                                                                                                                                                                                                                  <  1:718475/162‑6 (MQ=255)                                                                                                       aTCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatat                                                                                                                                                                                                     >  1:122295/1‑173 (MQ=255)                                                                                                       aTCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCAtta                                                                                                                                                    <  1:378780/224‑2 (MQ=255)                                                                                                       aTCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGtt                                                                                                                         >  1:289089/1‑251 (MQ=255)                                                                                                       aTCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCGGCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGATTGTTTCTTAAATGATGCAGTGGCAGTAaataatat                                                                                                                                                                                                     <  1:796776/175‑3 (MQ=255)                                                                                                         cAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGttt                                                                                                                                                         >  1:784040/1‑217 (MQ=255)                                                                                                               gCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCattattattgcc                                                                                                                                            <  1:57228/224‑4 (MQ=255)                                                                                                               gCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGTTAGGGAACa                                                                                                                 >  1:73577/1‑251 (MQ=255)                                                                                                                        gTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAGTGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATc                                                                                                                                                                                     >  1:365875/1‑174 (MQ=255)                                                                                                                        gTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGt                                                                                                                                                                                        <  1:865588/171‑1 (MQ=255)                                                                                                                                                                                                         aaTCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGTTAGGGAACAGCGTTTGATACACAATCTAGTCTTTGCTTTAGTCAAGAATAACGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTaacac                       <  1:276481/251‑4 (MQ=255)                                                                                                                                                                                                          aTCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGTTAGGGAACAGCGTTTGATACACAATCTAGTCTTTGCTTTAGTCAAGAATAACGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAAcaca                      >  1:702045/1‑251 (MQ=255)                                                                                                                                                                                                                   aaGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGTTAGGGAACAGCGTTTGATACACAATCTAGTCTTTGCTTTAGTCAAGAATAACGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAACACATCCCCtgat             <  1:590839/251‑1 (MQ=255)                                                                                                                                                                                                                      cgaaCGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGTTAGGGAACAGCGTTTGATACACAATCTAGTCTTTGCTTTAGTCAAGAATAACGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAACACATCCCCtgatga           <  1:243384/248‑1 (MQ=255)                                                                                                                                                                                                                        aaCGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGTTAGGGAACAGCGTTTGATACACAATCTAGTCTTTGCTTTAGTCAAGAATAACGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAAc                               >  1:664567/1‑228 (MQ=255)                                                                                                                                                                                                                              gATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGTTAGGGAACAGCGTTTGATACACAATCTAGTCTTTGCTTTAGTCAAGAATAACGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAACACATCCCCTGATGATttagttaa  >  1:675143/1‑251 (MQ=255)                                                                                                                                                                                                                              gATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGTTAGGGAACAGCGTTTGATACACAATCTAGTCTTTGCTTTAGTCAAGAATAACGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAACACATCCCCTGATGATttag      <  1:558830/247‑2 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             ATTCATCTCTACTCTATTTACAGCCGTTGCTAACAAACAGTGGGCAGTTAAGTGAAAAGGCCAAAGTTCTACAAGCAAGAAATTTAAATCAACTCAAGCGCTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCCAAAAACGGGGAAGTTTATAATCTTCGTGGTGTGACTTATGCCTTAATGGGAAATCCTAATAAAGCCAACGAAGATATAAATAAAGCCCGCGAGTATTTTTTAAACGATGCTGTAGCGGTAAACAATCTTGCTATGCTAAGTATCGTCAACGGAGATTATCGTAATGCCGTTTCATTGTTATTACCTCAGTATTTAAATGGGGTTAGGGAACAGCGTTTAGTACACAATCTAATCTTTGCTTTAATCAAAAATAATGACATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAATACATCTCCTGATGATTTAATTAA  >  1/1486316‑1486787

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