Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,486,587	C→T	N172N (AAC→AAT)	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,587	0	C	T	100.0%	102.3 / NA	39	N172N (AAC→AAT)	unknown	Flp pilus assembly protein TadD, contains TPR repeat
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/24);  total (15/24)

TGCTAACAAACAGTGGGCAGTTAAGTGAAAAGGCCAAAGTTCTACAAGCAAGAAATTTAAATCAACTCAAGCGCTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCCAAAAACGGGGAAGTTTATAATCTTCGTGGTGTGACTTATGCCTTAATGGGAAATCCTAATAAAGCCAACGAAGATATAAATAAAGCCCGCGAGTATTTTTTAAACGATGCTGTAGCGGTAAACAATCTTGCTATGCTAAGTATCGTCAACGGAGATTATCGTAATGCCGTTTCATTGTTATTACCTCAGTATTTAAATGGGGTTAGGGAACAGCGTTTAGTACACAATCTAATCTTTGCTTTAATCAAAAATAATGACATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAATACATCTCCTGATGATTTAATTAATGCTCTCAAAAAAA  >  1/1486343‑1486801                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                       tGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatg                                                                                                                                                                                                                  >  1:23537/1‑248 (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: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:151741/251‑1 (MQ=255)                                  ccAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGt                                                                                                                                                                                 <  1:117556/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:27851/251‑4 (MQ=255)                                                         ttAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCattattattgcc                                                                                                                                                          <  1:589520/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:802420/5‑251 (MQ=255)                                                                 acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTAtt                                                                                                                                                  >  1:266487/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)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                       TGCTAACAAACAGTGGGCAGTTAAGTGAAAAGGCCAAAGTTCTACAAGCAAGAAATTTAAATCAACTCAAGCGCTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCCAAAAACGGGGAAGTTTATAATCTTCGTGGTGTGACTTATGCCTTAATGGGAAATCCTAATAAAGCCAACGAAGATATAAATAAAGCCCGCGAGTATTTTTTAAACGATGCTGTAGCGGTAAACAATCTTGCTATGCTAAGTATCGTCAACGGAGATTATCGTAATGCCGTTTCATTGTTATTACCTCAGTATTTAAATGGGGTTAGGGAACAGCGTTTAGTACACAATCTAATCTTTGCTTTAATCAAAAATAATGACATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAATACATCTCCTGATGATTTAATTAATGCTCTCAAAAAAA  >  1/1486343‑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: ‑