Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,486,578	A→G	V169V (GTA→GTG)	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,578	0	A	G	100.0%	144.0 / NA	51	V169V (GTA→GTG)	unknown	Flp pilus assembly protein TadD, contains TPR repeat
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/32);  total (19/32)

CCGTTGCTAACAAACAGTGGGCAGTTAAGTGAAAAGGCCAAAGTTCTACAAGCAAGAAATTTAAATCAACTCAAGCGCTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCCAAAAACGGGGAAGTTTATAATCTTCGTGGTGTGACTTATGCCTTAATGGGAAATCCTAATAAAGCCAACGAAGATATAAATAAAGCCCGCGAGTATTTTTTAAACGATGCTGTAGCGGTAAACAATCTTGCTATGCTAAGTATCGTCAACGGAGATTATCGTAATGCCGTTTCATTGTTATTACCTCAGTATTTAAATGGGGTTAGGGAACAGCGTTTAGTACACAATCTAATCTTTGCTTTAATCAAAAATAATGACATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAATACATCTCCTGATGATTTAATTAATGCTCTCAAAAAAA  >  1/1486339‑1486801                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                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:625293/250‑1 (MQ=255)                    ggCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTAt                                                                                                                                                                                                    <  1:826444/250‑1 (MQ=255)                        gTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatggct                                                                                                                                                                                                               >  1:362741/1‑229 (MQ=255)                        gTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATcatta                                                                                                                                                                                               >  1:704594/1‑247 (MQ=255)                        gTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATcatta                                                                                                                                                                                               >  1:535635/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:557515/250‑1 (MQ=255)                                           attCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGc                                                                                                                                                                            <  1:858571/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:266487/5‑251 (MQ=255)                                                                     acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTAtt                                                                                                                                                  >  1:802420/5‑251 (MQ=255)                                                                     acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTAtt                                                                                                                                                  >  1:572086/5‑251 (MQ=255)                                                                      cttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTAttt                                                                                                                                                 <  1:829723/248‑1 (MQ=255)                                                                         aaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTaaa                                                                                                                                              <  1:844866/251‑1 (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)                                                                                                                                                                                                        atattaacaaagcacgagaatATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGTTAGGGAACAGCGTTTGATACACAATCTAGTCTTTGCTTTAGTCAAGAATAACGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAACACATCCCCTGATGATttagtta                 <  1:888312/229‑1 (MQ=255)                                                                                                                                                                                                                        agaatATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGTTAGGGAACAGCGTTTGATACACAATCTAGTCTTTGCTTTAGTCAAGAATAACGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCaagaaaa  >  1:358053/5‑248 (MQ=255)                                                                                                                                                                                                                         gaatATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGTTAGGGAACAGCGTTTGATACACAATCTAGTCTTTGCTTTAGTCAAGAATAACGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTaacac                                     >  1:33921/4‑209 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                CCGTTGCTAACAAACAGTGGGCAGTTAAGTGAAAAGGCCAAAGTTCTACAAGCAAGAAATTTAAATCAACTCAAGCGCTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCCAAAAACGGGGAAGTTTATAATCTTCGTGGTGTGACTTATGCCTTAATGGGAAATCCTAATAAAGCCAACGAAGATATAAATAAAGCCCGCGAGTATTTTTTAAACGATGCTGTAGCGGTAAACAATCTTGCTATGCTAAGTATCGTCAACGGAGATTATCGTAATGCCGTTTCATTGTTATTACCTCAGTATTTAAATGGGGTTAGGGAACAGCGTTTAGTACACAATCTAATCTTTGCTTTAATCAAAAATAATGACATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAATACATCTCCTGATGATTTAATTAATGCTCTCAAAAAAA  >  1/1486339‑1486801

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